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ON THE REALIZATION OF MOTION 
ALONG A GIVEN TRAJECTORY 


E. A. Barbashin 


(Sverdlovsk) 

Translated from Avtomatika i Telemekhanika, Vol. 27, No. 6, 
pp. 681-687, June, 1961 

Original article submitted December 12, 1960 


Methods are given of selecting system parameters which enable the approximate realization of a given 
process, Conditions are derived for which the problem can be reduced to a problem of the theory of 
RMS approximations, 


Let us examine the system of differential equations 


dz, ; 
eet fd ee Tar CYA, (C1, - +5 Cyr Yrre ees Ym t) (i =14, 2,..., a), at 


where Cy,..., Cy, are the constant parameters, and y,,..., Ym are functions of time t, 


Together with the system (1) we assume over the interval 0 =t =T, a trajectory xj =4(t) (i = 1, 2,..., mM) 
which, in general, does not satisfy the system (1). Our task will be to find the parameters c,,..., Cy or the functions 
Y(t), - ++» Yrn(t), Such that the solution x;(t) of system (1), defined by the initial conditions xj(0) = x (the differences 
x} - ¥j(0) are sufficiently small) be an approximate realization of the motion along the given trajectory xj = @ (t) at 
0<t<T. Aswill be shown inthe following, it is advantageous to take as a measure of approximation the root-mean- 
square value of the error at which the given trajectory satisfies the system (1). Thus, we shall select the control param- 
eters Cy,..., Cy (task A) or the control functions y;(t),... + Y¥n(t) (task B), or, finally, both the control parameters 
and control functions (task C), such as to have a minimum RMS error involved. Problems of the aforementioned type 
for certain particular cases have been solved in the works, 


The present work substantiates the proposed procedure for the nonlinear system (1), while in the case where the 
functions  (Cy,.... Ck» Yis--+» Ye t) are linear functions of the parameters c;,..., Cy and functions y;,.... Ym» 
it presents efficient techniques for the selection uf optimum parameter values, 


1. Substantiation of the Method 





Substituting in system (1) the variables 2; = x; ~ p,(t) (i = 1, 2,...,M), we get the sytem 


dz ’ 
oi =A ler+Wilt),--- 5 Za (Es VA (Car - = + + Cyr Yar +s Yur 2) — Ole). 


Introducing the denotations 
ri(t) = 9; (4) — A lVi(t), - -- > Ba» Ob, Zyeay - - -r2ny = 
=f ia+(t),.-.-, 2+ 4%, (8), t}—f,1¥i(),---, B@, t) @=1,2,..., 7), 


we write the new system in the form 


dz 
ay = 2, (ey. oes Zar LV AG (Care ++ Cyr Yareees Ys —7, Ce), (2) 


where it is readily conceived that Z;(0,..., 0,t)=O(i=1, 2,...,M). 








The problem then reduces to the approximate realization of motion along the trajectory z, = 0 (i = 1,2,...,n) 
ato st sT. 


Let us assume that the functions Z;(%,..., 2p), t) are continuous for all arguments, and that they satisfy the 
Lipschitz conditions 


f 2 


|Z, (Zi -- +> Zr t)—Z,(2,..., 2%, IL YS [2,—2]. (3) 


k=1 


il 


We shall not impose upon the function @(c;,..., Ck» Yir--- *Yme ©) any specific limitations, except those 
which are necessary for the existence and uniqueness of the solutions of system (2) in the region D. 


Following the course initially set, we select c,..., Ck, Yx(t).-- + ¥yp{t), such as to have the line z; = 0(i = 1, 2, 
+++,» M) Satisfy the system (2) with a minimum RMS error. This means that the selection of the control parameters or 
control functions should be such as to provide for an approximation of the function rj(t) by the functions ¢(c,... Che 
Yis++++Ymet) at a minimum RMS error in the interval 0 =t =T. 


Let us note that with the transfer from system (1) to system (2), the magnitude of the deviation of the actual 
from the given trajectory is determined by the deviation from zero of the solution z(t) of system (2), defined by the 
initial conditions z,(0) = Zz, where the numbers z are sufficiently close to zero. 


Let us introduce the denotations: 
T 


A? (t) = SIO (Cry +s Cys Yise-+s Ym t)—r, (i), H? = \ A? (t) dt. (4) 
i=1 0 


Theorem 1. Let z(t) be the solution of system (2), defined by the conditions z;(0) = Z; 2:1 <6 (i =1,2,..., 
n). If, at 0 st =T, this solution does not surpass the limits of the D region, then there takes place the estimate 


7 g@LT _.4\"/s 
|z, (t)|<det7 + (4) H (§ = 4,2,..., m). (5) 


The estimate (5) shows that for a sufficiently small H, the trajectory in the interval 0 =t =T can be approxi- 
mated with any degree of accuracy in the sense that the maximum of the deviation of the actual from the desired 
trajectory can be made arbitrarily small. However, if T is sufficiently large, the deviation of the actual from the 
given trajectory may become substantial, since the growth of the maximum in the general case may proceed accord- 
ing to the exponential law. Hence, it is important to define the conditions which provide for the absence of such an 
undesirable growth in error, Let us formulate the said conditions, 


ws ttt 
hk? = sup \ A? (t) dt, 
0<t<T—1 ¢ 
(6) 
where A(t) is determined by the equality (4). 

Together with the system (2) we examine the system 

dz; 

Gr = Zi (tay + + +1 Ens t) (i=14,2,..., 2). (7) 


Theorem 2. If the null solution of system (7) is uniformly asymptotically stable ([4] p. 32), then, for any €> 0, 
the positive numbers n and 6 can be taken such that at h < n, the solution zj(t) of the system (2), defined by the con- 
ditions 2,(0) = 2{,| z{|<5, satisfies the inequalities 

la(t)|<ce (i=1,2,..., x) (8) 
over the entire range of O=t =T (0< T=). 


The theorem 2 is in itself one of the modifications of the theorems of stability at constant disturbances [5]. It 
is clear that in the present case we speak of the disturbances which, in distinction from [5], are limited by RMS values 
rather than by mean values, However, the fulfillment of the inequality 


588 








+1 (+4 tf 

\ A(t)dt< ( \ A?(t)dt ) 

i i (9) 

which results from the Buniakowski -Schwarz inequality, shows that the conditions of the theorem 2 provide for the ful- 
fillment of the conditions of the theorem (1.1) in (5). The validity of the theorem 2 ensues thus from the aforemen- 
tioned theorem. Let us note, besides, that the theorem 2 can be simply proved in a straightforward manner from the 


first scheme in work [6]. This proof is given in the appendix in order to provide a rough estimate for the numbers 5, 
n, which ensure satisfaction of the inequalities (8). 


Then, emphasizing that always H 2=h, it is possible to make the following formulation, 


Theorem 3, If the null solution of the system (7) is uniformly asymptotically stable, then, for any € > 0, the 
positive numbers n and6 can be taken such that at H< n, the solution z;(t) of the system (2), defined by the condi - 
tions (0) = z?,| zi|< 5, satisfies the inequalities | z(t)| < € (i = 1, 2,...,m) over the entire range of 0<t=T 
(0<T =). 


Hence, in the case of uniform asymptotic stability of the null solution of system (7), a sufficiently accurate RMS 
approximation of the rj(t) functions will provide for satisfactory precision in trajectory approximation. 


Finally, we note that if the right sides of the system (7) are independent of t, the asymptotic stability is uniform 
[7]. If the asymptotic stability of the null solution of system (7) is provided by Lyapunov’s function, then this stabi- 
lity is uniform [8]. 


2, Selection of Optimum Parameter Values 





We shall demonstrate on certain particular cases the methods of selecting the control parameters and control 
functions for which the quantity H has its minimum value. 


First, we consider the case where function g has the form 


i= 2 Coxe (2) (i= 1, 2,... mde (10) 


k=1 


The methods of selecting the control functions y,(t) (task B), as well as both the functions y;(t) and the parame- 
ters cj, (task C) for which H has its minimum value, are given in [1]. In the present work, we solve the problem of 
determining the optimum system of Cj,parameters, i.e., task A. 


Obviously, in our case 


Tn ™m 
H? = \ >) (> CikYk (t) — 74 (t)) ‘ae. 


i=) k=1 


A necessary and sufficient condition for obtaining a minimum value of H is to select cj, such that for any i the 
value of 


Hi? = b Cinye (t) — ri w]} dt 


k=1 


be minimum, The validity of this statement follows from the consideration that whatever the selection of the param- 
eters Ciy, ..., Cim, it does not affect the quantity H,,, provided k # i. 


In this way, the problem has been reduced to the simplest problem of the theory of root-mean-square approxi - 
mations, Following the known rule of this theory ([9] p. 396), and assuming the set of functions y,(t),..., Ym(t) to 
be linearly independent in the interval 0 st =T, we find the cj, parameters from n sets of the form 


Ciy (Yj, yj) + eee + Cim (Ym; y;) = (Ti, y;) ( we 1,2,...,8 ). 


jut, &..588 (11) 


T 
Here, we introduce the designation (y,, yj) = \ Yx (t) y; (t) dt. 
0 
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If we consider the vectors Ch Cape ++ +» Cny) (k= 1, 2,..., m), then the n sets of the form (11) can be written in 
the form of a single set as follows 


© (¥1, yj) +... + &m (Ym, yi) =(r, y;) (j = 4, ee m)» (12) 
where (r, y;) represents a vector with the projections 
T 


rie) ys (t)dt (i =1,2,..., A). 


If we find the c;, parameters from the sets (11), then, according to [9], we have 


Hi= r (y,, «0 oe Haas 4) 
. r (y1, i Ym) 





where both the numerator and denominator contain the Gramm determinants of the corresponding sets of functions, 


Hence, 





H? = >) Hi = -— T (y+ + ++ Ys 72)» 
x . r (yi, eees Ym): 2 (v1 y r) (13) 


If the set of functions y,(t),..., Ym(t) is an orthonormalized set in the interval 0 =t =T, then from (11) and 
(12) we get 


Cx=(r yx), H=(rud—D ty P=Dewod—-DWD he 


k=l] i=j i=] k=1 


In the case where the set of functions y,(t),..., Ym(th- ++ is a complete set, H can be made arbitrarily small 
by selecting the number m sufficiently large. 


3. Selection of Optimum Initial Conditions 





Let us now examine the case where the functions ¢, are given by the formulas 


m 
— c t | as ee Sis 
% x wu) «i n) wo 


AS an example for this case may serve the following problem. Let us examine the system of differential equa- 
tions 


a I(x, t)+A(t)z, <F . B(t):, 
(15) 


where X is an &-dimensional and 2 an m-dimensional vector, A(t) is a matrix having n rows and m columns, B(t) is a 
square m-matrix, f(x, t) is a vector with the projections fj(x;,..., Xp, t) (i= 1, 2,...,m). 


Let the vector function x = p(t) be given. To be found is a vector function a(t) defined by the system 


dz 
a =B8B (t) Zz, (16) 
such as to have function x(t) defined by system (15), approximate the function x = ~(t) over a certain segment. 


Let Z(t) be the fundamental matrix of system (16). If at) is the solution of system (16) that is defined by the 
initial condition At) = c, then we have 2t) = At)Z *(to)e. System(15)now reduces to the system 


d 
re = f(x, t)+Y¥(t)e, a7 


where Y(t) = A (t) Z(t) Z~* (to). 
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In this way, the problem evolves down to one of selecting the optimum vector @(¢,..., Cm). If now we write 
the system (17) in the scalar form (1), then for the functions g;, we get the relations (14). 


Thus, let g; be given by the formula (14); to be found are such values of the parameters q,..., Crp at which 
the quantity 


H? = (3 [S cay (t) — re )] dt = [ ey; (0) —r (e)] ae 
0 i= j=1 0 j=t 


assumes its minimum value, Here, through y; is denoted the vector with the projections y4;, Ygj,.-.+ Ynj- According 
to [9] (p. 396), the numbers cy,..., Cc, are to be found from the system 


C1 (Yas Yi) +--+ + Cm(Ym> Yi)=(0, sj) (jf = 1, 2,..., m), (18) 


where the following designations 


T T 
(oyd=\n@ysOadt, Oy) =lrOys Cae. 
are introduced, 0 0 


If Cy,..-+» Cm have been found according to (18), we have 


I ews + os Se 
H? — 1 m 
PY qo + + +> Vou) 





where both the numerator and denominator contain the Gramm determinants of the corresponding sets of vector func- 
tions (or more precisely: the denominator contains the determinant of system (18), while the numerator contains a 
determinant constructed in an analogous manner.) 


If the vectors y; form an orthonormalized set of functions in the interval 0 st =T, i.e., if (Jj. ¥j) = 0 atiwj 
and (yj, yj) = 1, then we have 


T m 
c3=(", yj) (= 4, 2,..., m), H? =\ r*(t) dt — J} et. 
° j=1 
APPENDIX 


Proof of theorem 1, Let x;(t) be a solution of systern (7), 2;(t) be a solution of system (2), where | 2 (te) ~ Xj (to 
<6 ati=1, 2,..0, 2. 





Then it is obvious that 


Lay () — a (18+ \ 1% (4, vee Sy t) — Zi (By + os Sqr OI + 1Q— Dt 


Since, according to (4),| 9; - tj | < A(t), then allowing for the conditions (3), we have 


n 
| 2, (t)— (0) (|< 3+ \ [z> tp (1) — 4 (14.4 (t)] a. 
t, k=1 t 
Introducing the function u(t) = max | z(t) - xi(t)| , we get u(t) d+ {nLu (t) + A (t)} dt, 


thence, by virtue of lemma 1 of [10],it follows that 


t 
u(t) << bette) 4. ett) \ elt) A (4) dt. 
i, 


Applying the Buni akovski-Schwarz inequalities to the integral on the right side, we finally obtain 
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w(t) Sento 4 [ 4 ]" \ ar(at] (19) 
0 


If we take the null solution of the system (7) as the solution x;(t), and assume z= T, ty = 0, then we obtain the 
estimate required. 


Proof of theorem 2, According to the definition of the uniform asymptotic stability for the number € > 0, the 
positive number 5< €, T may be taken such that from | xj(t,)|< 5 (i = 1, 2,..., ) there follows that 





6 
le ()I<% for t>t , lalt+TI<z- (20) 


To may be assumed an integer without impairing entity. 
Let z,(t) be a solution of system (2) defined by the initial conditions i;(t)) = xy(t))(i=1,2,... ,n). 
Att) St st) + Ty, according to (19), there takes place the inequality 


opts "% 
1) —2,()1<4( \) atemat) |, 


where 


In this way, accounting for the conditions of the theorem, we have 
| 4 (t)—2,(t)| << ATER for tot <to+To. 


Let n = $5a"1,2; ifh < n, we have 


') 
14 (t)—2, (I< > for to St [to 4+-To. 


From the inequalities (20) and (21) it follows that the solution z;(t) of system (2), which satisfies the conditions 
| zi(to)| <5, satisfies also the conditions 


le;(t)i<e for to<St<to+To , 12; (to+To)|<4. 


Taking the number t, + Ty as a new initial moment of time, and setting forth analogous considerations, we 
come to the conclusion that at t=t) we have | z(t)|< € (i =1, 2,...,n). Obviously, ty might have been taken 
equal to 0 already at the initial point of the discussion. 


Annotation, Let the null solution of system (7) be stable according to the exponential law, i.e., let all solutions 

of system (17) satisfy the inequalities 
ja (thl< xoBe~*('—t) (i= 1, 2,..., a), 
where x, = max | x;(t»)| , B, and a are positive constants, 
i 

It is readily seen that in proving the theorem 2, it may be assumed that 5 = €/2b and Ty = 1n2B/a, on the 

basis of which A and n are readily computed, 
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THE ANALYTICAL DESIGN OF CONTROL SYSTEMS 


Ya. Kurtsveil’ 


(Prague) 

Translated from Avtomatika i Teletmekhanika, Vol. 27, No. 6, 
pp. 688-695, June, 1961 

Original article submitted January 12, 1961 


A mathematical theory is developed for the derivation of an analytical method of designing controls 


for linear systems with constant coefficients, where this method corresponds to a certain given quadra- 
tic-form optimizing functional, The converse problem is also solved, 


1 The Direct Problem 








We consider a closed control system with generalized coordinates ny, Mg,..., n,, anid control coordinate €, 
described by the equations 


m= >) Syn + m& (i =1,2,..., 2) (1.1) 


j=1 


(bj,j, and m, are constants, Let V(m,..-, Mn» &) = Dyan + & dYam+ak?, 
i,j t 


with aj; = ajj being a given quadratic form. We will solve the tollowing problem (a): for the given initial values 
Ni0e+++* Mp9 We will try to find the functions €(t), n;(t),.... N p(t) in the interval [0, ©), satisfying the equations 
(1.1), and the conditions 


1 (0) = Mo, limm(t)= 0 (=—4,2,..., 2) (1.2) 
and minimizing the integral Th 
J&) =(Vim(,.-- me, & Oat. (1.3) 


0 
(It is assumed that the minimum of the integral )(€ ) is finite.) 


Remark, If V is a positive definite form, then the problem (c) is equivalent to the following problem (c"): 
for given n49,-++» N no to determine the functions £(t), n,(t),...» Nn(t) for t@0, «], such that they satisfy the sys- 
tem (1.1), ng) = ngo (i = 1, 2,..., 0) and such that the integral J(€ ) takes its minimum value, 


eo 
It is obvious that jim inf Eni (t) =0, Since J <o , we alsohave fet) dt<e, 
0 
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f Enj(t) does not tend to zero, then there exists a sequence t,, -+@ with Eni(t,) = =c,c< 0(k=1, 2, 3,...) 
and it wl ows from the formula for the variation of parameters than on some intervals (t, 5, t, + 5), 5> 0 ck =1, 
2, 3, » Eni (t) > c/2. The integral \(€ ) therefore diverges. 


It will sometimes be convenient to use vector notation for the system (1.1) and the form V: 
n= Bn+ mg, 
V (m, &) = (An, n) + & (a, 0) + 8?. (1.1") 
Here (An, 7) is the scalar product of the vectors An and n. 


In the whole of the present article, we will assume that the vectors 


m, Bm,..., B"~*m are linearly independent (1.4) 


and that 


V (mn, ) > BE*, B>O (1.5) 
for any n and €. 


We assume that we already know that for any values of n49,..-» Nno the problem (a) has the solution €(t), 
Ma (t),.-+» Np_(t). Then, as in [1], from the calculus of variations or from the maximum theorem (see [2]), we can 
deduce that there exist functions A,(t),..., Ap(t) such that 


E(t) = > > [mjA; (t) — ain: (¢)], 
i 








(1.6) 
hi = Dain; + a8 — Dj djihi. (1.7) 
If we substitute (1.6) in (1.1) and (1.7), we i 
n= >) (65-7) Nj +e sae Ro, 
j 
fue (0%) w+ 3 (2 mit) 4, ii 


j 


A. M. Letov in [1] carried the reasoning further, using certain assumptions concerning the matrix of the system 
(1.8). We will show that these assumptions can be avoided, starting only from the fact that for any values of ny9,..-., 
Nno 4 Solution of the problem (&) exists. 


Let X be the set of those points (n19,.-++ Mno» A10+---» Ano) Which are initial points for solutions of the system 
(1.8) thattend to zero for t -»«, The set X is evidently a linear subspace in Egn. If the problem (a) has a solution, 
then there exist values of Ayo, ..., Ano, such that the solution &.., nj(t),... Aj(t), ~~~)» 74(0) = nig» Aj(O) = Amo of 
the system (1.8) satisfies the condition nj(t) +0 fort «(i =1, 2,...,m). Since €(t) satisfies 1.6 and Fe %(t)dt <0 
[see (1.5)}, then &(t) > 0°and Emjr j(t)> 0 for t 2 . 


It is not difficult to establish that if a certain component ¢;(t) of the solution of a system of linear differential 
equations with constant coefficients tends to zero for t +, then ¢;(t) 0 fort-+«,** Since joa A; j(t) is a com- 


ponent of such a system, obtained from (1.8) by using a suitable linear transformation of coordinates, then e mjX fo 
~> 0 for t=», We can write the last n equations of (1.8) in the vector form 


A= An— 3, (a, 0) — BY + 3 (m, d) (1.9) 
(B* is the matrix transpose of B). 


*Since see NG(t),-++» Aj(t),... is a Solution of the system (1.8), ee E(t) = tk e nis sin @ jt +5j) + Z(t) where £ 
is real, k=O, €;# 0,04 * 0; fori *j, sind; #0, ifa;=0 andte"z(t) 0 fort +, It follows from this that 
r< 0. 


* "The function ¢;(t) can be of similar form to that of €(t). 
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We now form the scalar product of (1.9) with m 


(m, 2) = (m, An) — >-(m, a) (a, n)— (Bm, 2) + 35 (m, a) (m, 2) 
and so (Bm, A) + 0. 


We conclude in a similar way that (Bm,\)° +0 fort», If is it already known that for some integer $8 we have 
(BSm,) 0 and (Bm, \)’ +0 fort + then when we form the scalar product of (1.9) with B*m we find that (B**4m, 
d) 0, and, as above, we conclude that (B**!m, A)’ +0 fort, It has been proved that (m,A) +0, (Bm, A +0,... 


(B2-!m,) +0 for t ©; and since the vectors m, Bm, ..., B"~!m form a basis in En, then we also have A(t) 0 for 
t > oo. 


We have therefore proved the following: if the functions m(t),..., M(t) of the set (m(t),..., m(thAg(t),..., 
Ap(t)), where this set is a solution of the system (1.8), together with the function &(t), satisfying (1.6) give a solution 
of the problem (a), then (Mo,...+ Tpor-++»Atov---» Ang) EX. Since it has been assumed that for any values My»,.... 
Nno the problem (@) has a solution, then it follows from the above reasoning that dim X = n (where dim X denotes the 
dimension of X). 

Z, Sy 

We write the matrix Q of the system (1.8) in the form Q = (5, _ 2) where S, and S, are symmetric 
ae ag of order n, Z is a matrix of order a and Z' is the transpose of Z, The matrices Q, s transpose Q' and 
# Ss, ~) have the same invariant polynomials, For the matrices Q’ and & So, is 2) this follows from 


the following elementary fact: for a given matrix (gi;) (i, j= 1, 2,...,8) and fixed k and 2 we set hy; = gj;, ifi =k, 
j=t,ori> k,j >Z, and set hij = -gjj if the converse is true; then 


det (ij) = (— 1)'* det (g4;). 


= Z, S; = 
An interchange of rows and columns in the matrix Q* transforms it to the form ( Se, ry this means that 


the matrices ( . de “ ') and ( a > ) have the same invariant polynomials, It follows from this that the 
oe 2 / 


elementary divisors of the matrix Q over the complex-number field, relative to the eigenvalues p = 0, occur in pairs 
(4—p)", (4 +p). 


The space Egy is the direct sum of invariant,cyclic spaces I, (see [(3]), Here the minimal polynomial p, of the 
space I, is either a real elementary divisor of the matrix Q or the product of pairs of complex-adjoint elementary di- 
visors, Let p, be a root of the polynomial p,, The space X is the direct sum of those Is's for which Rep, < 0, Since 
elementary divisors occur in the pairs (1 — p)*, (A + p)*, then to each space I, such that Rep, < 0 corresponds 
a space I,, such that Rep, < 0, Consequently, dim X = n (and there are no eigenvalues for the matrix Q with zero 
m real parts, 


“* Since dim X = n, then for a given (me,.... Mno), ther is a unique set (Ayo,.... Apo), such that (No, ... , Mino» 
Aios-++) Ano) EX. Thus the functions m(t),..., git) from the set (nm, (¢),..., Mn (t), Aa (t),..., An (2), 
- when this set yields a solution of (1.8), together with a function &(t) satisfying (1.6), form a solution of the problem 
e (cx) if and only if (To, -- +» Nnor Atoy- ++» Ano) EX. It thus follows that the solution of the problem («) is 
of determined uniquely. This can also be proved directly, starting from (1.5), 
are The mapping (N0,-+-»Mno)—>(As0,--+ Ano), is evidently linear and ) jo is a definite linear function of 
the variables Mo,...,. Mp9» and so [see (1.6)] — (0) = p> PiNio» where the pj are suitable constants, 
al i 
m- Let E(t), m (t),-.-,Mn(t) be @ solution of the problem (c) for given mo,..., Igo and 7 < 0. We set io = 
t)> ng™), Since the functions § (¢-+ 1), M%,(¢-+),.-.,%n(¢+ 1)  givea solution of the problem («) for given 
Tho» eee No we have 
& (t) = >) pi 15 (*). 
: (1,10) 
We thus arrive at the conclusion that the optimal control is given by formula (1,10), If (1,10) is substituted in 
: (1,1), we obtain the system 
al 


T = D>) (Ou + mips) nj, 


‘ (1,11) 
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which has an asymptotically stable, trivial solution, 


From the linearity of the mapping of (mo,.... Mp9) Onto (Ayo,...+ Ano). We obtain the following result: let 
E(t), m (t),.--, Mn (2) or & (t), 1; atthe Mn (t) be.solutions of the problem (c) for the initial 
values M%o,...+%p9 OF Thos...» Npos then E(¢) + §(¢), mi (¢) + th (t),...,Mn(t) + Mn (2) also is a solution 
of the problem (c) for the initial values 1) + %,.... 1p + p- 


It remains to determine under what conditions the problem («) has a solution, We denote by X,¢ the algebraic 
sum of the root subspaces of the matrix B corresponding to the eigenvalues with zero real parts. 


Theorem 1, Let the quadratic form Wm,..., Ip,6) satisfy (1,5) andlet Wm,..., Ip, &) > Oif 


(t---+Mn)EXo, Sin? +F?>0. (1,12) 


3 


Let functions &, %,...» Mp exist that satisfy (1.1) and(1,.2), Then the problem (a) has a solution. 


Proof, After a suitable change of variable, the system (1.1) can be written in the form 


Y, =By, +18, ye= BeyotleE, ys = Bsys + ls. (1,13) 


Here yy,1; ye. 4: yg. lg are vectors; By, By, By are square matrices having eigenvalues with positive, zero, or 
negative real pari:, respectively. It is obvious that X9 is the set of those points (y;, yg, ys) for which y, = 0 and ys =0, 
Let En, Yin» Yen» Ysn (0 = 1, 2,...) be an optimizing sequence for the problem (a), (En) tends for n + to the 

oO 


e e 


least possible value for J(&), It follows from (1.5) that \ 5a (t)dt ck < oo. We set =, (t) = \ En (t) dt. 
From this and from ° ’ 
t 


tin (t) = €°* (yon (0) + |e 7#* Lib (t) de) = €”* (yin (0) +e 7** Led (1) 


: (1.14) 


it follows that En (t), Yin (t)- are uniformly continuous and bounded functions in every interval [0, T], 
T<, If we extract a suitable sequence, we find that =, (i)—>=(t), yin(t)— yi (4) for n +e 
uniformly in every interval [0, T]. 


From a theorem of F, Riesz, it follows that the derivative + = (t) = & (t) exists almost everywhere and 
foe] 


\e (t)dt <k (see [4], Ch, IX), Taking the limit in (1.14)*, we obtain 
0 


t t 


ys(t) = e#*(y (0) + \ePU:dz (x)) =e (ys (0) + Je? (a) ar). (1,15) 


0 


This means that €(t) and y,(t) satisfy the system (1.13), It is easily verified that 


co oo 


\Y (us 2), yo (O), yo (0), & (t)) dt < tim int |V (yan (0), yan (0), Yan (8), En (2) at. 


0 0 
It remains to prove that y;(t) + 0 for t +, 


since | e®* x, || > 8, |z1| for any 2, (8, >0, t>0), then for y;p(t) we have the representation 


co 


yan (t) = — [oP En (t) dt. 


: (1,16) 


B,(t—t) 
& af 


co 
Taking the limit in (1.16), we obtain y, (t) = —{ L,& (t) dt and so y,;(t) + 0, since 
t 


* We should note that the variation of the functions (1) in any given interval [0, T], T < @, is bounded by a 
constant independent of n. 
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co 


\e (t) dt <k. It follows in a similar way from (1.15) that y¢t) + 0 for t +0, 


0 


Since y,(t) + 0 and ys(t)-+ 0 fort +, then from (1,12) and from the formula for the variation of parameters 
we obtain the results that J(€) diverges if y(t) does not tend to zero, The theorem is thus proved, 


We now write the system (1,1) in vector form: 


1 = By + mk. 


As is known (see [5)), (1.4) guarantees for any initial value np» the existence of functions &(t), and m(t)satisfy- 
ing (1,14), (1.2), and the condition (€) <=, 


Summing up, we obtain the following theorem. 


Theorem 2, Let the conditions (1.4), (1.5), and (1,12) be satisfied, Then the problem (a) has a unique solu- 
tion €(t), m(t),..., ,(t). This solution satisfies the condition (1,10), where the constants pj are independent of 
the initial values Mo,.... Ino 


2. The Converse Problem 





In this section, we obtain a positive answer to the following question, Let the system (1.11), for suitable values 
of pj, possess an asymptotically stable, trivial solution, Then is there a quadratic form V such that (1.10) is a solu- 
tion of the problem (a) ? 


We note that the functions § (¢) = >) Pini (t), M (t),.--,%n(t) —_—‘ form a solution of the problem («) for 


the equation (1.1), the quadratic form V, ahd the initial values Tho.-+++ Mp9 When and only when the functions 
C(t) =0, m (¢),.--.%n(¢) ‘form a solution of the problem («) for the system 


. = b ; j j 
nh Xi ig + mips) Hj + mo (2,1) 


and the quadratic form W (m,,..., na, 6) =V (m, oes 9 Nas >) Pi -+ ¢) for the same initial values. 
j 


We now assume that the trivial solution of the system (1,11) is asymptotically stable. We will use vector no- 
tation, We write the system (2,1) in the form 


n= Dn+mt, D= (bi + mpi), (2.1") 


and the form W in the form W (4, ¢) = (Cy, n) + O(c, n) + 702, where C = (cis) Cj = Cjj. 


We will try to find a positive definite form W, such that £ = 0 is a solution of the problem (a) for (2,1) and 
W, for any initial values, We write the solution of (2.1') in the form 


t 


n (t) =e? nM + hahha mt (t) de. 


0 


(2,2) 


Since |e” |? =>) ((e"");;)"< k,e~",8>0, the substitution of (2.2) in the form W, after some elemen- 
ry 


j 
tary transformations, yields 
foe] 


W (n(t), £(t)) dt =| (Co'no, ey) dt + 


0 


ot? 8 


co 


+4. \ (c + 2{e”*Ce™* dom, eno) (t) dt + 
0 


+ (w ( eP tt—*) m{ (t) dx, t (t)) dt, 
0 0 } me 
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where D’ is the transpose of the matrix D. Since W must be positive definite, then y > 0, (Cy, 7) is a positive defi- 
nite form, and the relation 


c= — 2\e”*Ce?* dom 
0 (2,4) 


is a sufficient condition for ¢ = 0 to be a solution of the problem for (2,1) and W for any Np. But (2.4) is also a neces- 
sary condition, as can be easily seen by substituting ¢¢(t) for ¢(t) and taking the limit for « +0. It is evident that 
W is a positive definite form if (c, n)*< 4y (Cn, n) for n #0, We take any positive definite form (Cn, )and 
determine c from (2.4), Then there exists yg ='0 such that W is a positive definite form when y > yo. When we 
carry out the transition from the form W to the form V, we obtain a positive definite form V such that the solutions 

of the system (1,11) and € (¢) = bs p33 (t) form a solution of the problem (a), In this way we have obviously ob- 


? 
tained all the positive definite forms having the desired properties. 


3. The Case of Bounded Rate of Control 





Let V be a positive definite form in the variables ™,..., Np» &y» &g. We consider the following problem(8): 


for given initial values To,...,+ "no, and &», we will try to find functions &(t), n,(t),..., np(t), such that the con- 
co 


ditions (1,1), (1.2), &(0) = &» are satisfied and the integral /,(§) = \V(n (t),..-, M(é), § (0), : (t)) dt. 


0 
takes its minimum value. If we set Np = & (see [6]) andg = €, we reduce the problem(8) to the problem (a) 
for the system n 
n= yy bisnj + MiMn41, 


j=1 (3,1) 
Thi = 
and the form V (M4, -- +» Mn» Mn+» 6) (see the note in section 1), 
Let (1.4) be satisfied. We denote by B the matrix of the system (3.1) for = 0, B= "1 sy and set ffi = 
0 , 
=| 9 | Then the vectors 
1 
om weer == ome —-. n—t1 
m, Bn=(7), Bin = ("S).--+ B'n = (? ee 


are linearly independent, The following theorem follows from Theorem 3, 


Theorem 4. If (1,4) is satisfied and if V is a positive definite form, then there is a unique solution for the 
problem(s), In addition, the equation 


E =>) Pini t+ Pniré 
j=1 (3,2) 
is satisfied (the constants pj are independent of Mpo,.... no» &o). 


As a corollary, we obtain the results that if (1.4) is satisfied, then there are values of 28) 21,2,...,n+1), 
such that all the eigenvalues of the matrix 


bn, sey ban; Mn (3,3) 
Pis++-> Pn» Pnt+i 
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3) 











have negative real parts. It follows from the results of Section 2 of the present article that if all the eigenvalues of 
the matrix (3.3) have negative real parts, then there exists a positive cefinite form V such that the solutions of the 
equations (1.1) and (3.2) are solutions of the problem (6 ). 
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In order to minimize a quadratic function in the presence of random noise the authors propose a method 
involving extrapolation with respect to three values at equally remote points. The steady-state error 
and the rate of search convergence are computed; the optimal relationships between the search parame- 
ters are found. 


1, Introduction 





In this paper we investigate the operation of an automatic optimizer which searches for the minimum of a 
quadratic function 


y(z)=az*+br+c, a>), (1) 
where a,b,c are the unknown constants when a random error is present at the output of the measuring device, The 
latter statement means that for a specified argument x the measuring device of the optimizer determines the quantity 
¥ (x) = y (x) +2. (2) 
rather than the quantity y. 


As to the error z, we assume merely that it is identically and independently distributed for each measurement 
and that its mathematical expectation Z , its dispersion Dz, and its third and fourth order moments are finite. 


The search method described below for determining the minimum of y differs from the step and gradient meth- 
ods [1,2], Our method uses the idea of extrapolations and is based on the following preliminary concepts, 


If the error z were to be absent, then the abscissa of the minimum for (1) 


a b 
min = — 3 








could be found by extrapolation with respect to three values y(x) at arbitrary equally remote points xX9—h, Xo, X» +h 
(the number h shall be called the trial step). 


In fact, using the substitutions 
Yo = Y (20) = ary + bry + ¢, 0 =Y(Xoth) =a(xo+h)? + b(%+h) +e 
and specifying the working step by the formula 


h > — ye 





A =—=-5 ’ 
on 2 ys yt — ye 
it is easy to prove that 
fot b 
y=H+h=—sz 


It is natural to attempt the application of formulas of the same type to the function (2), but in that case the 
quantity x, will be random; thus not only it but its mathematical expectation x, may differ from x,,;,- However, 
for the natural assumptions made relative to the error z we should expect that X; will lie closer to x in than Xo; if 
this is so, then it is logical to repeat the process with the same trial step beginning at x;, etc. 


m 


These concepts lead to the following preliminary search scheme. 


The point xX» and the trial step h are chosen arbitrarily, and the transition from xp to xp. (which we shall call 
one cycle) is achieved as follows: we measure the values 


Yn =Y(an—h), Yn=Y(an) Yet=Y(ta+A), 
and then compute the working step according to the formula 


— — - se ee 
“2 yr4yt—2y os) 


n 





and assume 


Inti = In + An. 


We find, however, that in this preliminary form the search algorithm is practically useless for the following 
reason, 


The numerator and denominator of the right side of (3) are computed with reference to the function (1) and 
yield the respective results 


Yn —yt = —Aahz,—2bh, yn + yt —2yn = 2ah?, (4) 


whence it is evident that when there is no error z the denominator is always positive. In the same formula (3) the 
denominator may prove to be negative or equal to zero due to the random nature of the quantities A this leads to 
false or impossible values of the working step, Evenif we discard the latter values by agreeing to omit one cycle in 
this case, the search remains supersaturated with random false and excessively large working steps ("hunting" occurs) 
and computations will show that the process diverges, 


In order to avoid this situation, we can choose a certain precautionary number k > 0 in advance and agree to 
perform the working step (3) only in the case where either the denominator of (3) or its absolute value prove to be 
greater than k, In the contrary case we also have a choice: we either omit the cycle, or perform the working step 
(3) while replacing the denominator by k, In all there can be four variants of the search scheme, 


Investigations similar to those performed below demonstrate that the optimum results are obtained for the fol- 
lowing variant: 


y,-—Yf 
tnt = Int A, if Yo+Yt—2Y,<k, (5) 








all 








— h ro -V¥z : ns + ‘ 
tnt tnt ee Ve bE OY, ) 








n 


It is this variant of the search algorithms which we shall investigate. It is clear that this variant consists of a 
combination of the extrapolation and gradient methods. The random values of the argument x, X,,... obtained by 
means of (5) cause the function yto have the random values y,, yz,.... Moreover, within each cycle we encounter 
the additional values 


ya = y (anh). 


The exact value of the minimum for the function (1) is ymin = ¢-(b*/ 4a), and the purpose of the optimization 
consists of approximating the quantities y, = Yn - Ymin» Yh = Yn ~ Ymin *© Zero in the probability sense as n increases, 
Since y,. y,* = 0, an appropriate measure of their over-all deviation from zero, for example, if the arithmetic mean 


1 a , , 
= -_ (Yn 4 Yn + Yn) 


We shall study the following probability problems associated with the series {u,} = Uy, Ugs-s.- 
a) computing the mathematic expectation U, and the dispersion Dun; 


b) finding the convergence conditions for the series {Up} and {Du,} and computing their limits U = lim Up 
and D = lim Du,;: 


c) determining the convergence rates for the series; 


d) finding the relationships between the characteristics enumerated here for the search quality and their depen- 
dence on our arbitrarily chosen parameters— the trial step h and the precautionary number k; 


e) determining the allowed requirements relative to the optimalness of the choice of parameters and determin- 
ing the possibility of realizing these requirements. 


The quantity D is called the steady-state dispersion, and the quantities U and Us, = lim Yn are called the 
steady-state search errors. In fact, we are interested only in U, but many details of the investigations are associated 
with Us». 


The series {Up} and {Duy} prove to be convergent at the rate of geometric progressions that have the posi- 
tive denominators A and L, respectively, The basic part is played by A which shall be called the convergence coef- 
ficient, Its closeness to zero denotes rapid convergence, and its closeness to unity denotes slow convergence. 


We shall enumerate the final results. 


1, The convergence region for the process is shaded in Fig, 1 and does 
not contain a series of values for k that are close to zero, Thus the precau~ 
tionary number cannot be made arbitrarily small (at least not when using the 
trial step h). 


2. A simultaneous reduction of the convergence coefficient A and the 
steady-state error U is impossible: if h and k are chosen so that A - 0, then 
O--e; if U0, then A*1, 











0 h 


3. In the convergence region there exists an optimal curve which has 
Fig. 1. The convergence region the property that from any other point in the region it is possible to move to- 
and the optimal curve: 1) optimal ward the curve in such a way that A decreases and U remains the same, or so 
curve; 2) k = 2ah" that U decreases and A remains the same. 


For movement along the optimal curve in the direction of an increase 
in the trial step h we have A + 0, but U-»«; in moving in the direction of a decrease in h we have U + 0 but A +1, 
The steady-state dispersion D tends towards zero both for h + 0 and h +e, 


The optimal curve has the shape shown in Fig. 1. Thus, on the optimal curve there exists a positive minimum 
for k, and k increases without limit both for h + 0 and h + *, In other words, it is not practical to perform this pro~ 
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cess with a precautionary number which is too small, and when resorting to large or small trial steps (for accelerating 
or, on the contrary, for refining the search process) we should in both cases use sufficiently large precautionary num- 
bers. 


All the auxiliary computations are cited at the end of the paper in Appendix I and Appendix II. 


2. Basic Definitions 





Instead of the arbitrary parameters h > 0 and k > 0 we shall henceforth use the following quantities* that are 
associated with them in a single-valued manner: 


| k 
sad 2aht and B= 2ah* 1 
(6) 


with a range of variation 
a> 0, B>—1. (7) 


We shall reduce the algorithm (5) to a convenient form of notation. Assume Zp, Z;* are the values of the error 
for the measurement of y at the points xp, x, +h, respectively, so that 


Ya=Yatta Ya=ya tse. 


We shall assume that 
, b ’ ' b* 
= t—fmin= tts, ¥ =¥—Ymin=y—(¢—z) 
and 
Sn = O(n + tn — 2p), t=V £er—2). (8) 


When these definitions are used we obtain 


y’ = az (9) 
instead of (1) and algorithm (5) becomes 
. B , th 
tH = TpA*t+ TTF for S,<B, 


> s > t (10) 
m= TE7 + TT; 8 m>B 
i.e., after obvious transformations which are made on the basis of (4). 


Since the previous values of x and y are not encountered again, we shall omit the primes to simplify the nota- 
tion (i.e,, we shall assume x‘ = x, y' = y. 


Assume the function g{s) is defined for all s = 8. We shall introduce the definition 


g(B) for s<B, 


[g(s)lp ot ae)= {FE Hg gO a1) 


Relationships (10) can now be written in simpler form: 


*We should remember that a is an unknown, but fully defined positive constant. 











Sy by 
fayi = F + el Zn + [i Flo . (12) 


Assuming that the density distribution for the error z is specified, it is easy to find (cf, the computations in 
Appendix I) the density P(S,T) for the joint distribution of the random quantities 


Sq = tn + 2¢ —2zp, T, = tn —tay (13) 


which, in view of the propositions governing z,is independent of n; then we can also find the density p(s,t) of the 
joint distribution for sy and ty. Here P and p are related by a simple function (cf, Appendix 1), In particular, if we 
use the definition 


co oo 


\ 1™P(S,T)dT = Pm (S), \ t™ (s, t)dt = pm (8), (14) 
—o0 -—oo 
then 
| 1 
| Pe)= Pelz), Prl)=arPa(z)s al) = ear Pa(Z)> (15) 
these functions are nonnegative and integrable over the entire axis -« < s < +, and 
Pi (Ss) = ps (s) =0. (16) 
It is evident that pgs) is the density for s,, whence 
J pos ds = 1. an 


Then,we shall define the density of the random quantity xp in terms of f(x). From this definition and from 
(9) it follows that 


foe] co co 


| fn(z)dz=1, a | a%q(z)dz= ym, a | af (z) dz = yf, = 


—os —oo —oo 


Since the errors z,, 2 arise for measurements performed after the value of xp has already been found and by 
definition do not depend on the abscissas of the measurements points, it follows that the quantities (8) and x, are in- 
dependent (whereas Xn4,, of course, depends on Xp, sp andt,), Therefore the density of the triplet (xp,sp, tp) will be 


/n (x) p(s, t). (19) 


From algorithm (10) it is required to derive the mathematical expectations Up and the dispersions Du,. On the 
basis of (9) we have 


an ae on 2 
Yu = GIy, Yn = a(z, +h)’, (20) 


whence 


‘a 2ah* i 
ln = (Ya + Yn + yn) =Ynt jz = Ya t+ 3° 


This simple relationship permits us in many respects to reduce the study of the series {up} to a study of the 
series {y,}. 


In particular, from this we have 



















- 1 
in =Ynt ay? Duy = Dyn, 


and therefore 


U=U.+4, 
and 


D = lim Du, = lim Dy,. 
n--co n--0o 


With respect to computations of the dispersions, we note also that 
Dyn _ yn = (Yn)?, 
whence 
D= Vi a Us, 


where 
n-+0o 


Thus, the problem reduces to computing the quantities yp and ¥;, from formula (20) and the algorithm (10) or 

(12); then 
U,=limg,, V3=limy?. 
n->0o n--co 
and finally 
| 
U=U,+%) D=Vi— Ui. 
All of the enumerated quantities will evidently depend on the arbitrarily chosen trial step h and precautionary 


number k; i.e., in terms of the definitions (6) they will depend on the arbitrary parameters @ and 6. A study of this 
relationship is a significant portion of our investigation. 


3. The Recursion Relationships for the Mathematical Expectations and Dispersions of yp 





From (12) and (20) it follows that 


s? 5, 
von ellasep],2+2[ aie | oot er] = eee on 
Therefore, in view of (19) 


co 


In = ( \ ( W (a, s, t) fn (x) p(s, t)dxdsdt. 


Performing the integration while taking (15), (16), (18), and (21) into account, we obtain the recursion relation- 
ship for the mathematical expectations y,, 


Ynti — Aj, + B, 


(22) 








on- 











Analogously, we obtain the relationship for Yn 





yas = Ly + Mj, + N, 








(25) 

where 
i= | [arta], mtoee os 
M = 6a \ [aar |, Ps (9) as, (27) 


The recursion relationsjips or finite-difference equations (22) and (25) make it possible to compute the mathe- 
matical expectations y, and Yn» and then their limits for n +, In fact, solving Eq, (22) in the conventional manner, 
we find: 


z B _ B 
Yn = 5,A" + 7—G (6, = J.—;—z): if A+i, 


Jn = J, + mB, of Aw, 





(29) 


From this it follows from the positive nature of A and B [cf, (23), (24) ] that the steady-state error 
U,=limg, 
n->oo 
exists and is finite only when the condition 


aA< tf, (30) 
is satisfied; in this case 


B 


U,=;— > 0. (31) 





(The special case 0» = 0 can easily be proven to lead to the same results), 
We shall proceed to relationship (25), 
We shall substitute the solution (29) which we have found within it while making use of (31); as a result we ob- 


tain 


yaaa = Ly + Mo,A* + MU, +N. 








Now it is also possible to solve this equation: 








ya = tol” + Moy 5-4 if Leet, Le, 
yh = tol” + Mon” + “20d 8 iff L=A<A, 
yi = y+ Mo. ~—*+ +n(MU,+N), i < Sind 


(the case L = A>1 is not of interest to us), The initial value 7) can easily be determined, We can again prove that 
the condition bounding* the series { yf} is the condition 


L<i, (32) 


which also automatically assures its convergence; in this case we have 


Ve = limy, = SPet* (33) 
n-+00 od 
and therefore the steady-state dispersion is written as 
MU,+N 
D=Vi—U3 = — ts _ vi. (34) 


All of the quantities determined here are functions of the arbitrary parameters @ and 8 ; this is directly evi- 
dent from formulas (23)428) if we take (15) into account. Therefore we shall write A = A (a, ), U = U(a, 8) ete, 


4, The Allowed and Optimal Values of the Arbitrary Parameters 





All of the computations and the missing proofs for this section are cited in Appendix IL 


1, The region of allowed parameter values or the region of convergence for the process [satisfaction of ine- 
qualities (30) and (32)] is bounded only from below by the curve 8 = 8(&) lying in the band 


where 6(a) is a nondecreasing function of « and does not coincide identically with 
—1 (Fig. 2). In view of definition (6) this is result 1 (cf, section 1); Fig, 2 becomes 
Fig, 1, 


2. The level lines for the function A(a, 8) are arranged in the manner shown 
in Fig. 3, In particular 





a) small values of A are possible only for small @ (i.e., from A — 0 it follows 
that a-— 0, and if 6> 0 then it is also true that 6 +0: 


b) moreover, if 6 +, then A(a,8)- 1 uniformly with respect to a, 
The level lines for U(a,8) are arranged as shown in Fig, 4. In particular: 
a) if a — 0, then Ula, 8) uniformly with respect to 6; 


b) small values of U(a,8) are possible only for large ; i.e., from U-+0 it follows that a+ ., and even B+« 
and 6/a + «, 


c) each level line U extends at both end in the direction of an unlimited increase of 8. 


The enumerated facts make it possible to draw a number of conclusions. First of all, result 2 in section 1 de- 
rives from them, In fact, from A — 0 it follows thata -— 0 and therefore U + e; however, if U + 0, then 8 — © and 


therefore A+ 1, 


* For the dispersions, it would have been permissible not to require any more than satisfactory search when the other 
conditions were satisfied. 
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Thus, a simultaneous decrease of A and U is impossible. Therefore, the problem of involving the optimal 
choice of the parameters « and 6 (i.e,, the problem involving the optimal values of A and U) can be posed only in 
the sense of finding the relative minimums of one of these functions for fixed values of the other. 


3. When the problem is posed in this manner it cannot be solved, and therefore we come to the concept of an 
ptimal curve in the sense of result 3 in section 1, 


fA! 


0- 




















ST 


Fig. 3. Fig. 4. 


In order to prove this we employ the following reasoning. 


Irrespective of which of the two functions A and U is assumed fixed, the relative extremum for the other func- 
tion is found from the equation 


aA @A 
“0a + =OB 
au eu |=" (35) 
“Oa OB | 


which in expanded form (cf, Appendix II) is written as 


B © B ph 
a \ pa(s) ds | Ps =B | polsyds | —B+-a\ wear? (L) as]. (a6) 





This equation defines a curve on the plane (a,8) which a priori can degenerate into an empty set or can decay 
into several branches (into the line for relative minimums, relative maximums, etc,), However, none of these 
branches intersects the a axis, since for 6 = 0 the right side of (36) goes to zero and the left side is positive, 


In geometric language, Eq, (35) expresses the condition governing the coincidence of the tangents to the level 
lines for A and U at their common point 


dA /0A au. au 


~ Oa / BB — Oa / OB" 


Therefore, Eq. (35) can be solved graphically by superimposing Fig, 3 onto Fig, 4 and then taking the “line of 
contacts" for both families, From Fig. 5, it is not difficult to foresee that such a line exists and is precisely the line 
of relative minimums, Moreover, if the line of contacts were to split up into several branches, then one of them 
would also be the line of minimums, We shall provide a more thorough justification for this convenient reasoning, 


From (35) it is evident that both functions participate equally in the graphical process, We shall take an arbi- 
trary level line Ula, 8)=C. Both of its ends extend in the direction of an unlimited increase of 6 under conditions 








where A - 1, and since A < 1 in general in the region under study, it follows that on the line U = C the function A 
has a minimum point, For continuous variation of C this point describes the curve shown in Fig, 5 (or coincides with 
one of the branches shown in Fig. 5 in the case where it splits up); this point is optimal in the sense of result 3 in 
section 1, In fact, by the very graphical process which we have used it is possible to move onto the curve along the 
level line U from any other point in the region under study and thus decrease A while retaining U constant. 


In view of the equal weight given to A and U in this process, 
they can be interchanged in the reasoning above. On this basis, it 
follows in particular that points on the optimal curve exist not only 
on each level line U but also on every level line A, In other words, 
the optimal curve intersects all the lines of both families, Here 
small values of U cause it to go in the directions © +6 ,8-—+« and 
even 6/a +; small values of A cause it to go in the directions of 
a-»0 and 6-— 0 since it cannot intersect the « axis and remains 
above that axis, 


Thus, the shape of the optimal curve actually corresponds to 
that shown in Fig. 5. 





Fig. 5. 1) Optimal curves Assume the equation for the optimal curve is 8 = &), It has 
just been demonstrated that (a) > 0 and Ha)/a-—+« together 
with a, Making use of formula (6), we obtain the equation for the 
optimal curves in the form 


k= TOO) or = k= 2ah? [1+ O55 5)]. 


Now it is clear that k > 0, and that for h + 0 or #(i.e., for a «© or 0) we have k — «; this is the last postulate 
in the result 3 of section 1 [relationships (6) cause Fig. 5 to make the transition into Fig. 1} 


APPENDIX I 


We shall study the random quantities (13), In view of the proposition cited in section 1 relative to the error z, 


the quantities Zn Zao n=1,2,... are independent and have the same density distribution ¢€ ); here 


co 


| c"e@at<too (m=1. 2, 3, 4). (37) 


—oo 


Assume Q(S,T) is the joint distribution for Sp and T,; i.e., it is the probability for the simultaneous satisfac- 
tion of the inequalities 


a t+et—2s=S <8, 2 —zt=T,_<T. 


In that case 


Q(S, 7) = i © (0) 9 (Gs) 0 (Ee) dE at a, ori 
3 


where the region G of the space (¢ , £1, €2) is described by the inequalities 
Sitle—2<8, GA—&<T. 
For the substitution of variables 


» 
2 


t=n- 











as 


ate 











with the Jacobian 


9G, 0a, Ge) _ A 
0 (n, , N2) 2 
these inequalities become 1, <T, "Nz < S. 
Thus, (38) can be transformed to 
‘ T Ss co 
n n 
QS, T) => ( dm \ dM ( p (n—*)¢ (n+-t)e (n—-t) an. 

—co —0o —oo 


From this we find the density distribution for the pair (S,, Tp) 


ee: ee 


(39) 
We note that this function is even with respect to T, and therefore‘ the functions 
co 
P_ (S) = \ T™P(S,T)aT 


are identically zero for odd m. However, if m = 2,4, then P,,,(S) are integrable over the entire axis since the inte- 
— = 
grals ( P,,(S)dS can easily be reduced to linear combinations of the integrals (37) using (39) and (40), 
Having determined the distribution function Q(S, T) for (13) it is easy to find it for (8) as well, It is obvious 


that 
q(s,t) =P {s, <s, <= P{as, <s, ye r.<th=o(+, V =.) : 


where P is the probability symbol, 
In this case the density for the pair (sp, tp) is 








yet (EM, (s 9/ By) 
Ps. = 359 =aSeT a @ ~~ uva\a' a‘): 


This relationship between p(s, t) and AS, T) leads in obvious fashion to (15) and to (16), since we already know 
that P,(S) = Ps(S) = 0. Note that if g{¢ ) = 0 outside [a, b] then for P,,,(S) such a segment will be represented by 
[-c, c], c = Ab-a), This can easily be derived from (39) and (40), 


APPENDIX II 


We must now perform a detailed investigation of the functions (23)(28), We shall begin with certain general 
propositions, 


Assume that on the semiaxis, (0, ©],the function g(s) is continuous and has the limit ge), and that P(S)=0 is 
integrable; assume further that p(s) = @/o)P(s/a), Then 


co co co co 
tim { ¢ (*) p(s) ds = g(0){ P(S)ds, lim | g(s) p(e)ds = g (00) P(S) aS. 
0 0 re 0 

















Proof: We have 


g(s) p(s) ds = \ g (aS) P(S)dS, 
0 


oc" 8 


where in view of the boundedness of g(s) the latter integral converges uniformly relative to 0< a< e; from this the 
validity of both transitions in the limit follows. 


Corollary. If the same conditions are satisfied over the entire axis, then 


lim { g(s) p(*)ds = (u+-v) 6 (0), 
e+e <, 
Pow \ g (8) p(s) ds = jag (— 00) + vg (00), 


where 
0 co 
a y P(S) dS, v= ( Pisyas. 
—co 0 


Assume g(s) is piecewise differentiable and that sg‘(s) is bounded; assume further that 


co 


R (a) = { g(s)p(s)ds, Ri(a)= | 6) P(£) ae 


—oo 


Then 


4 co 
R’ (a) =—- J 86’ (6) p(s) ds, 
—oo 


(41) 
, i ° » 
R, (a) = [2 (a) + { sg )P(£) as}. a 
— 4 
Proof. We have 
co 
R (a) = ( g (aS) P (S) d5. 
—oo 
Formal differentiation yields 
co 
R’ (a) = { Sq’ (aS) P (S) ds. ae 
—oo 


Since the function Sg(aS) =, sg'(s) is uniformly bounded for the variation of a over any finite segment 0< o% 
S&S Qg, it follows that the latter {te gral converges uniformly over this segment; therefore, formal differentiation 
is valid for any & >0, In order to derive (42) it is sufficient to note that Ry = oR, 


Note that in definition (11) we have 


Sg (—0)=g(B), gg (00) = g (00) 
and 





ds = 


agp (s) + for s < B, 
g’ (s) for s>B. 

















and (26) we find 





= S { Po (8) de. 


directly. 


In order to compute the derivatives with respect to a, we make use of (41) and (42) while taking the defini- 


tions (15) into account, As the results we obtain 





co 
4 2s 
ec. (1 + 2) Po (s) ds, 
~ 


1, 
a] 


‘ oo 
Ly fy fi ay Po) ds 


From this it follows in particular that 
A,>0 , L,>0, 


As ° Le have the sign of B. 


The computation of the derivatives U = +s + = is now obvious. 


We establish the fact that 


A(a,B) << VL(G, 8) » 








Be |, =m 2 <2, 8) <limL = 1—p [4 — 
0 a—+0 a-0o 


where definition (11) is at the left and 


a a 
OF Bp <7 P< oe oy Bp 
a a 
Hs 2B +1 < U0 (a, B)< oe a7 


are at the right. 


= Va tpa(S ) a), 


mim 4<A(@, 8) <lima=1—p[y_7 Pe 
0 a0 a--oo 


In the latter expressions the left sides should be replaced by zeros when 6 < 0, 


We shall now compute the partial derivatives A, B, L with respect to8 and with respect to a, From (23), (24), 


(44) 


(45) 


(46) 


(47) 


(48) 


(49) 























We also establish the fact that 


psa for OB<1, 
M (a, 8B) = 2 
fase Mrs ae for B>0; (50) 
a? 
N (a, B) << ps a+py° (51) 


In formulas (48)-(51), u; = const > 0, 


Proof, The expressions for the functions are taken from (23)-(28) while taking (15) into account, and as a rule 
we shall use definition (11), Note that (23)-(28) are nonnegative. In order to prove (45) it is sufficient to make use 
of (17) and the Buniakowski inequality 


A= \ i ral V pols) Vv(s)as< 


—oo 


i 
co ” co _ 
<V j [rial noe V po(s)ds = VL. 


—oo 














? 
Further, if g(s) = Ges then it is obvious that 


p? 


0 for B<0, 
=r 1,- | 


p? 
a+ for B>0, 


and Bt 
&,(— ©) = TBP 8p (0c) = 1. 


From this we obtain: 


’ * Rt us B2 
lim A(a, 8) =| 7s) im A.B) = a +¥. 


Since in this case we have 


co co 
pt+v= \ Po(S)dS = \ po(s) ds = 4, 
—oo —co 


in view of (15), (17) and A(a, 6) increases with respect to a (44), we arrive at (46), In completely analogous 
fashion, we prove (47), The estimates (48) are obtained directly from (24), It is obvious that 





B B co co 
: ds<B a yes Fae —— . 
ater \” <?< aap FP eae) wer ) 


B 0 
since j pads > \ pads forB2=0,and Pa(s)= ae Ps (+) (cf, 15), we obtain the relationship (48) where 
—oo —oo 0 fo) 
m= | Pa(s)as, pa = | Pa (S) ds. 


—oo 0 
However, if 8< 0, then we can definitely substitute zero for the lower boundary of B. Now (49) derives from 
(46) and (48) since U= B/(1-A). Inequalities (50) and (51) are proved by the same method as that used to prove(48), 
but for the function #/(1+s)* we use the coarse estimate s*/(1+s)*< 1 for 0 = s < 1 and the more exact estimate 
#/(1+s¥'s 87/(1+8¥ for s= 8 =1. 





le 


48), 











We shall now study the level lines A and L, From (44) it is 
evident that for any level line AQ, 8)=C we have 


A d <0 for B>0, 


—-+" 3 = 00 for 8 =0, 


>0 for g<0. 


From this it follows that any level line or portion of the level 
line lying above the a axis is a decreasing function of a, and any 
level line or portion of a level line lying below the @ axis is an in- 
creasing function of o; if the line lies on the axis proper or intersects 
Fig. 6, it, it has a vertical tangent, 





The function A(a, 0), in accordance with (44) and (46), in- 
creases with respect toafrom 0 to l-y. Therefore for 0< C < 1-y the level line A(a,8) = C has one and only one 
common point a. with the o axis; here ac-> 0 together with C, However, if 1-y = C < 1, then no such point exists. 
The level line adherestothe 6 axis at the point 8 cwhich is determined from the equation lease 8° |- 

(i + 8)? 
[cf (46)]; from this it is evident that 6; also tends toward zero together with C, 
The line A(a,8) is in a special position. 


In accordance with (46) it has only one branch lying below the straight line 6 = -4, This fully determines 
the disposition of the curves in Fig. 3, 


If Ry(S) = 0 outside [-c, c] (cf. Appendix I), then the derivatives (44) may go to zero (for |6 |> ca, This has 
a negligible effect on the remaining reasoning and leads to Fig, 6 instead of Fig. 3, 


The investigation of L(a,8) is also based on (44) and on the estimate (47) that is fully anologous to (46); there- 
fore it yields the same results, But due to inequality (45) the region L< 1 is contained within the region A < 1, and 
this leads to statement 1 in section 4 where the curve a =8(@) is the level line L{a,8) = 1 


1 
We shall study the level lines U, From (49) it follows that for U = Up + 3a We have 
Mi to <U Bm +e, 
2B + 7 2B+4 + i 


where we should substitute zero for ps in the left side when 6< 0, 
Thus, the level line U(a,8) lies between the hyperbolas 


a he a oe 
Moma tae" Moat tag 





(52) 


1 
in the region 8 =0, and in the region 8< 0 it lies between the straight line — = C and the second of these hyper- 
bolas (and above the line A = 1 since U = © on the latter), = 


>. 


From Eq. (52) we conclude that as C decreases both hyperbolas (and together with them the entire level line) 
are shifted in toto to the region of arbitrarily large a and 6; this occurs in such a way that a/8- 0 (in view of the 
fact that the terms a/(28 +1) go to zero), 


We shall estimate the steady-state dispersion, From (34) we have 


i-L ’ (53) 
On the other hand D is nonnegative as the limit of the dispersions. If we move along the optimal curve to the 
origin, then L, M, N and Us go to zero as is evident from (47)-(51); therefore it is also true that D+ 0, If a +« 
along the optimal curve, thena/B - 0 (section 4), In that case we obtain 

























a ap? 4 a? 
De Halts 28+1 (1+ 6)* Ms TB = a ap? ‘ + a 
esd "R+T Wp + e+ ap Fi * ap + oR + AB a 
BL +B 


by intensifying inequality (53) by means of (47)-(51), 


The authors express their appreciation to A, A, Fel'dbaum for his statement of the problem and for his discus- 
sion of the results, 
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It is shown that the study of stability of periodic states in nonlinear pulse automatic systems can be 
reduced to the investigation of a linear pulse system with periodically varying gain or to an equivalent 
multiconnected linear pulse system with constant parameters, Thus the well-known stability criteria 
for pulse systems with constant parameters becomes applicable, The effect of random disturbances 
on periodic states is considered, 


The papers [1,2] dealt with the determination of periodic states, a period being a multiple of the sampling 
period in nonlinear pulse automatic systems (NPAS), The method of harmonic balance forms the basis of these 
papers, As shown in [2], by operating with discrete periodic functions instead of continuous [1], one is able to ob- 
tain the exact solution for periodic states with their half-period N = 1, 2, and also an improved approximate solu- 
tion for periodic states with half-period N= 3, However, the problem of stability of the obtained periodic states has 
not been tackled until now. In this paper, an exact solution to this problem is given, 


It has been shown that the study of stability of a periodic state can be reduced to the study of stability of the 
equilibrium position of a linear pulse system with periodically varying gain. The latter system can be made equiva- 
lent to a multiconnected linear pulse system with constant coefficients, The transfer functions of the multiconnected 
pulse system are obtained, to which the well-known stability criteria of ordinary pulse systems are applicable, The 
effect of random disturbances on periodic states is also considered, ; 


1, Equations of Nonlinear Pulse Automatic System 





We shall consider an NPAS consisting of a nonlinear pulse element (NPE) in series with a lumped-parameter 
continuous part (LCP); their block diagram is shown in Fig, 1, The equation of this system connecting the images 
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can be written as 


D {x (n}}=D {f (n}}—W* (9) D{® (z[n))}, (1) 


where D{x[n]}, D{/[{n]}, D{M(x[n})} denote two-sided discrete Laplace transforms of the error x{n] = 
x{n, 0], of the external force f{n) = f[n, 0] and of the output quantity from the nonlinear pulse element,respectively; 
x) is the nonlinear characteristic of the pulse element; W*(q) is the transfer function of the open-loop linear pulse 
system, that is of a system whose nonlinear pulse element has been replaced by a linear one, The expression 


D{x[n]}} = >) e—™z[n] (2) 


denotes the two-sided discrete Laplace transform? 


We assume that a symmetric periodic excitation ifn) of the half-period N is applied to the system, Such a 
discrete periodic excitation can be represented in the form of a trigonometric polynomial [2] 


aa N, . kr 
Finlk=> DP cer”. (3) 


k 
k=—N, 


We denote the periodic state response by x {n], where 


~ 1 Ns . gor 
a(n] => >; oe . (4) 
k=—N, 


We then obtain for an odd nonlinear characteristic to be used henceforth for the sake of simplicity, the formula 


~ 1 oy, jee 
O(zin)=> Sy be® . (5) 


k=—N, 


In the expressions (3)-(5),if Nisan odd number N, = N, but with N even we have N, = N-1, The dash in the sum- 
mation sign signified that the summation takes place only over the odd indices [2], The periodic state equation as- 
sumes the form 

D {& [n}} = D F{n}} —W* (q) DO @In))}. (6) 


The determination of symmetric periodic states was presented in [1,2], Below it shall be assumed that the 
period state is known, 


2. Equations in Terms of Deviations 











The ordinary approach is used in the study of the stability of periodic flv) 2h0 
states, Small deviations from the periodic state are considered, If, in due 
course, the deviation approaches zero, the periodic state is stable; if the devi- 
ation increases without bounds, the state is unstable, Equations in terms of 
deviations shall be formed, It is assumed that at a certain time moment a Fig. 1. 
small disturbance occurs in NPAS due to variation, say, in initial conditions, 
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* In the case of x[n] = 0 for n< 0, the two-sided discrete Laplace transform reduces to the one-sided one. 








This disturbance is fed to the input of a pulse element. Denoting the disturbance by fg{n], we shall write the 
external excitation at the input of the pulse element as 


f{n) =7{n) + fgln- (7) 





This additional disturbance f{n) causes a deviation £&[n] from the periodic state. Thus 


z(n] = z[n] + §[n}. (8) 


By substituting (7) and (8) into equation (1), and by subtracting the equation (6) of the periodic state, we obtain 
the equation for the deviation 


D{& (n}} = D {fg (n}} — W* (gq) PD (D(z [n) + E (n}) — D(x [n))). (9) 


Bu: for small values of the deviation £[n] we have 


@ (x [n} + &[n}) —D(e[n}) = O' (z [n}) E [n), 


PF an d® (x 
(x {n}) = 


x=sx [n] 


Consequently, the ‘deviation equation becomes a linear one with variable coefficients 


D{& {n}) = D (fy (n}} — W* (q) D(®" (x [n}) & [n}). (10) 


According to the Lyapunov method, the stability of a periodic state is determined by the stability of the system 
equilibrium position described by an equation in terms of deviation regardless of the neglected terms which corres- 
por . o a higher order terms of the expression on the right-hand side of (9)? 


But to the equation (10) corresponds a linear PAS with perio- 
dically varying gain k{nJ= (x(n), whose block diagram is 
f (n} Et PE PIz\n))-— LCP Cin) shown in Fig, 2, Thus the study of period state in NPAS has been 
reduced to the study of stability of a PAS with a periodically vary- 
ing gain (Fig, 2), 
































Pig, 2. 3. Equivalent Pulse Automatic Systems 


A direct investigation of a PAS with periodically varying gain 
is not possible because one is unable to find p{t [n}} from its equa- 
tion in an explicit form. However, the PAS with variable parameters can be reduced to one with constant parameters, 
This can be achieved in two ways. The first one is to replace a pulse element with its relative repetition period 
equal to unity and its variable gain, by a parallel combination of N pulse elements with relative repetition periods 
equal to N; these elements operate with a shift in relative time equal to unity and possess constant gains equal to 
©'[n), where n = a, 2,..., N-1 [3], 


A similar equivalent PAS was used by Pyshkin [6] in order to analyze the period state of a PAS with a width- 
pulse modulation. 


* Here the equations are understood to be difference equations; Lyapunov method was generalized to include the 
latter by E, Cotton [4], P. V. Bromberg [5] and others, 
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We shall avail ourselves of the second possibility. Consider x(n). Since oxjn) is an odd periodic discrete 
function, #(x{n)) must be an even periodic discrete function and it can be represented in the form of a trigonomet~- 
ric polynomial 







J ~ aft wa, : itn 4 a ing” 
O (z{nl) =O' (> SY ae "\=> DY due 
wile k=—M (11) 


k= 
where M = [N/2] is the integral part of N/2, and dey the complex amplitudes of the harmonics, are equal to 
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day = dyxe’ ™ re @’ (x (n}) e (Jk | <M) (12) 












. :9 N—1 eat 
diom = dsome’ + — 7 > QD’ (z[n})e"? (12") 


when N is odd, but sts 


















N—1 
diem = dom = or >) (In) eri (12") 


v=0 
when N is even, 





In the expressions (11), (12) and (12") the quantities ‘dy = dy and diy, = d, yy are real. 






By substituting (12) into the equation (10) and by inverting the order of the D-transformation and summation, 
we obtain 


jut a 
D{E[n}} = DZ (n}} — +W (9) 5 dy,D {e' * "§ {n)). (13) 
k=—M, 
By introducing the usual notation for the images 





B@=D€in)}), FH g)=Diyinp (4) 
and noting that, in view of the shift theorem [3] 
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Die * "gn == (g—j 7), (as) 









we can write the equation (13) as 


{1 


° . 4 . =- , = . Zhan 
={@M—-zW@) > d= G—ia)- (16) 


k=—M 


The equation contains the require image © (q) for different arguments (q-j) (2k/N). By putting in 
the equation (16) (q-j) (2km/N) where r = -M, -(M-1),..., °1, 0,1,..., M1, M for an odd N, butr= -M, 
-(M-1),..., 71, 0,1,..., M-1 for an even N, and in view of the periodicity 


& (9 + arn) = &*(q) 





having the relations 





e 2\k ° a. 
E*(q4/ a) = 2 (9+) Ital. *) when |ki+jri=N, (17) 
we obtain a system of equations from which the required images can be determined: 
For an odd N, that is N =2M+1 


© (0-17) + 2 (0-77) Dy dae QAM 
= F (q—j7) a ee Pn, Se fs cb, M1, Me 


and for an even N = 2M 


(0-1) +2 (C157) dae (oi 2GP*) + 


k=—(M—1) 
1 > ; ; . -2(M e° -2 
+3 W" (q) (dam-+dom) & (g — | -“AG*) = Fag — 137) 
(7 =—M, —(M—1),..., —1,0, 4..... M —1). 


Noting that in agreement with (12") 
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2 a 2 = dam 
and by denoting 
, d 
dom —_ aa ’ 
where 
, GPIB) ting 
dim = 5) O Gin) 
v=0 


is of the same form as the remaining complex amplitudes, ~ - write the equations in the case of even N as 


(0-9) +b) ST dawe(a—/2HM)+ 


=—(M—1) 
+7W ) dink ‘(q- j2uene) noi 2) (19) 
ene teh. cnc anh Bde co Mmm 4). 


The systems (18) and (19) of equations correspond to equivalent multiconnected linear pulse automatic sys- 
tems (EMPAS) with a certain number of pulse elements operating synchronously and cophasally. 


In view of the shift theorem for the independent variable q (see [3]), with the transfer functions w'(q — i=) 


are associated transfer functions W(q —j ae) » of continuous units; the block-diagrams of these EMPAS's can 
be represented in the form as in Fig, 3a, when N is odd, and when N is even as in Fig, 3b, 


Solving equations (18) and (19) for 8° (q), we obtain 


N—M—1i A* (q) " 
5° _ vy, et rR’ i al arn 
(9) 2 “ws ( W): (20) 
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where for an odd N, N = 2M+1 
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A4Y,M+1(q) are minors obtained by deleting the vth row and(M + 1)th column fromthe respective determinant, For 
the sake of brevity, one abbreviates W,,—W*(q +=). in the formulas (21) and (22), 


The concept of a transfer function of an EMPAS is introduced 


. _— A ai 
Ky, m+ (9) = ae (23) 
Then the equation (20) can be finally written as 
. + e . .2vn 
B@= > K.mulg§ (q—i-y . (24) 
v=s—M 


We shall now consider particular cases 


Where N = 1 we have M = [N /2] = 0, and consequently we obtain from (24) 





B (q) = Kw (9) % (@), (25) 
where, in agreement with (23) and (21) 
i 
Kvo(9) = 
o(9)= + 2 W* (q) (26) 
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Fig. 3a. 


Where N = 2, we have M = [N/2] = 1, and consequently we obtain from (24) 


8° (q) = Ki1.0(9) Ki (9 + ix) + Koala) % (@) (27) 


where in agreement with (24) and (22) we have 
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4. Stability Criteria of Periodic States 





For the stability of a periodic state it is necessary and sufficient that the deviation €(n] from the periodic 
state approaches zero in the course of time n, that is that the EMPAS be stable. 


The stability of an EMPAS, as it is known from the general theory of pulse systems [3], is determined by the 
position of the poles of its transfer functions K}, M+1(%- If the real parts of the poles are all negative, that is, if 
the poles are situated in the left half-plane of the complex plane of q, then EMPAS and thus the periodic state of 
NPAS will be stable. Ks: 


But the poles of the transfer functions Kis 160 coincide with the zeros of their denominators, that is with 
the roots of the characteristic equation 











An (q) = 0. (29) 


Thus the periodic state in an NPAS is stable, if the real parts of the roots of equation (29) are all negative, 


In order to investigate the distribution of the roots of the equation 4X; (q) = 0, the well-known stability criteria 
of linear pulse systems [3] may be applied. 


As seen from (25) and (26), it is necessary and sufficient for the stability of the highest frequency periodic state 
(N =1), that the equivalent linear system be stable, the latter obtained by replacing the nonlinear pulse element by 
the linear one with gain d, /2. 


When N = 2, the stability of an EMPAS is determined by using the usual stability criteria of linear PAS [3], 


5. Effect of Random Disturbance on the Periodic State 





We assume that a periodic state exists inan NPASand that a stationary random excitation f,[n) sufficiently 
smal} in its absolute value is applied to the input of the pulse element in addition to the periodic external excita- 
tion, Then the deviation from the periodic state €[n] due to the random excitation f4{n] will also be sufficiently 
small in its absolute value, The image [£(q) of this random excitation is determined from the previously derived 
equations (10) or (25), Hence, the conclusion is drawn that the study of the effect of random excitations in a perio- 
dic state can be reduced to the study of the behavior of a PAS with variable gain (Fig. 2) or of an EMPAS with con- 
stant coefficients (Fig. 3a, 3b). 


Consider EMPAS with constant coefficients, Turning from the images in (25) to originals, we obtain from the 
convolution theorem 


ln] = > >; wy, m—1(m) f.[n — m]}, 


v m=—oo 


(30) 
where v assumes integral values from -M to N-M-1, w,,, M-1[™] are pulse characteristics corresponding to the trans- 
fer functions K}, \447(q) and 


Qvx 


junj=e™ "4 inj =D{R (9 + (2*)} (31) 


4 


are random excitations applied to the inputs of the multiconnected pulse system. In formula (31) we understand 
Dp { } to be the inverse two-sided discrete Laplace transform [3]. 


As the system was assumed stable, the quantities Wy,M-1[ approach zero with the increase of n. We shall 
assume that fgn is a random interference such that its correlation function is Rg[n, m]. The cross-correlation func- 
tion of the external excitations is, by definition, equal to [7] 


co co 


R,,sin, m= \ | fen) fetm)p (frm faml) ded, = Ef, (nj, (mi), 
—2o —co (32) 


where p( ffs) is the two-dimensional distribution density, and f denotes the function conjugate to f, 
By substituting the valuesoff,, [n] from (31) into (32) we obtain 


j 
R,,,(n, m] =e = R,[n, m}. (33) 
We shall find the correlation function of the deviation € (n]) 


Ry (n, m) = E(& {n) E(m)). (34) 
By substituting (30) inio (34) and by inverting the order of averaging and summation, we obtain 
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co co 
R:[n, m| =>) >) Wy, M4114] >» Wy, M+11"q] E(f [n—m,], } {m—m,)| 


v.44 ™y=—0o M3 =—00 (35) 
or, in view of (32), we finally obtain 
na co 
R,{[n, m) =) s) Wy, M+1[™,] Zz Wy, M+1[M] Ry y(n — m, m— Ms}, (36) 


v,~2 m,=—oo m,=—0o 


where vand yp assume values from -M to N-M-1, 


The variance of the deviation which indicates the magnitude of the effect of random disturbances on the devi- 
ation is 


of = E(t" [n]) = Rz[n, m] (37) 
or,in view of (36) 
. ew e 
cG= > >> Wy, M+i1] > Wy, M+1 [M9] Ry. [n—m, n— Mg}. (38) 
v,u Mmy==—0o ™M,g=—0O 


When computing oF from the formula (38) as w, ,M+1[nJ approaches zero with increasing n, one only considers 
a finite number of terms in the internal sums. For the periodic state of N = 1, we obtain, by putting M = 0 in (38) 


SE= SY} Wola) SY) Woo (M2) Roo(n — my, n— my). (39) 


In this case the estimation of the effect of random disturbances on the periodic state is reduced to the study of 
the behavior of an ordinary linear pulse system randomly disturbed. 
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This article is devoted to the study of the amplitude quantization of a stationary random signal. It 

is shown that the quantization process can be considered as a process that adds stationary noise, which 
we will call quantization noise. We consider conditions under which the quantization noise can be 
assumed to be practically uncorrelated with the input signal and to have a spectral density that is con- 
stant over a very wide frequency range. Methods are developed for applying this theory to automatic 
control systems, 


Automatic control systems containing digital computation devices (DCD) are at the present being applied in a 
continually widening field. The basic advantages of such systems lie in the possibility of attaining practically unlimi- 
ted accuracy and interference rejection. The investigation of such systems, however, is very difficult, since in a con- 
trol system containing DCD the signal, in part of the channel, is in discrete digital form. The representation of acon- 
tinuous signal in a digital code is essentially nonlinear because of the quantization of the signal amplitude. 


There are several articles in the literature devoted to the study of amplitude quantization [1,2], The main re- 
sult of these articles is the treatment of the quantization process as one that adds quantization noise, where this noise, 
under certain conditions, has a practically constant spectral density, These articles, however, do not contain any ri- 
gorous theoretical basis for this result; in addition to this, they contain several erroneous conclusions. In particular, 
the quantization noise is assumed to be independent of the input signal. In the article by Widrow [2], conditions are 
derived, on the basis of a faulty premise, for which correlation in quantization noise is absent; the resulting condi- 
tions appear to be more restrictive than those obtained in the present article. Moreover, the articles referred to con- 
sider only the error analysis of quantization in closed coupling channels for pulse-code modulation, 


In the present article, we develop a theory for the amplitude quantization of a stationary random signal, and 
consider methods of applying this theory to the investigation of closed, digital, automatic control systems. 


The Problem 





The process of amplitude quantization is the nonlinear operation of representing a continuous physical quantity 
in digital form. The range of the input is divided into a series of adjoining intervals o;. All the values of the input 
that fall in a given interval o, are transformed into a certain integer that lies in the interval. A matheinatical model 
of a quantizer is constructed on the basis of the following assumptions. 


1, All of the intervals o; are of the same length (oj = o) and located so that the ends of the interval 0; coin- 
cide with the inputs (i-4)o and (i+4)o. 


2. If the input to the quantizer lies in the interval o;, then the output is io; thus the mean amplification of the 
quantizer is unity, 


3. There is no delay in the quantizer. 


Such a quantizer and its characteristics are shown in Fig, 1. Let the quantizer input be the known, two-dimen- 
sional probability distribution of a stationary, random process Wx(x;, Xz, 7) or the corresponding characteristic function 


Fox (Uj, Ug, t) = M (eT) — \\ efit Usx2) W 2 (Xi, Xq, T) dx, dx. 
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Then the two-dimensional distribution law for the output of the quantizer will be a distribution law for a two- 
dimensional random variable taking its values from the set of points with coordinates that are multiples of o. The 
probability density of this variable is the sum of terms of the form 


(i+) 0 (e+ )0 
W, (is, ko, t) = (x, — is, rz, — ks) ( ( W 2 (21, Xe, T) dx, de, (2) 
(i ~)eo (x- 1)s 


where i and k are integers, 


As can be seen from (2), the probability density W? is the 
product of periodically (with period o) repeating, two-dimensional 
x” 6 -functions with the function 




















X+ — Xy+ . 
2 2 
eal Q \ \ W 4 (2, Le, tT) dx, dx». 
x\— > Seo > (3) 
° -2o 
b When we consider the characteristic functions, it is not 
Fig. 1. difficult to obtain 


Foxe (Uy, Ug, T) = > > Fx (uy +i, ue + kg, t) 


i=—oo k=—oo 


sin > (us + ig) sin (ua + hp) (4) 





Flu +ig) — (us + kg) 


where ¢ = 2n/o,. 


Quantization Noise and Its Characteristics 





We now consider a quantizer as a combination in parallel of a linear amplifier with amplification coefficient 
unit and a nonlinear element with the sawtooth -shaped characteristic p(Fig. 2) 


We will give the name “quantization noise* to the random process defined as the difference between the input 
and output of the quantizer 


N= X*—X. (5) 


As we can see from Fig, 2, the quantization noise is formed at the output of the nonlinear element with charac- 
teristic p. 


























In the determination of the statistical characteristics of the quantization noise, we use the expansion of the 
function N = ~(X) in a Fourier series 





N=¥(X)= 2 y CN sin™ x, 
(6) 


Then for the autocorrelation functions of thequantization noise Rpyf T) and for the function R,,.( 7) defining the 
correlation between the quantization noise and the input signal, we have 























gyi 
Ran (t) = M(W (2) ¥ = > Yasin 2, sin»), we 
i=1k=1 
R 3 (—1)* 2nk 
xn (t) = M {x,¥ (2)] = = k M [2 sin <= 2,| e 
k=1 
(8) 
We now write the sines as differences between exponentials, and note that 
° . OF (uy, Ug, t) 
> 2x 
jM [zye?™*2] = Ou; adit (9) 
We finally obtain 
—_— tin oni k 
Rants) = $s 3S MMH eH, BH, 9), 
a (10) 
ko 
2nk 
Rea(t)= & 5 uh (uF: *) 
t\)= — 
an ( axel Ou, —_ (11) 
k#0 
We assume that the function F,,(u,, ug, 7) is finite, i.e, that it is identically zero everywhere outside the 
square |ur|< Mims | Ue |< em. Further, let the conditions* 
2 12 
Uim< =, tam < (12) 


be satisfied, 


Then both the expressions (10) and (11), for the quantization-noise autocorrelation and for the mutual corre- 
lation between the quantization noise and the input, become zero. Thus, when the conditions (12) are satisfied, the 
quantization noise is uncorrelated with the input signal, and is a stationary random process with uncorrelated values, 
In order to find the mean square quantization noise, we note that if the conditions (12) hold, then the moments of 
the quantizer output are equal to the sums of the moments of the output and the uncorrelated noise with character 
istic function (Fig, 3 a) 
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Fin (u) = ny (13) 
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To this characteristic function corresponds the probability density (Fig. 3 b) 


*In Widrow's paper [2], the more stringent conditions “im S = > Sen S = ‘ were obtained, 
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The mean square of the quantization noise is therefore 


co 
M {nj = | nW,(n)dn =<, 
12 
—co (15) 

It should be remarked that the conditions (12) do not hold in any physically realizable case, However, for the 
majority of probability densities occurring in practice (for example, for the normal distribution), the characteristic 
functions decrease rather rapidly, so that for sufficiently small steps, the quantization-noise autocorrelation function 
is very small for reasonably large values of 7, The spectral density of the quantization noise under these conditions 
is thus practically uniform in a very wide range of frequencies, its value being 


4n? ik — »—, t) dt 
k+0 


It becomes smaller when the characteristic function of the input signal decreases more rapidly, and in the limit tends 
to a function equal to zero for all frequencies, However, since the equality 


Sun (0) = Sm(0) = 2 SS HEH \ Fax (2% , 2 


1 f 2 
5s— \ Snn(@)do = Ran (0) = = , 
wind = 2 (17) 


always hold, the spectral density in this limiting case is the singular 4 -function on the frequency axis, 
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Fig. 3. 


Amplitude Quantization in Pulsed Systems 





Since the input data and the output results from DCD are not usually produced continuously, but only at defi- 
nite discrete times, it follows that in systems containing DCD, in addition to the quantization of the signal amplitude, 
there is also a quantization of the signal in time, The presence of time quantization makes such systems discontinu- 
Ous, pulse controls, 


In this connection, it is of interest to study the statistical characteristics of quantization noise subjected to a 
double time quantization, i.e. quantization noise that has passed through a pulse element [3], 


A convenient characteristic of a random process subjected to time quantization is the discrete spectral density, 
given by [3] 


S*(o) = >) R(m7’)e-irer, 


m 















where (m7) =R(t)/-—mr gives the discrete spectral density of the quantization noise the function 


. 3 
Snn (@) = Ran (0) = 'q3 ° (19) 


from which it follows that this discrete spectral density is practically constant over the whole frequency range, and 
that its value depends only on the size of the steps in the amplitude quantization. This makes it very simple to take 
amplitude quantization into account in pulse systems. 


The Case of a Gaussian Input Distribution 





To illustrate what we have said above, we consider the case when the input signal X to the quantizer is a sta- 
tionary random process with a normal distribution. It is known that the characteristic function of such a process 


F yx (uj, Us, t) = exp {- + (ui + u3) 8° + 2uyu,Rxx (x)1} (20) 


is not finite, but it decreases rapidly. If we substitute the characteristic function in the expression (10) for the auto- 
correlation function of the quantization noise, then after some transformations, we obtain 


wf et — is 
| (t) ars >) al ae 5 “k tiny ike (*) - 
i=1 k=1 
Bd? 4 dee, 4n®x? 
_ x: are 6 sinh — Rxx (t) + 
— 4yit ei is 
‘ aS XN 5 nal ah “(i Fin it kPa (t), (21) 
i=1 k=1 


where 6 =o" /6” is the depth of the quantization and 5” is the mean square of the input signal. 


If 8 « 1, we can neglect the double sum in (21), Further, if we replace the hyperbolic sine by its asymptotic 
expansion, we obtain* 


_ wrens xx(*)) 


Ran (t) = BB 
( Da Dae (22) 


In particular, for T = 0, we obtain the known expression for the mean square of the quantization noise 


BPS 1 Bex? ot 
Ran (0) = 553 2) a = ont © = TD: 


n=1 


(23) 


Let Rf 7) =( & sin 4)44, corresponding to a uniform distribution in the interval A of the spectral density of 
the input process power, Then the spectral density of the quantization noise will be 


4n*x? sin tA 


San (@) = 88S) : Oe a 74’ cos wrt dt. 


mo 


=. (24) 


When(sint4)/7™ is expanded in a series and only the first two terms are retained (which is permissible since 
the function(sinrA)/TA oscillates and the error is to a large extent cancelled out), we obtain 


*This result was first obtained by Bennet in {1}, The further calculations up to (26) are the same as those in the book 
by B, R, Levin [4], 
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cos wt dt 
n=1 0 (25) 


or, after integration 


_ Bo 


Ea E54 >t =. 
n=] (26) 


Snn (@) = 


It is clear from this formula that the spectral density of the quantization noise can be assumed to be constant 
up to the frequency 


A Vor | 
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a = 1,9 


(27) 


where y is the relative error in percent, The boundary frequency p», as can be seen from the formula (27), increases 
as the spectrum of the input signal widens and as the depth of quantization6 decreases, It should be noted that the 
condition that 6 be small is satisfied if the dynamic range of the signal exceeds by several times the size of the 
quantization step o, 


The spectral density of the quantization noise is given by the expression 


Snn (@) = San (0) = 8 By) 1 _o@ VBS 


2A 2n - sae 


n=l (28) 


The mutual correlation between the quantization noise and the input signal is obtained from (11) and (20) af- 
ter some transformation, and is 


co anh? 
Ren (t) = 2Rzx (T) >> (— 4° e 6 ; 
The (29) 


It therefore follows that it does not exceed in absolute value the number € for 
2x? 


B= oe 
§ (30) 





An Application to the Investigation of a Digital System 


We consider the digital-system scheme shown in Fig. 4, Here CD is the converter for transforming continuous 
quantities into digital form, the block diagram of which can be given in the form shown in Fig. 5, 


Following Ya, Z, Tsypkin in [5], we use the term “limiting 
pulse system" for a system differing from a digital system in the 
absence of amplitude quantization. Then if the quantization steps 
are small enough, the digital system in Fig, 4 becomes a limit 
pulse system, Acting on this system is not only the useful signal, 
——{ cD + lip) +-——y-~ but also two noise components, uncorrelated with the input signal, 

and due to the amplitude quantization in the input cirauit and in 
the feedback circuit, The discrete spectral density of each of 
these noise components is constant and equal to’/,,0*. In the in- 
CD vestigation of this system, all the methods for the analysis and 
synthesis of pulse systems under the action of stationary random 
perturbations were used, 
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Now, let some interference as well as the useful signal act at the input. It is not difficult to see that the addi- 
tion of this interference can only broaden the dynamic range of the input signal, This means that the spectral density 
of the quantization noise remains unchanged, and that the error caused by the action of the interference on the digital 
system is the same as the error that would be caused the the action of the same interference on a limit pulse system, 
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SUMMARY 


1, The amplitude-quantization process of a stationary, random, input signal can be considered as the process 
of adding to the input signal a certain equivalent noise, called the quantization noise [1,2]. 


2. The autocorrelation function of the quantization noise and the mutual correlation between the quantization 
noise and the input signal can be determined from the formulas (10) and (11) by the known methods of characteristic 
functions or from the probability density of the input signal. 


3, If the amplitude quantization steps are small enough in comparison with the mean square input signal, then 
the quantization noise can be assumed to be uncorrelated with the input signal. The exact conditions for the validity 
of this result are given by the inequalities (12), 


4, With the same condition that the quantization steps be small, the spectral density of the quantization noise 
can be assumed to be constant over a wide range of frequencies, and its magnitude can be determined from (16). The 
discrete spectral density of the quantization noise depends only on the size of the quantization steps, and is equal to 


5. For sufficiently small amplitude-quantization steps, a digital system reduces to a limit pulse system, at 
the input of which acts not only the useful signal, but also the quantization noise with spectral density */,,0°. 


6. If there is interference acting at the input of a digital system as well as the useful signal, then the error 
caused by this interference is equal to the error that would occur in a limit pulse system under the influence of the 
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EXTREMAL RECEPTION OF SIGNALS 
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Systems are considered in which the received signal is compared with the signal from an automatically 
adjustable internal generator. The adjustment of the generator is carried out by using an external con- 
trol system which operates on the minimizing principle of a criterion of differences between the genera- 
ted signal and the received one, The equations are obtained and the dynamic process is analyzed in a 
continuous system of this type. 


Receivers of various kinds of signals* usually belong to systems requiring some kind of adjusting operatic... 
Hence it seems natural to apply the principles of automatic self-adjustment to signal receivers, 


There are various possible ways of implementing a self-adjusting receiver. As in the case of ordinary hand 
adjustment, it is possible to have an automatic adjustment of passive networks of the receiver based on the principle 
of reaching the extremum of certain indicators of the quality of reception, 


In this paper an account is given of another principle of attaining extremal 
reception which is as follows, 


f.+ f The received signal contains a useful component f, and interference fp. In 

the general case, the noise of interference can be assumed to be an arbitrary random 

time function, Only at the synthesis stage of an optimum self-adjusting system is it 

3 necessary that the statistical characteristics of the interference fy be known, The 

Lay general form of the useful component f, is assumed to be known beforehand; it is 

Fig. 1. characterized by a number of paremeters x,; such as amplitude, carrier frequency, 
duration and phase of the pulse train, etc, Thus f, = f,{t, Xq, Xg...- Xgm)- 












f, 
(i r, 









2 





The initial information on the useful component of the received signal gives the composition of the parameters 
Xs only; the accurate values of these parameters at a given moment of time are not known,” The received signal 
is compared with the signal f,, from an automatically adjustable internal generator (Fig. 1). 


The construction of the internal generator is such that the form of the signal fc = fft, Xe Xge--- Xm) as well 
as the adjustment parameters x; of the generator correspond, at least approximately, to the form and the parameter 
constructicn of the signal useful component so, that the differences f,-fg can be made sufficiently small by select~- 
ing the values of the parameters x; within the admissible range of f,. 


The difference between the received and the generated signal proceeds to the deviation estimator, block 1, 
The block output is a criterion of the current deviation of the generated from the received signal, This quantity is 
fed to the optimizer 3, a system of extremal regulation which implements the search and leads to an adjustment of 
the generator 2 such that a deviations criterion is minimum, The output signal of the system of extremum reception 
is either the signal of the adjusted generator or only the adjustment parameters x;. 


A general outline of automatic approximation using a self-adjusting device was briefly developed in [4], The 
results of the experimental examination of one case of approximation with the aid of an automatic optimizer were 
given in [5]. 


In this work, continuous systems of extremal reception are investigated, the study being based on the above 
mentioned principle where the search proceeds by the gradient method, 


*The term “receiver” is used in the text in its wider meaning and indicates any receptor of information. 
**A signal with all the parameters fully given does not carry any information, 











Block Diagrams and Equations of the System 





The block diagram of a continuous system of extremal reception is shown in Fig. 2, Insofar as the search for 
the extremuin is concerned, the system is in every way analogous to the multidimensional continuous systems of 
extremal regulation with synchronous detection [1-3]. 


The search oscillations 6x; correspond to the adjustment parameters of the generator 2. If the extremum of 
the criterion is not attained, these oscillations produce synchronous variations in the magnitude of the criterion, The 
latter are extracted by means of synchronous detectors, thatis, amplifiers to which the search signals are sent as re- 
ference signals. These signals thus obta'ned pass through filters Wfj and integrators, and produce variations in the 
generator adjustments as long as the synchrom-u; components of oscillation appear of the deviations criterion, that is, 
until the extremum of this criterion is attained. 
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One can cite the scheme shown in Fig, 3 as an example 
LZ; Gs, ; of an adjusting generator. Here, there are generators with 
T constant adjustment producing m constant signals {), fgo..-- 
«+ «fg { “carrier” functions of the generator), Along the 





signal paths one finds the multipliers with gains as adjust- 

Fig, 2. ment parameters, The generator output signal is obtained by 
summing the signals after they have passed through the mul- 
tipliers, Thus the generator signal becomes in this case a 
linear functiou of the adjustment parameters 


m 


Ill, 2, %,..., 2m) = bs Xi fe; (¢) (1) 


i=1 


Before considering the equations corresponding to the system, one should note that different functions may be 
taken as the deviations criterion, The use of a squaring device as a deviation block driterion corresponds to a quad- 
ratic function 


U = (fy + fy — f?. m 
2 


The use of a full-wave rectified linear detector as a deviation block criterion corresponds to the absolute va- 
lue of deviations 


u= |f + fy —fel- (3) 


In all cases, the deviation criterion can be represented in the form 
u =Fif (t, Tsiy+++ » Ism) +n (t)—Jg(t, Tyee »Xm)] = 
=F ify Ol + 2fn () HAW Ax + >) fy (t) AnAz; 


i=1 i, jm (4) 
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Here Axj = x,; ~x; are the deviations of the signal parameters of the generators from those of the useful compo- 
nent of the received signal; the functions f(t), ftv) = fit) depending on time are close to the functions f,, fc: 


For many deviation criteria (that is, for many functions F) and for a number of functions f, and fg. the term 
®t, Xo,» +++» Xem) in the expression (4) is of a higher order of smallness in respect of x; than the preceding terms, 
Thus, for example, if the adjustment (modulation) of the generator is used, expressed by the formula (1), with an 
analogous modulation of the useful signal 


) (t, Tyin-++s gn) ves > tala (t) 


i=] 


and if the quadratic criterion (2) is selected, then the term © vanishes altogether as 


u= fi, (t) + 2f, (t) df(t) Ax +3 hij (t) Ax Ax;, 


(5) 
where in this case f(t) = f¢Av), fi ft) = fgj(t) gir). 
The deviations Ax; can be represented as consisting of two parts: the search oscillations 6x; and the deviations 


* 
Axi of the outputs from the integrations 


Az; = Ax; + 52. 


In accordance with the latter, the expression (4) assumes the form 


UF (fy) + 2fy (t) 3} f(t) Aai + 2f, () Salt)Oa+ 5 fas(t) AaiAzs + 


i=1 i=1 i, j=1 


+2 5) fg (t) Aeidas + SD) fj (t) dxdx; + ©. 
net hss (6) 


The output from the ith amplifier (Fig, 2) is equal to uéxj, Having passed through the filter whose transfer 


function is W,,(D) and through an integrator, the signal udx; _ represents the output from an integrator which, in ac- 


cordance with out notation, is X,j-AX}. Thus 


Px) We ,(D) 


or %i— Ax; = D 





(udz;) 


DAx; = — W, « (D) (ubx;) + Dai. 


Using the expression (6), we obtain a complete system of equations for extremal reception 


DAx; = —W 5 (D) 25) Ay (t) Axb2; 82; + fn (t) Se (t) da, 82, + 


v, j=l ver] 
+ 2fy (t) > fo(t) Aarda;, + S) fyj(t) AxpAzjdx, + 
v=] ¥, j=1 
+ >) M(t) da bajbx; + F [f, (t)] 6x; + O82; + Dry (i = 4, 2,...,m). 
v, j=1 (7) 


These equations are nonlinear in the general case, the coefficients being ordinary (fy(t), f yjit)) as well as ran- 
dom (fp(t)) functions of time, In the general case, it is not possible to integrate these equations as was not possible 
with other extremal regulations [2,3], Nonetheless, under some conditions, justifiable from the practical point of 
view, it is possible to obtain an approximate solution in a comparatively simple way. 








Quasistationary Extremal Signal Reception 





A behavior of considered systems, ruc: that their parameters xj, Axj vary slowly compared with ft), f(t), 
5xj;, shall be called quasistationary or, in other words, we shall say that tie spectrum of the functions Ax; (t) are 
sufficiently narrow. 


The lower bound of the admissible spectrum width of the functions Ax; is, of course, determined by the spec- 
trum width of the signal parameter xsi, i.e., by the spectrum width of modulation, The latter, in turn, depends on 
the quantity of information transmitted in a unit of time. In a quasistationary state, the spectrum width or the spec- 
trum “lines” of the signal useful component are known to be less than that of “carrier” frequencies or search frequen- 
cies, 


In a quasistationary state, the passband of the filters Wr; is made to be narrow compared with the signal or 
search frequencies. 


Quasistationarity of the behavior is the first assumption that we have made in order to be able to analyze the 
process, 


The second assumption, closely related to the first one, is as follows. We represent the time functions entering 
the equations (7) in the form 
fog (t) BarjSxxq = fj (t) bxjdx; + Ej, 
fy (t) fu(t) Bax; = J, (t) J. (t) bx,62; + Not, 
fn (t) fu(t) 8x1 = fy (t) fo (t) bx + Sa, 
fv (t) B74 = Ju; (t) O45 + 805, 


fj (t) Sxvdx ja, = fy; (t) d2,b2j521 + Ty 50. 7” 


Here the bar denotes the time average (constant components) of the respective quantities. 


This condition implies that the spectra of the “rapidly varying” quantities £, ,¢, €, 7 are outside the trans- 
mission bands of the filters Wfj, and the components are completely suppressed by those filters. 


As far as the functions €, «, and T only are concerned, this condition can always ben ensured for a quasista- 
tionary state if an appropriate choice of the frequencies of search vibrations is made. 


But with regard to the functions n and ¢ in which the interference ft) is present, the above condition can be 
ensured when the interference spectrum is found within the spectrum domain of the carriers of the useful signal fy(t), 
this being, however, the most important case, In other instances, one can obviously reduce it to the latter by intro- 
ducing a sufficiently narrow-band filter at the system input, 


We shall also assume that the search oscillations take place about zero and that they are orthogonal 
[623 for v=i, 
\ 0 for V=i. 


It is beyond the scope of this work to analyze fully the compatibility of these assumptions. For a number of 
specific instances, these assumptions have proved compatible (see the example cited below). 


bz; = 0, 62,6264, =0, b2,d2; 


f, (t) f(t) 
I sj = 2 3 
—&—-—_— 264; W; ,(D) D 
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If our assumptions are valid, the equations (7) are transformed into 


DAx; + 2d23W44(D) S' fa (t) Ax, = 


v=] 


= Dis —W 4i(D) [2h (t) fi (t) dx? + F Uf, (0) Ox; + OBz)). 


(9) 


The spectrum of the quantity F[fp(t)}6 x; lies, in view of our assumptions with regard to fp(t), outside the trans- 
mission band of the filter Nfj. 


The quantity 6x; can certainly be neglected for small deviations Ax;, in view of its higher order of small- 
ness, 


Therefore, in accordance with our assumptions, the equations for a quasistationary extremal reception degen- 
erate into the following system of linear equations with constant coefficients: 


DAx; + 2823W 1 (D) 3 Fa) Ax, = Dag, — 2824W 44 (D) (fy fi) 
v=] 


(é == 1, 2,..., m). (70) 


The block-diagram corresponding to these equations is shown in Fig, 4, where hfs | is the matrix of coeffi- 
cients, 


It is apparent from the equations (10) and from the bloclediagram that the processes of extremal adjusting of 
the generator are, under our assumptions, analogous to the processes in a linear multidimensional follow-up system 
whose input quantities are the parameters of the useful component of the input signal x,j, and whose outputs are the 
parameters of the generator signal x;. 

The role played by interference is taken here by the functions fp(t) f(t), This enables one to perform the ana- 
lysis and the synthesis of extremal reception in accordance with out assumptions on the same lines as the well-known 
methods for stationary linear systems, 

The stability of the extremal reception is, in view of (10), determined from the characteristic equation 

r - - 
=m I + os Sima 
2dx°W 5 (2) 
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We shall say that the “carrier functions” f(t) are orthogonal if the conditions are satisfied 


W for i= v, 
0 for i=v. (12) 





fu(t) = 
When the carrier functions are orthogonal, the system (10) of equations breaks up into m independent equations 


DAz; + 2822 fiW (D) Axi = Day — 2622W f(D) {/,, (t) fa(t)) (13) 
(én 4,2,...,@) 


and the equivalent multidimensional follow-up system (Fig. 4) also breaks up into m independent channels. 


The synthesis of the investigated extremal reception system consists in finding transfer functions of filters such 
that they ensure the lowest level of errors Ax; for the known characteristics of the random function f,(t) and of ran- 








nr — 


dom or ordinary functions x,(t), It is customary to formulate the synthesis in this way and it is not necessary to dwell 
on it here, 


As in other extremal regulations problems, these can be investigated in a sequence of approximations (of in- 
creasing complexity), starting from the simplest case of a quasistationary state [equations (10)] to the general case of 
forced self -adjustment (equations (7)], Such an attempt was made in the paper [6] for one specific case of extremal 
regulation, The development, however, of some more powerful methods of analytic investigation of extremal regu- 
lation processes is a different problem. 


It would be useful to clarify the important question of noise stability of quasistationary extremal reception com- 
pared with other well known types of reception. 


Let the useful component of the received signal be of the form 
m 
hy = >} Xvi v(t), 
v=] 
where the carrier functions f,At) are known without error. 


One is able to apply successfully to such signals a correlation mode of reception or a synchronous detection 
method known as a “noisestable” method, In order to have such a reception, the signal f,+fp is multiplied by f(t), 
and afterwards a filter is applied, Thus, in case of correlation mode of reception, one has 


r= Wi (D) Shh +h Oh 9]. 


v=] 


As previously, transmission band of the low frequency filters W, can be adopted such that they will suppress 
completely the variable components mj; of the functions 


ROK) =AOKO + ru. 


2m Wo(D)] Sint + by (OKC) 


v=] 


In=hOKO® ( =1,2,...,m). 
(14) 


When comparing these equation with those of quasistationary extremal reception in (10), it is easy to establish 
that with a specified relation between the transfer functions Wfi(D) and WD), the two systems of equations are essen- 
tially the same. 


Thus, if the carrier functions are orthogonal [as in (12)], in order that the equations (10) and (14) be identical 
it is sufficient that the condition 
28a? W ¢, (D) 


W;(D) = —= ; 
1 +282? f,, We, (D) 





(15) 


be satisfied, that is, that it is enough that the transfer functions of the filters correspond to the operators of the closed 
and open loop single-circuit systems, 


Hence, the conclusion that potential noise-stability of an extremal reception is equivalent to that of a correla- 
tion reception, The correlation reception can, however, take place only if the carrier functions f(t) are known accu- 
rately, or more precisely, when these functions are implemented in the receiver. 
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An extremal receiver operating as a self-adjusting system does not require a preliminary reproduction of the 
carrier functions. This is a special feature of extremal reception and constitutes an advantage over correlation recep- 
tion. An example will be given in order to illustrate this special feature. 


Example, Let the useful component of the received signal take the form 


f, = Asin (wt + q). 
We discuss two cases, 


1, The carrier frequency w is fixed and known exactly, The information is transmitted in the form of ampli- 
tude A and phase ¢ subjected to small variations not known beforehand, 





We put the useful component of the received signal in the form f, = Asin (ut+¢) = xg sin ut+xXg cosut, with 
Xy = Acos¢, X = Asing being the transformed parameters of the received signal, and f,(t) = sinut and ft = cosut 
the carrier functions, 


The formula representing the self-adjusting generator of the extremal receiver shall be taken as 
ig = 7 sin wt +4- x2 cos wt, 
where xX, and X, are adjustment parameters, 


With the quadratic deviations criterion (the block 1 in Fig, 2 now is a squaring device) the function u becomes 


u=(/, +f, — fg)? = (Az: sin ot + Az, cos wt + 7, 2, 


The search oscillations are considered harmonic with two different frequencies «4, and Us: 
Az, = Az: +6,,8in@t, Az= Az, + 5, Sin Get. 
Then the equations (7) of extremal adjustment assume the form 


DAz; = — W «,(D) (32, (2Az; sin* wf sin @¢ sin @,t +- Az Sin 2@¢ sin @ei sin w,t + 
+ Az, sin 2ot sin t sin w;t + 2Az, cos* wt sin @st sin @;t) + 
+ 2/,,5,, (sin wt sin @¢ sin @,t + cos wt sin wet sin @;t) + 2f,,9m (Az; sin wt sin @,t + 
te Az, cos wf sin @;t) +-6,, (Az;? sin*w?é sin @,t +- Az: Az, sin 2@¢ sin w,t + 


+ Az,? cos? wt sin w;t) + 83, (sin® wt sin*® @¢ sin wt + sin 2ot sin @t cos @,¢ sin @;t) + 
+f 5, 8in @t] + Dz, (i= 4, 2). 
If the frequencies 44, u, and w are selected such that the combination of frequencies 


%+0,@+01, Of, 204+0;,, 20+, @40:+@:,, ©420;, ©4202, 20 + O + O,2@ + 20, 20 + 2s, 
as well as uy, and Us are all outside the passbandof the low-frequency filters Wa and Wg, then in a quasistationary 
state, all the components in the square bracket are very nearly suppressed except 


i . ; : ‘ 2 
and the equations become 2 Az; Om» fn %mSinwt (for i=1), f, of, coswt (for i= 2), fy Sm Sin @, t, 
« pont . 
Dhz, + > 85, W¢, (D) Az, = Dr, — Wg, (62, f, sin ot +3, f£ sin w,t), 
"as ue A 
DAz, + > 5;,W 5, (D) Ax, = Dag, — W 5, (82, f, cos ot +-6,, f2 sin st), (16) 


which is in full agreement with the equations (9) (the carrier functions satisfy here the orthogonality conditions (12), 
and @ = 0), 
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If the spectrum of the interference is close to that of the useful signal, then the last components in the right- 
hand sides of (16) simply vanish and the equations of the quasistationary extremal reception assume the form 


‘4 
2, = W,(D) (4 Ts) —I, sin wt), y= Ww (D)(- Tso — Sn cos ot) é 


Exactly the same results can be obtained for correlation reception, if the received signal xysinut+xgcosut+f, 
is multiplied by the carrier functions sin«t and cosut, and then sent through the narrow-band filters W, and Wg. 


Because the correlation reception is more straightforward, it would be futile to use extremal reception, The 
extremum reception in a slightly different form should, however, be used when the carrier frequency is variable or 
unknown, and the correlation reception would not now be applicable. 


2. Carrier frequency w variable and not known precisely, Let the useful component of the received signal be 
as before, that is f, = Asim ut+e) = Asing,, , =wt+g with the frequency unknown and in addition subjected to 
some random variations. 





We take the frequency “G, or more precisely, the current phase ¢ = “t+, and also the amplitude Ag: fg 
= Agsingg as adjustment parameters of the internal generator; this is ensured by an appropriate construction of the 
generator different from the one in the previous case, We assume that during the automatic adjustment the ampli- 
tude difference AA = A - AG, and also the difference of the current phases Ay = 9, - % are small. Then 


—{ =Asing, — A-sing,= AAsin 9, + Acos@, Ag. 
Ys 3 4,5i9& 5 $ 


By denoting Ax, = AA and Ax, = AAg and with other conditions being the same (harmonic search oscillations 
of AA and of AAg, quadratic criterion), it becomes obvious that the dynamics of this process is identical with the 
one we have just considered, 


In particular, one can write down the equations analogous to the ones in (16): 
4 
[p 4. > 8. Wey (D)| AA=DA— &. Ww fi (D) (fg sin @t), 


4 
[p +> on W fe ‘D)| AAg = D (Ags) — 82, W 5, (D) (fn 008 @t), 
where ¢, = ut, 


It can be seen from these equations that when the transfer function Wy(D) has no poles at zero, that is when 
the corresponding filter contains no integrators in its single-circuit scheme, the system of extremal reception is static 
with regard to the change of frequency w w Dg, of the received signal. In order that it be astatic in frequency, it is 
necessary that the filters contain an integrator, As regards the change in amplitude A, the system is astatic when even 
the transfer function Wa(D) has no pole at zero, 


In conclusion, it should be said that the signal reception is also solvable in the case of no a priori complete 
information on carrier frequencies, if appropriate follow-up systems are used, 


Compared with the latter, the investigated method is certainly more flexible and versatile which, generally 
speaking, is inherent in self-adjusting systems. 
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CHARACTERISTICS WITH THE HELP OF A GENERALIZED FAMILY 
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pp. 739-747, June, 1961 

Original article submitted December 22, 1960. 


A method is presented for plotting amplitude and phase frequency characteristics of electric networks 
for which the complex transfer coefficients can be represented in a canonical form, Tables for 28 
circuits and the formulas which connect the circuit parameters with the parameters of a generalized 
family of logarithmic characteristics are given. 


Amplitude and phase frequency characteristics* are the most important characteristics establishing the basic 
properties of electronic amplifiers. Knowledge of AFC and PFC are especially important for measuring and evaluat- 
ing (operating amplifiers) in which a large amount of negative feedback is used so that the problem of ensuring sta- 
bility and a broad passband becomes complicated, Determining AFC and PFC takes much time during the design of 
such amplifiers because their transfer functions are generally of a high order, and a calculation must be repeated 
several times before a desired shape of the characteristics is successfully attained. 


It is much easier to use the method of logarithmic frequency characteristics [1 to 4] which permits the AFC 
and PFC of elementary sections to be constructed with the help of a set of prepared model curves, from which, in 
tum, the characteristics of an entire network are obtained by a simple summation. However, the method of loga- 
rithmic frequency characteristics as developed, was designed chiefly for the analysis of automatic regulating systems, 
and for the calculation of the elementary sections which are most prevalent in such systems, In electronic amplifiers 
a diversity of types of elementary sections are encountered for which families of characteristics or model curves are 
required, In addition, the association between circuit parameters and section parameters often proves to be so com- 
plicated that it is difficult to trace the direct dependence of the shapes of the AFC and PFC on the circuit parameters, 


Presented below is a method of plotting AFC and PFC which permits a substantial simplification of the comput- 
ing process so that the AFC and PFC of almost all the simplest networks encountered in amplifiers is reduced to only 
two generalized families of logarithmic characteristics. Moreover, each characteristic is determined by one parame- 
ter « which characterizes the attenuation of the transfer coefficient of the section for a frequency variation from 0 
to*. The parameter « is physically descriptive, and is especially convenient for the synthesis of correcting networks. 


1, Generalized Families of Logarithmic AFC and PFC 

In Tables 1 to 4 there are given the basic circuits of the elementary RC and RL networks most prevalent in 
electronic amplifiers. All the networks given are related to a class of networks whose transfer coefficient modulus 
changes monotonically as a function of frequency, but not faster than 20 db/ decade. The circuits of Tables 1 and 
2 have a negative attenuation (the transfer coefficient increases with frequency) and have a leading phase angle. 
The circuits of Tables 3 and 4 have a positive attenuation and a lagging phase angle. 





It can be shown that the complex transfer coefficient of any circuit of the first group (a> 1, Tables 1 and 2) 
can be reduced to the following canonical form** 


*Hereafter for brevity we will designate the characteristics referred to by AFC and PFC respectively. 
**It can be easily shown for each circuit that by substituting in (1) the values of X,a, and T from the tables given, 
we will get an expression for K(jw) in the general form as written. 












1 
+ jX 
K (jX) = K (00) S—*__., (1) 


+iX 








a—i 





where X = wT is a dimensionless variable; w is the circular frequency; T is the equivalent time constant which is 
given in terms of circuit parameters in Tables 1 to 4; a is a parameter characterizing the change of | K(jX)| for a 
change of w from 0 tos, 
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TABLE 1 
Circuits K (0) E (co) « r 
RC, om 
C1 Ci (Ri + Rs) a—4 
Cy+Cs| Re(Ci+ Cs) — RC R, (Ci + Cs) 
— "™! RC, — RC 
1 1+ Re CaRs 
C, 
Cy 
0 fore) R(C, + C 
R; C, Ci + Ce ewe 2) 
C; 
“3 
R 0 1 0° CiR, 
2 
C; R; CR, ua 
Ry C; R, Ra (Cy + Co) a—4 
Ci+Ca| Ait Re} Ci (Ai + Aa) = GR C; (Ri + Ra) 
T? ™~ "9 CR. — CR 
C 
o— 
if Ci Cy 
2 6+ Gy + C: 1 i+ Ci CR; 
os Lemma, 
Cy R, 
R: 
% . Ri +R ~ CiRs 








(2) 








From (1) we obtain 

















j 2 
| K (7X) | = K (00) (saa) + 
(sz) +> 
tng => a -§ ’ 
G@—ipt * 


where ¢ is the phase of K(jX). 
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Circuits K (0) K (co) a T 
R,C (+ 1) = 
Ry Cy C; (Ry + R,) mt a 
hy + Ry Cy, + Cs R; (Cy + C.) ait Ry (RC, —_ R,C,) 
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4 0 0 CR; 
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RC; (=—1)= 
sé Ci Ry Riz (Ci + Cs) ~ (a 
% | G4+Ca| Rit Ral Ci(Ai + Ms) _ Ca (RiCs — RiCs) 
T% Ci+ Cs 
Ry 
Rs R; 
f, 1 Rit} Rts ort 
T%, 
o—— 
C, R, 
| Ci Ci 
[? Ci + Ce 9 ° Cok Ci + Cs 








By changing the value of the parameter a, it is possible to construct generalized families of logarithmic AFC 
and PFC for the circuits of the first group (Tables 1 and 2) which are shown in Figs. 1 and 2 for K(e) =1. Inasmuch 
as the quantity lgK(©) is a component on a logarithmic scale, the initial reading along the axis of ordinates is shifted 
by a quantity lgK(@) which can be taken into account with the help of the K(e) scale shown in Fig. 1. 


As will be shown below, the generalized families of characteristics make it easy to plot AFC and PFC for any 
circuit of Tables 1 and 2, as well as for any complicated network that can be represented by a series connection of 
such circuits which are separated by detecting elements* (for example, electron tubes), 


Ww sare 


*In a number of cases it is permissible to break circuits down into sections when they are not separated by detecting 
elements provided that the output impedance of a preceding section is at all frequencies much smaller than the input 
impedance of the following section. 
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For the second group of circuits (Tables 3 and 4) having a positive attenuation (a<1) it is convenient to use 


TABLE 4 
























































Circuits K (0) K (00) e T 
L(t _,) 
Ry Ll, I4(Ri + Re) Ay \a *) 
A+ Aki ht+h| &A(h+) — (a/ Ri —Ty/ Ra) Ri 
Ri + Ry 
Ye ek ee 7 
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Rt R 0 0 R ar 
R, 1 2 it MM, 
“; 
1 0 0 3 
k, Ra 
a(t 1) = 
Ll, R,_ | Re(ly4+1,) Ai \a vd 
+1.) Rit Re) Lh(A1 + Re) = (3 -— +) I, 
Ry Ry/in+l, 
Fee oer Ra L 
Ri + Rf Ai + Ry Ry 
= 0 0 pie 
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the following expression for the transfer coefficient 


whence 


K (jX) = K (0) 


| K (jX)| = K (0) 


1+ jx 





i—a 





1+/x 


i 


i—a 











(5) 
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Fig. 1. Generalized family of logarithmic AFC 
of the first group of circuits (a> 1, Tables 1, 2). 
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Fig. 2. Generalized family of logarithmic PFC 
of the first group of circuits (a > 1, Tables 1, 2). 


Tables 1 to 4, 
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_ a ys 
“+ (1— a)? (7) 


Families of generalized logarithmic characteristics for the second group of circuits (Tables 3 and 4), plotted 
according to formulas (6) and (7), are shown in Figs. 3 and 4. 


It is easily seen that the families of characteristics 
shown in Figs. 3 and 4 differ from the families of Figs. 1 and 
2 only in the direction of reading along the coordinate axes, 


It can be shown that the expression for the complex 
transfer coefficient (5) transforms into expression (1) upon 
replacing the variable X by -/X"* [in this case & = K()/K(0) 
becomes equal to K(X'=0)/K(X'=e) i.e. Ya'}. This 
implies that for the family | K(jX)| only the direction of 
reading along the X axis is changed, while for the family 
y(X), in addition, the sign of g is changed because ¢(X) is 
an odd function. 


Thus, for a useful plot of the logarithmic AFC and 
PFC of any of the 28 circuits given, it is sufficient to have 
either the families of Figs. 1, 2, 3 or the families of Figs, 1, 
3, 4, i.e. to unite Figs, 2 and 4 by drawing additional scales 
in the opposite direction. 


We note that networks containing both L and C cannot 
be calculated from the families derived because the maxi- 
mum rate of attenuation of | K(jX)| with a change in fre- 
quency for such circuits exceeds 20 db/ decade, and their 
complex transfer coefficients do not reduce to the forms of 
Eqs. (1) or (5). 


2. Method of Plotting AFC and PFC with the 
Help of a Generalized Family of Logarithmic 








Characteristics 

The logarithmic characteristics shown in Figs. 1 to 4 
are plotted as functions of the dimensionless variable X =wt 
= 2n fT. For practical plotting of Afc and PFC it is desirable 
to have the characteristics as direct functions of the circular 
frequency w or frequency f expressed in cps. Inasmuch as 





lg X = lg 2nT +tan/, (8) 
the logarithmic characteristics are shifted along the axis of 
abscissas by the quantity lg 2"T, which is equivalent to a 
shift of the origin of coordinates by—lg21T. Allowing for 
this, it is possible to plot beforehand a scale of time con- 


stants that show where the origin of coordinates should be located along the axis of abscissas for a given value of T. 
Such scales are laid off along the axes of T in Figs. 1 to 4. 


There is presented below the sequence for constructing logarithmic . AFC and PFC for any circuit given in 


First, from the formulas given in atable, T, a, K(0) and K(«) are determined, and two coordinate grids are 
prepared as shown in Fig. 5° (separately for AFC and PFC). Then the coordinate grid of AFC is superposed on the 


*The scale of the coordinate grid must conform closely with that shown in Fig. 5, seeing that it agrees with the scale 
of the curves of Figs, 1 to 4, The length of the axes depends on the particular case. It is expedient to make them 
cover a considerably larger range of variation of f, K and ¢ than is shown in Fig. 5. 
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corresponding generalized family of AFC (Figs. 1 or 3) so that the axis of ordinates will stand on the mark of scale T 
corresponding to the time constant T® of the circuit, and the axis of abscissas is fixed at the mark corresponding to 





the given value of K(e) or K(0) along the vertical scale. After this, the curve for the required value of a is plotted 
on tracing paper. In analogous fashion, the PFC is plotted, but for this case the axis of abscissas on the tracing paper 
must always coincide with the axis of abscissas on Figs. 2 or 4. 
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Fig. 3. Generalized family of logarithmic AFC 


of the second group of circuits (* < 1, Tables 3, 4). 
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Fig. 4. Generalized family of logarithmic PFC 
of the second group of circuits (a< 1, T ables 3, 4). 





Fig. 5, Example of a coordinate grid for plotting 
AFC and PFC 


After plotting the characteristics of all the series-con- 
nected sections, their simple sum is taken, and the AFC and 
PFC of the comlete network being considered is obtained. 


The process of plotting AFC and PFC can be made 
even easier by preparing two sets of model curves, 


The accuracy of the AFC and PFC plots is determined 
chiefly by the accuracy of plotting the characteristics of the 
generalized family (or model curves), and also by the errors 
which sometimes occur during the analysis of a circuit into 
elementary sections. As experience has shown, for a plot of 
logarithmic characteristics by this method, the total error 
can be held to a value not exceeding 0.5 db for the AFC and 
2to % for the PFC. 


To the advantages of the present method besides high 
accuracy must be added its simplicity. For a comparatively 
large class of networks there are used only two families of 
characteristics, which significantly speeds up the process of 
plotting AFC and PFC, In addition, the use of the descrip- 
tive parameter @ = K(#)/K(0) significantly helps in the se- 
lection of parameters for the various sections during the syn- 
thesis of a complicated circuit for a given AFC and PFC. 


The method presented was developed chiefly for the analysis and synthesis of AFC and PFC of electronic amp- 
lifiers for which the equivalent circuits can be represented in the form of a series connection of elementary sections 
as given in Tables 1 to 4. However, the advantages of the method enumerated above apply to the calculation of any 
other electrical networks which can be represented as a series connection of circuits of the class under consideration, 
especially to the calculation of correcting networks for systems of automatic regulation and servo systems. 


"in 


case the required value of T lies outside the scale of Figs, 1 to 4, the T scale should be extended in the 


necessary direction with the adopted logarithmic scale. 





SUMMARY 


1, The method presented permits the reduction of calculations for amplitude and phase frequency characteris- 
tics of a whole class of elementary sections (28 circuits) to the selection and copying of the required characteristics 
from a family of generalized logarithmic characteristics which are calculated beforehand. Because of this, great 
tapidity of synthesis, while retaining accuracy, is ensured. 


2. The method set forth is most convenient for the analysis and synthesis of the AFC and PFC of electronic 
amplifiers, However, it can be successfully used for the synthesis of electric correcting networks of automatic regu- 
lation systems and servo systems. 
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The concept of recursive event, formalizing the concept of "an event, effectively described by a 
verbal description ,” is introduced, The problem is posed of recognizing those recursive events which 
may be represented on a finite automaton. A theorem on the algorithmic insolubility of this problem 
is proved, At the same time, the concept of representation of a recursive function in an automaton is 
introduced and a criterion for recognizing representable functions is given. 

As this group of problems is being considered in “Avtomatika i Telemekhanika” for the first 
time, certain well-known concepts and facts are also explained. 


1, Recursive Events 

In the notation used in [1, 2], the symbols pj from the alphabet {p}={po, p,,..., Pr-s denote the input states 
of a finite automaton. By definition, the number of states is finite and equal to r, In the present article we consider 
only nonautonomous finite automata such that r = 2. 





An arbitrary sequence composed of a finite number of symbols pj, for example 


Ps ’ 
PePiPs» 


PsP1PsP1PsPa 
etc., 


will be termed the input sequence, The number of symbols constituting the sequence will be termed its length. 
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Let in any way, in an infinite ensemble { E} of all conceivable input sequences (each of which has arbitrary, 
but finite, length), a finite or infinite subset { S} of input sequences be distinguished, Let now there appear in suc- 
cession certain input sequences from {E}. We say that an event S has occurred if the sequence appearing belongs to 
{ S} , and that the event has not occurred if the sequence appearing does not belong to {S}. If the set {S} is finite, 
the event S can be given simply by counting all sequences which enter into { S}. If the set {S} is infinite, this 


method of prescribing the event is excluded, and the set must be defined in some other way. Thus, for example, the 
events: 


1) there has appeared the sequence containing in any arbitrary position the symbol ps; 
2) a sequence longer than three has appeared; in which the fourth symbol is p3; 


3) there has appeared the sequence terminating with the symbol pg, and containing nowhere the symbol p, — 
although they occur during the appearance of sequences from infinite sets, are fully characterized by the above 


verbal description, and it is always possible to recognize whether the event has occurred in the appearance of any 
sequence. 


In general, to prescribe an infinite set {S} of sequences signifies giving such a “verbal description.” as will 
permit for arbitrary preassigned input sequence, to answer the question, does it belong to the set { S} or not ? The 
existence of such a verbal description testifies to the existence of an algorithm for recognizing the membership of an 
arbitrarily appearing sequence in a given set {S}. In this sense the giving of an algorithm for such recognition is 
the natural and most general “language” for the effective prescription of infinite sets {S} , i.e., for prescribing events. 


The purpose of the present work is to clarify to what degree a finite automaton can realize such recognition. 
To obtain an answer to this question it is necessary to formalize the intuitive concepts introduced above of “verbal 
description” and “recognition algorithm." A formalization of the concept “algorithm,” as is well known, is given 
by the theory of recursive functions and Church's thesis. 


Let us assign uniquely to each sequence a certain arbitrary number, such that to different sequences there 
correspond different numbers. It is possible, for example, to use the following coding of sequences by numbers, 
written in the r-nary system. We write the number 1 and immediately after it from left to right we write in sequence 
the indices of the symbols p; in the prescribed sequence, reading from right to left. Thus, for example, with r= 5, 
the sequence P2PoP;Psz is put into correspondence with the number 13102, while the sequence p2P3P9P4P;, with the 
number 114032, in the quinary number system, In this coding system, to each sequence there corresponds a single 
number, which in the r-nary system has the first digit 1, and vice-versa, to each number, having in the r-nary system 
the first digit unity, there corresponds a single sequence. Numbers whose first digit is not equal to unity do not cor- 
respond to any sequence. Therefore, to the universal set {E} of all conceivable sequences there corresponds the set 
A(E) of all r-nary numbers beginning with 1, while to the set {S} , a certain subset A(S) of these numbers. 


This problem of recognizing the membership of a given sequence in the set of sequences {S} is now substituted 
by the problem of recognizing the membership of the corresponding r-nary number (beginning with 1) in the set of 
numbers A(S). 


Let us consider the integer function $x) defined on the set of r-nary numbers beginning with 1, and taking 
the values 1 for all numbers in A(S), and 0 for all remaining numbers. If there exist an algorithm for recognizing 
membership of the sequence in the set {S$}, the values of the function s(x) can be calculated for all x in A(E). 
In this sense the function #s5(x) is calculable. 


In agreement with Church's thesis, asserting the partial recursiveness of each integer function, calculated for 
values of an argument from a certain set, the function $5(x) is partially recursive,* defined in A(E). 


*In the theory of algorithms an integer function is termed generally recursive (or frequently, simply “recursive”), 

if it can be obtained from certain more simply given functions (identity functions, constant functions, etc.) only by 
means of certain defined operations (superposition, primitive recursion, etc.), The function ~(x), defined in a cer- 
tain subset of integers, is termed partially recursive if there exists a recursive function ¢{x), taking on the samie 
values as ~(x) for values of the argument taken from this subset. Church's thesis declares the equivalence of the con- 
cepts: recursive function, calculable function, algorithm (cf. (3, 4). 
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Let us construct the integer function ¥g(x), already defined among all integers, the values of which coincide 
with ¥g(x) for x€A(E), and is equal to 0 for all other x. 


In regard to the existence of an elementary algorithm for testing the membership of the number x in the set 
A(E) (writing x in the r-nary number system and determining the presence of the 1 in front), and the partial recur- 
siveness of s(x), the function ¥¢(x) is recursive. From this it also follows that the function ¥,(x), taking the value 
1 for x€A(E) and 0 for the remaining x, is recursive. 


A function given for all values of an integer argument, and taking the value 1 in a certain set {M} of ar- 
guments and 0 outside it is termed the characteristic function for {M} sets of numbers whose characteristic functions 
are recursive are termed recursive sets. 


We term the event S recursive if the set A(S) is recursive. The universal event E, consisting of the appearance 
of any arbitrary sequence of symbols from {E}, is recursive, since the function ¥p(x) is recursive. 


Above it was shown that in the case where the event S is such that there exists an algorithm for recognizing 
the fact that event has occurred (i.e., the membership of any given sequence in the set { S}), the characteristic 
function ¥g(x) of the set A(S) is recursive, i.e., A(S) is a recursive set. 


Thus (this is a direct consequence of Church's thesis:), each event for which a “recognition rule” can be stated, 
is recursive, and vice-versa. 


The concept of recursive event is the broadest and most natural formalization of the intuitive concept "event 
which can be effectively described," “event for which a recognition rule can be given,” etc.* 


2. Representation of Events in an Automaton and the Problem of Recursive Event 
Representability 








Les us consider a finite autoton A (cf. [1, 2]), described by the recursive relation 


“Pti = F (xP, p?), p=—0,1,2,..., (1) 


where F is a single-valued function, put into correspondence with the pair of symbols p; (from the input alphabet 
{p}) and x; (from the alphabet of states { x}={x9,. . ., %q—1}), a symbol from {x}. In the set {x} we 
define a certain subset M. We shall assume that automaton (1) begins operating at the time p = 0 in the initial state 
»(0) = »°€{x}, under the action of the sequence p°p'p”..., developing from cycle to cycle. In each given cycle 
p the sequence of input symbols appearing from the start of operation of the automation to the present cycle, has a 


finite length p + 1. 


Let us assume now that an event S is given, in the sense in which this term was used in Sect. 1. It can happen 
that up to the time p the input symbol sequence belongs to the set {S}, i.e., that at this instant "the event has 
occurred." The event S is termed representable in a finite automaton if there exists such an automaton A, initial 
state %° and set of states Mc {x}, such that the sequence of input symbols up to the time p belongs to {S} if and 
only if the automaton is in one of the states belonging to M at the time p + 1. 


The automaton A can have an output decoder 
a” = @ (x), 


which, for all symbols yj EM, puts out »? = 1, while for the remaining symbols, -,P = 0. Then the appearance of 
the symbol 1 at the output of this “automaton with decoder” indicates the occurrence S. 


The problem naturally arises of determining the class of events representable in a finite automaton, and their 
discrimination from events which cannot be represented in any finite automaton. 


* If the requirement of effectiveness of describing the event is not put, it is possible to prescribe the events in the 
still broader language of recursive-denumerable sets. In such a description, however, the algorithmic unsolvability 
already extends to the problem of recognizing whether an event has occurred with the appearance of some given 
sequence, and, a fortiori, and problem formulated below of recognizing the representability of events is algorithmic- 
ly unsolvable. 
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Kleene [5] has introduced the concept of regular event, and has proved that for an event to be representable 
in a finite automaton, it is necessary and sufficient that it be regular. 


In the terminology of the present work,” the following definition can be given for regular event. 
Three operations on sets of sequences are considered. Let there be given the sets of sequences {S } and {Q}. 
A. The union | = S|jQ__ is the set containing all sequences occurring in at least one of the sets S or Q. 


B. The product W = S x Q is the set containing all sequences obtained by the adjoining to each sequence from 
S at the right one each of the sequences from Q. 


C. The interation W = S* is the set containing all sequences obtained by adjoining at the right to each sequence 
from S any number of arbitrary sequences from S itself. 


The set consisting of an arbitrary sequence of length 1 (a single symbol p;), is termed an elementary set. By 
definition, elementary sets, and all sets which can be obtained from elementary sets by any finite number of the 
above three operations, are regular. 


If the set {S } is regular, then the event S corresponding to it is also termed regular. 


It can be considered that an event is given directly in the language of regular sets. In this case, by virtue of 
Kleene's theorem, no question arises as to whether it car be represented by an automaton. However, an event at 
the input can be given in any other manner, for example, by a verbal description, giving the instruction for recog - 
nizing the occurrence of the event. To answer the question whether this event is representable in an automaton, it 
is necessary to clarify if it is regular. Thus, the problem arises of defining the regularity of recursive event. It is 
easily proven that any regular event is recursive. But the inverse statement is not true, since recursive irregular 
(i.e., nonrepresentable!) events are known. 


The problem of present interest may be formulated as follows: does there exist an algorithm permitting, for, 
an arbitrary recursive event, an answer to the question of its regularity (i.e., its representability)? We term this 
problem the problem of recognizing the representability of recursive events. 


3. The Algorithmic Unsolvability of the Problem of Recognizing the Representability 





of Recursive Events 








In this section we shall be concerned with the phenomenon of algorithmic unsolvability. A problem is termed 
algorithmically unsolvable if the nonexistence of an algorithm solving it can be proven. By virtue of Church's thesis, 
the concept “algorithm” corresponds to the concept “recursive function.” This signifies that in each concrete case 
all possible initial data and results of solving the problem are coded by integers, and the existence of an algorithm 
solving the problem corresponds to the existence of a recursive function, which takes on values coding the results of 
solving the problem when the argument takes on values coding the “initial data." Therefore, the problem is al- 
gorithmically unsolvable if such a recursive function does not exist. 


Theorem 1. The problem of recognizing the representability of recursive events is algorithmically unsolvable. 


* Kleene has operated with sets of input tables, but to each table there may correspond a sequence of input symbols. 
For Kleene these sequences have developed, as it were, "from right to left," beginning from the instant p 
+ P(p—2), p(p—1), Pp (P), 


in contradistinction to the present work, where the sequences have developed “from left to right" 


p (0), ep (1), p(2)... 


In Kleene's terminology this signifies, in particular, that we are considering primitive events, and the defini - 
tion of regular operation changes somewhat (see below). A similar presentation is given in [6]. 
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After the present authors gave a proof Theorem 1, it was found that the statement contained in this theorem 
was given earlier without proof in 5° of B, A. Trachtenbrot's work [7]; we shall, therefore, in what follows term this 


theorem Trachtenbrot's theorem. 


To prove Trachtenbrot's theorem, we shall first formulate a more restricted problem than the problem of recog- 
nizing the representability of recursive events, and we shall prove that even this narrower problem is algorithmically 
unsolvable. This in itself will prove the theorem. Let there be given a recursive function p(t), defined on a set of 
integers, and taking on values from the finite set of numbers {0, 1, ...,.r—1}. We consider an automaton A with 
input alphabet {po,.... pr—;}- Of all possible input sequences we shall note especially the sequence 

Poy? 
Peco) Poxay? 


Poo) Proc) Poe) 

etc. 
where py,;y is a symbol from {p} whose index coincides with the value of the recursive function ¢(t) att =i. We 
consider the event S%, consisting in that up to the given instant there has appeared at the input to the automaton any 
one of the above-defined sequences. We say that the automaton represents the recursive function g(t) if it represents 
the event s?. In other words, the automaton represents the recursive function y(t) if at each instant of time p the 
state of the automaton %(p) belongs to a given set M if and only if the indices of all input symbols appearing be- 
tween t = 0 and t = p are sequences of values of the function ¢(t). 


Obviously there exist recursive functions certainly representable (for example, any periodic function, since for 
if the event S? is regular) and certainly nonrepresentable (for example, the function a(t), taking the value 1 for 
t= 1 and 0 in the remaining cases). The problem therefore arises of recognizing the representability of recursive 
functions. It is easily seen that this problem is a particular case of our general problem of recognizing the represen- 
tability of recursive functions. 


The function ¢(t) will be termed “finally periodic” if there exists a pair of numbers r and T, such that for 
arbitrary t= 1, the condition 


e(t+T)=9(t), 
is satisfied, i.e., the value of the function for t = r repeats periodically with period T. 
Theorem 2, The function ¢(t) is representable in a finite automaton if and only if it is finally periodic. 


We note that Theorem 2 permits constructing an unlimited number of examples of nenrepresentable events. 
What is essential for us here is that this theorem reduces the problem of recognizing the representability of recur- 
sive functions to the problem of recognizing whether a given recursive function is finally periodic. But since even 
this more restricted problem is algorithmically unsolvable,” the proof of Theorem 2 completes the proof of Trach- 
tenbrot's theorem as well. 


We now proceed to prove Theorem 2. 


Proof of Theorem 2. We shall first prove that if the function ¢(t) is finally periodic, it is representable in a 
finite automaton. For this purpose we prove that the set {S?} of sequences for whose appearances the event occurs, 


is regular. 
Let to the period of the function ¢(t) for t = r these correspond the sequence of input symbols Phy? Py 2"? 
Prep" Then the set {S”} can be represented in the form 


S°=(p. 4) X Poy % ae X Pores) x (p,,* roe X Pap) * [p,,U (P, * P,,) U ts U (Py, * roe P,,)! , 
i.e., it can be expressed in terms of elementary sets of sequences, using the regular operations, Q.E.D. 


We shall now prove the converse statement, i.e., that if the function ¢(t) is representable in a finite automaton, 
it is finally periodic. Let A be a finite automaton having q different states, representing the function ¢(t), with initial 
state %» a property of the state M. 


* The proof of this fact is contained, for example, in Chapter 10 of [8], where it is carried out for the case where the 
function g(t) takes on two values; the extension of this proof to the general case is obvious. 
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Let us consider the input sequence of arbitary length p > q: 
Pato) Prony * * * Poca * * * Prey (2) 
and the sequence of automaton states corresponding to it 
x% (1)%(2)...%(¢+ 1)... (p+ 4), (3) 


satisfying the recurrent relationship (1). 


Since automaton A has only q different states, in the sequence 


% (1) %(2)...%(¢ +14) (4) 


of length q + 1, at least two symbols will coincide, i.e., there exist numbers a, B(a< B, 1< a<qti,i<c 
<B<q+1) , such that ; 


% (a) = x (8). (5) 
Let us consider the function 
-_ t _ ” 
om ={e0 for O<t<a, ) 
@(¢{+B—a) for a<t< oo, 


whose values for t= 0, 1,...,a—1 are defined by the values (0), @(1),...,@(a—1), , while for t= a, 


a+1, a +2,..., by the values (8), p(B + 1), p(B+ 2) +--+» ie., the sequences (a)m(a+ 1) .. .p(a +s) 
and @ (8) p(B + 1)...@(8+ s) coincide for arbitrary s. 


We shall now consider the sequence of input symbols 


P Foy Pay * * * PSa—1) PHay* * * Popp * > + PSipy (7) 
and the sequence of states 
xx (1)...%(a—1)x(a)...*(p + 1) (8) 


corresponding to it, by virtue of (1). 


Since the first a symbols in sequences (2) and (7) coincide (since by definition $(t)= g(t) for t < a) and, 
in addition, the first symbols in (3) and (8) coincide, then for t = a %(t)= »(t) and, in particular, taking into 
account % (at) = % (a) = x (8). 


It is now easy to see that for arbitrary s the sequence 
% (a) x (a +41)...%(a+s) 
coincides with the sequence 
%(B)% (8+ 1)...%(B + 8), 
since both of these sequences are the result of the operation of automaton A for identical input sequences 


Poa) Pati)’ * * PSats)’ 
Pap) Po(p+1) 7 Po(p+8) 


and identical initial states (a) = %(B ). 
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Since automaton A represents the function ¢(t), then for arbitrary p in the sequence (3) %(p)€M, while by 
virtue of the above proof, only those symbols are contained in (8) which are contained in (3), it follows that X(p) €M. 


This signifies that the event S? occurs both in the appearance of the sequence (2) and in the appearance of 
the sequence (7). But since the set {S%} contains only one sequence of each given length p, it follows from this 
that (2) and (7) fully coincide for arbitrary p, i.e., 


 (p) =9(p). 


For p <a, this is automatically satisfied by virtue of (6), while for p= a this is necessary, in order that 


(Pp) = 9(p + B—2). 


This signifies that the function ¢(t) is periodic in the last analysis with period T = B- a, Q.E.D. This com- 
pletes proof of the theorem. 
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CANONICAL METHOD OF SWITCHING CIRCUIT SYNTHESIS* 


A. Sh. Blokh 


(Minsk) 

Translated from Avtomatika i Telemekhanika, Vol. 22, No. 6, 
pp. 756-764, June, 1961 

Original article submitted November 17, 1960 


A new method for the synthesis of contact networks is described. Upper limits are given for the 
total and mean numbers of contacts. 


The automation of relay control-circuit synthesis requires the development of simple and improved synthesis 
algorithms. 


In this article we present a method for the synthesis of switching circuits which we term canonical. The 
canonical method for electronic circuits has been described in [1]. The basic concept of the canonical method is 
the canonical table of Boolean functions. 


* The basic results of the article were presented 19 November 1959 at the Section of Cybernetics and Mathematical 
Logic of the All-Union Conference on Computation and Computers. 
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We remark that the canonical table permits the unified treatment of the methods of universal multipoles, the 
method of stages [2], the graphical method of V. N. Roginskii [3], and the construction of twopoles by the method of 
O. B. Lupanov [4]. 


Assignment of the conductance function of a two-terminal contact network. The structural conductance of a 
two-terminal contact network is defined by a Boolean function of binary arguments y = f(xy, %, ..., Xp) 
We shall adopt the following convention for describing Boolean function. We put the set x, = G4, % = Qy,...,X_ = @, 
into correspondence with a number N-= @, + 2@)_4+... + gate, It is obvious that 0 = N =2"-1, Let us define 
a row of 2" points. To each point we assign an ordinal number between 0 and 2"~*. Over the point with index N = 
=, + 2ap-14+... + qn-tq, we write the value of the function y = f(x, %,....%p,) forx,=,i., f( a, 
%, ...» &,). For example, for the function y = XyX%_ + Xy% we obtain the row 0110, We omit the indices of the 
points, as they are defined by the positions of the points in the row. 





Synthesis by the canonical method is carried out in two steps. We first construct the canonical table from the 
Boolean function, and then realize from it the circuit. 


Canonical table. The canonical table for a Boolean function y = f(x,,X2,...,.X,) has a triangular shape and 
contains (n + 1) rows. 





The top nth line contains 2" indices, the bottom, a single index. The top row constitutes the row of the func- 
tion to be realized. We now describe the method for constructing the canonical table. Assume that the kth row read 
from below, has already been constructed. It contains 2k indices. We divide the kth row from left to right into pairs. 
Below each pair (a, b) in the kth row we put a number c in the (k—1) th row, where to different pairs (a, b) we assign 
different indices c, to identical pairs, identical indices c. 


To the pair (a, a) we assign the index a. If in the (k—1)th row it is found that o = c, and c; is to the left of 
©, then c, is marked in some manner to indicate that in further synthesis o, will be ignored. 


A table constructed in this manner will be termed canonical. For the function y = xy% + %% we have the 
following table: 


10 O01 
2 3 


The two-terminal contact network, corresponding to a canonical table of (n + 1) rows, is realized in the form 
of a series connection of n blocks rp, where each block r, contains the contacts of only a single relay X;. The in- 
dices of the kth row, as the outputs of r,. To obtain the circuit of ry, it is necessary to connect the index c, defining 
the input and standing under the pair of indices (a, b), by the contact x; to the index (node) b and by the contact 
xX}, to the index (node) a. If a = b (consequently, c = a), the index c is short-circuited to the index a. If a = 0 (b=0), 
index c is not connected with a (b). To obtain the outputs of r,,, it is now necessary to connect identical indices 
(nodes) of the kth row to a single output. The basis for the method is given in the Appendix. 


Example 1. Realize the Boolean function 


Y = LyL_Xy + LyQFq + TyTy%y + TH TeTs- 
For this function we obtain the following canonical table: 


01101001 


The top row gives the function itself. The circuit is shown in Fig. 1. 
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Undefined states. In the canonical method of synthesis it is easy to calculate the number of contacts neces- 
sary, from the canonical table. For this it is necessary to determine all distinct indices contained in the table, begin- 
ning with the (n—1)throw. Two contacts correspond to the table entries a, b, c, for which a #b, ab # 0; a single 
contact, if a # b, ab = 0; in the remaining cases there are no contacts. 





Therefore, if in the table of the function there are undefined states, it is necessary to complete the function 
f(x;, Xg, ..+» Xp) in such a way as to obtain the smallest possible number of different indices in the canonical table. 
Naturally, we try to have identical indices in the lower rows of the canonical table, since in this case the upper rows 
of the table will have the smallest number of distinct indices. 


Let us give one possible method for completing the function. We divide the top row of the table into two parts, 
and attempt to complete the function to make these parts identical. We then divide the upper row into four parts 
and attempt to complete the remaining undefined state to have the smallest number of distinct quarters, then into 


eight parts, etc. 
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Example 2. A Boolean function in the arguments x,, %, Xs, X4 is given by the row-110- 1-111 - 0° 011. 
Dots denote undefined states (constituents). Using the above completion method, we construct the canonical table: 


-110-1-111-0-01 1 
wee 6 Ba 4 
ny XS 4 
5 6 
7 


The corresponding circuit is shown in Fig. 2. (In this case index 5 could be short-circuited to index 3 due to the 
presence of 1 in the pair). 


Permutation of arguments. The realization of the twopoles depends on the order selected for xX, X, ..-» Xp» 
i.e., on the choice of auxiliary function N = x, + 2xX,_, +... + 227 4x. 


In this function the order of arguments defines the coefficients with which the x; enter into the function N. 





Let us consider the changes in the top row with transposition of the arguments x), , 1, Xk +2- 


Consider one of the nodes in the kth row and the nodes above it in the (k + 1)th and (k + 2)th rows. In the 
(k + 2)th row we have four nodes (Fig. 3). With transposition of x;+;, Xk+2 the nodes interchange position, ag, as, 
since these nodes correspond to the values 1,0 and 0,1 in the arguments of x, , y, X;42- From this it follows that if 
the indices of the (k + 2)th row are divided into quarters, corresponding to the indices of the kth row, with trans- 
position of x; + 4, Xk +2 the inner pairs of the quarters exchange positions. Consequently, to obtain the ap a 
tions of the upper, initial row, with transposition of x, +4, X; + 2, it is necessary to divide the top row into ok large 
parts (number of indices in the kth row), then divide each part into four smaller parts, and interchange the second 
and third small parts in each of the larger parts. 
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In the synthesis method described here, the order of arguments has been strictly observed, which would corre - 
spond to circuits consisting of the series connection of individual blocks of contacts from single relays in a defined 
order. In a certain sense, this could be considered a defect of the method, since it would 
limit the class of circuits by which a Boolean function is realized. It is true that well- 


Nn known fairly effective methods for the synthesis of bridge networks [2, 3] have this defect. 
q; "ey, In this canonical method, it is possible to remove this restriction. For this, it is sufficient 
; A, 1; to transpose in the top row not all parts, but only individual parts corresponding to indi- 
‘ey vidual nodes of the kth row. The validity of this statement follows from the above con- 
A siderations, illustrated in Fig. 3. 
Nf! This generalization exceeds the limits of the previously considered circuits and can 


lead to circuits with a sma!ler number of contacts. In this sense the method is better. But, 
Fig. 5. like each improvement, it complicates the method of partition, which, however, can also 
be programmed, 





Synthesis of (p, 1)-poles. The contacts of a (p, 1)-poles is given byp Boolean functions yy = y4(Xq, %,..-» Xp)-++s 
Yp= f,(Xy, Xe, . + + » Xp). It is necessary to construct such a multipole, that the conductance between the ith input 
and output be f;(x;, Xg, . . . » Xp). The method of constructing the (p, 1)-pole is analogous to the method of construc- 
ting a twopole. 


The canonical table is constructed analogously to the above-described case for p= 1. In the top row we write 
the p2” values of the p Boolean functions in n arguments. In the bottom, zero row, we obtain p indices, corresponding 
to the outputs of the circuit. Realization is carried out in the usual manner. 


In a particular case, where the fj are pairwise disjoint, i.e., fjfj = 0 for i # j, it is possible to employ a single 
table. In this case, we put in the top row not 0,1, but the indices of the functions 1, 2,3, ..., p and 0. For exam- 
ple, for the functions y, = XyX_Xg + XyXpXg, Yo = XyXeXe + XpXpXe, Vg = XpRXy + XyXpXy + XyXoXy, we Obtain the following 
table: 


The table is synthesized in the usual manner. 


Introduction of auxiliary arguments. Introduction of an auxiliary argument reduces the volume of the canonical 
table. This reduces the required memory capacity, which in turn leads to reduction of machine time, and permits the 
synthesis of circuits with a larger number of arguments. The auxiliary argument can sometimes be found directly from 
the verbal description of the Boolean function. This, in a certain sense, characterizes the invariant properties of the 
given Boolean function. For example, for a symmetrical Boolean function it is possible to take as the auxiliary ar- 
gument the ordinary sum of arguments. 





In the general case the finding of an auxiliary argument for a Boolean function is a problem of linear program - 
ming. 


Let there be given a Boolean function of the auxiliary integral argument y = f(s), s = 9(X,, %, . - « » Xp)- 


The general principles of synthesizing the circuit do not differ from the above synthesis method: we construct 
the canonical table and then synthesize it. 


Let us give the method of finding the canonical table. 
Since s is a function of binary variables, it is representable in the form of a sum of monomials with integral 


coefficients. Let 5 — Y} 446;, , where the aj; are integral coefficients while the o; are monomials in x;, X;. 
f=1 


We put - 
Pe se = Di 046, 5, = 0 (A <k<m). 


i=1 


. We write in the top row of the table, in increasing order, all values 








Sm = 8. Below this row we write the values of sy —, also in increasing order. In the same way we construct all re- 
maining rows. In the last, lowest row there will be a single value sp = 0. In the table thus constructed we join 
each index y of the kth row by two line segments to the indices y and y + a, 4 ; in the (k + 1th row. Above the 
indices %, of the top row we write the corresponding values of the functions. The previous indices sj can now be 
struck out; they were necessary only for plotting the line segments. At the nodes of the table obtained we put, ac- 
cording to the above method for constructing the canonical table in the general case, new indices. Thus, we ob- 
tain the canonical table. The technique for constructing it is very simple. In fact, it is not necessary to write the 
indices s,. It is only necessary to write the values of the coefficients a. 


The synthesis of the canonical table is basically carried out as before, except that if previously, on the right- 
hand segments, we marked x; and on the left X,, we now, on the segments to the right with a; > 0, and the left, 
with a, < 0, form a series network corresponding to the monomial 0}, while on segments at the left with a, > 0, 
and the right with a; < 0, a parallel network, corresponding to &},. This method is easily extended to circuits of 
other logical elements. 


Example 3. For the function 


-| 1 for s=0,1,3, 
0 for s= —1,2, 


where s = -xX, + 2X»_ + Xs, we have the table: 


0 1 4 0 1 
er a ZO NZ 
2 \ 3 2 
ee 4 

4 5 


\/ 
6 


The corresponding circuit is shown in Fig. 4. 


The canonical synthesis method, inclusive of partitioning, has been programmed for the "URAL." We have 
obtained canonical tables for Boolean functions of seven and eight arguments. 


On the syntheses of twopoles. We indicate a slight extension of the canonical method for the synthesis of 
twopoles. 


If the canonical table of the Boolean function is divided into two parts along the (n—m)th row, in which there 
ate p different indices, we then obtain in the lower part the canonical table for a (1, p)-pole in (n—m) arguments 
while in the upper, P “small” tables for the (p, 1)-poles in m arguments. If the “small” tables are pairwise disjoint, 
i.e., the functions d by them are pairwise disjoint, then these p tables can be reduced to a single table [see the 
special case of synthesis of the (p, 1)-pole]. 


Example 4, Let there be given the canonical Table 1. 





10010100 
2 330 4504 
4 5 
6 


The first row defines two disjoint canonical functions with apices 4,5. We reduce them to the single canonical 
Table 2. For this we transpose the index of the apex of a given table to the locations of the unities in the top row 
of the given table. After this the table is constructed in the usual manner. 
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For the initial function 10010100, we obtain as a result Table 3. The corresponding circuit is represented in 
Fig. 5. 
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If only a certain set of “small” tables is pairwise disjoint they can be reduced to a single one synthesized separately 
from the remaining tables. This extension of the method is efficiently combined with expansion of the initial table. 


Expansion of the table is its representation by r such tables, such that the Boolean sum of their top rows be 
equal to the top row of the given table. This corresponds to representation of a Boolean function by the sum of r 
Boolean functions. 


We now describe a special form of expansion which leads to the construction of twopoles proposed by O. B. 
Lupanov [4]. 


The top rows of 2™ indices 0.1 in the "small" canonical tables above the (n—m)th row of the given table are 
separated intor parts with t indices each, except a possible last part with a smaller number of indices. It is obvious 


that 
r= [7 |44e 


We expand a given table in r tables My, Mg, . . . , Mr. To the top row of Mj we transpose the jth parts of the 
"small" tables, the remaining indices of the row Mj will be 0. It is obvious that the “small” tables with differing 
Mj are pairwise disjoint. Therefore “small” tables, each taken from a different Mj, may be reduced to a single 
table. There will be not more than 2° such reduced tables, i.e., not more than the maximum number of distinct 
"small" tables in the individual M;. Synthesizing the reduced tables of m arguments and r tables in n—m arguments, 
we obtain a realization of the prescribed Boolean function in the form of the construction proposed by O. B. Lupanov. 


APPENDIX 
Basis of the canonical method. For the conductance function y = f(x, %, . . . , X,) we introduce n + 1 inter- 
mediate integral arguments, satisfying the following relations: 





ro = const, 

ry = Gi (Fo, 2), 

re = Po (ri, 72), 

rn = Pn (rr %n)s 


4=Tn, 


(1) 


where the y; are single-valued functions. 


The two-terminal contact network corresponding to the function y = f(x,, %, . . - » Xp) is realized in the form 
of the series connection of n blocks, where each block r;, contains the contacts of only a single relay X;. Let us 
consider the synthesis of a single block. 


Here and below the square brackets denote the integral part. 
































The function r, = 9, &—1» Xk) defines the kth block. Let the function Ty, take on the values aj, ag, ..., am 
while the function rj, takes on the values by, bz,...,by . The values aj, a, ..., a, are interpreted as the input 
nodes of the kth block, while the values b;, bz, . . . , by , as the output nodes of the kth block. 


The conductance between the input a; and the output b,, is taken equal to unity with x, = $4, if the equality 
by = (54 44) is satisfied; in the contrary case the conductance between the input a; and the output b,, is taken 
equal to zero. The input node to the first block will be ro. 


Since there are two indices in the top row: 0, 1, in the following block there will be two outputs. In fact, the 
zero Output is not realized. 


To each node aj we assign a binary variable a; according to the following rule: If r,_, = a;, then a = 1; 
if rt, , # aj, then aj = 0. Analogously, for the output nodes b,, we introduce the binary numbers 6,,. Let us find 
the relations between a,6, x,. We put ry, =b,, and list all values of a; for which the equalities 


RP; (@;, 1) ss... Px (4; 1) = b,, 
VP, (a; O)=...= Pr (4; 0) oF b,. 
are valid. 


From this it follows that 


By = ty (4, +... + 04) + xy (aj, +... +5.) @) 


Iq 
To construct the connection circuit in the kth block, it is necessary to connect the node b,, to the nodes ai 


eres i by the contact x; and to connect the nodes a -»@; to the node b,, through the contact X,. 


py ig 
In this realization 6,, is the conductance between the input to the circuit and the node By. 


This is obvious for the first block. Let the statement be valid for the a;. In this case, there follows from 
relations (2) the validity of the statement for the 6, as well. 


Since the ¢), are single-valued functions, then, taking r,— , we find only two indices b,, y Py,» such that 
by = (aj 1), by, = %_ (aj, 0). Consequently, 


ey «8, + 2,8, (3) 


From this we find that each node a; is connected by the contacts x, and X, to not more than two nodes by, 
and b,,, which is the guarantee that sneak paths not be formed. 


It is easily seen that the canonical table defines a system of ¢), functions. In fact, the indices of the (n—k)th 
row are the values of the argument r;, and the correspondence between the pair (a, b) in the kth row and the indices 
c in the (k—1)th row signify that ¢}(c, 1) =b, 9,(c, 0)= a. Obviously, there are g, single-valued functions. We 
note that r,, takes on only two values: 0, 1, while rp takes on the single value, given in the lowest row of the table. 
By the construction of the table, y = rp. 


To synthesize the circuits we do not require passage to the functions gy ,,. Synthesis is carried out directly from 
the canonical table. 


Since the indices of the (k—1)throw denote the inputs to the kth block, while the indices of the kth row signify 
the outputs of the kth block, then, employing formula (3), the following rule can be formulated: The node of the 
(k—1) th row standing below the pair (a, b) in the kth row is connected by the contact x, to the node b, and by the 
contact x,—c to the node a. If a=b, then c = a and the node c is short-circuited to the node a. Identical indices 
of the kth row are joined in a single node. If a = 0, then c is not connected to a; the same is valid for b = 0. 





Generalization 1. Let there be in the top row of a canonical table in n arguments not one nonzero index, 
but q distinct indices. This signifies the case of the (1, q)-pole, when its Boolean functions fj(xy, %, . . . » X,) 40 
not have common constituents, i.c., fj(xj, %, . . « » Xp)fj(Xy, Xe, - - - » Xp) = 0 for i #j. In this case the construction 
of the canonical table and its synthesis do not in any way differ from the case of a twopole. This follows immediate 
ly from the fact that, beginning with any arbitrary intermediate row of the twopole table, there exists a canonical 
table and the corresponding circuit of this type. 
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Generalization 2. For the p Boolean functions y, = f;(xy, X2, ...»Xp)--- Yp= fy (X, %, - + + » Xp) Of the 
(p, 1)-terminal contact network the system of intermediate arguments (1) is constructed analogously to the system 
(1) for a single Boolean function. In this case, however, rp takes on p distinct values and r, = z, where the function 
z= F(T, Xy, Xe, . . - » Xp) Satisfies the relationship yj = F(i, xy, %, .... Xp). Graphically this system of functions is 
constructed in the form of p canonical tables for p Boolean functions. 





Upper limit for Lp(n). Let Lg(n) be the minimum number of contacts sufficient for the realization by the 
canonical method of an arbitrary Boolean function in n arguments. Let us consider the canonical table for the Boolean 
function f(xy, Xg, . .- » Xn). In the synthesis of the canonical table, from a set of identical indices in a given row a 
single index is joined to two indices of the following row. Therefore the number of contacts is equal approximately 
to 2R¢(n), where R¢(n) is the number of distinct indices of the canonical table, but does not exceed this quantity. 


Let us estimate R(n) = sup R¢(n) over the set of all Boolean functions of n variables. We consider the (n—m)th 
row of the table. In this row there are 2"—™ indices. Therefore the number of different indices R(m, n) in the 
(n—m)th row is defined by a 2™-place binary number in the top row. From this it follows that R(m, n) also does 
not exceed 22™) the number of distinct 2™ -place binary numbers. Thus, R(m, n)<2"-™, R(m, n)< qt 


It is evident that the number of distinct indices R(n), for arbitrary whole m,, satisfies the inequality 


n ™, ® 
© 1 » SY ars 1 9n—m 
R (n) = >; Rim, n)< 2° : >» ~ ° 


m=1 m=1 m~1 


In the first sum the terms increase as the square. Therefore, the first sum is less than twice the largest term. 
Consequently, 


R(n)y< 2.227" 5 on—m, 


Taking m, = [loge(n—a@ log,n)], where @ > 1, we obtain 











yn y Aa 
R(n) <2— - i . 
n = (an — a loge ») 
Hence 
gn+2 9n+2 
L ro Ser he 
6("< n—alogen * 
As n— , we have 
g*t2 
L("< = -(1+e) (e-+0), 


In the case of a (p, 1)-pole we have the inequality: 


Rim, n) << 27", Rim, 2) < p2"—", 


from which 
R (ny < 9.92" 4 pa*-™ 
Taking n = [logs {n + logs p — a loge (n + log: p)}| , we obtain, as n-> « 
gn+2 


Le (Pp, ny< P 


n+ lop re) (e@-—+ 0), 































For linear s = asx, i anil ant", aj = 0, it is possible to give a better estimate R(m, n) = pSp_m, in 
place of the previous R(m, n) < p2"~"". As far asconcerns the estimate R(m, n) < 2” , for slowly increasing co- 
efficients aj we can substitute the improved estimate R(m,n) < g5n—Sn—m_ 


Estimate of the mean number of contacts. Let Hp, ; denote the class of functions of n variables taking on r 
values equal to 1, and (2"—r) values equal to 0. It can be assumed that u = r/2" is the probability that a given 
constituent enters into the expansion of the function. 





Let us divide the row of 2" values of the Boolean function into 2"~™ parts with 2™ binary values in each 


part. The probability of two such parts coinciding is equal to P,, = (u? + v2", where v= 1-u. Consequently, 
in the (n—m)th row. We shall carry out the calculation as follows. To each index in the row we assign 1 if all in- 
dices in the row to the left are distinct from it; in the contrary case we assign 0. Both of these values will have a 
certain probability. For the first index the probability of 1 is equal to 1. For the second index the probability of 1 
is equal to Qr = 1—Py, for the third index G,,, etc. Therefore, the mathematical expectation for the number of 
distinct indices in the (n—m)th row will be 


M (n,m) =1-+Q), 4.2.5 g2” ~m -- 1 
m =a a. . 


Thus, in the (n—m)th row the mean number of distinct indices is bounded by the number a . With 
(i +v 
increasing m this number increases. 


But N (n, 1.1) does not exceed the total number of indices in the (n—m)th row. We thus obtain for M(n, m) the 
two inequalities: 





M (a, m) < oat. M (n, m) rat _, 
(u? 4 22)? 


For the mean number of distinct indices M(n) we have, for arbitrary whole m,, the inequality 


my, 


n n 
M (n) = >) M (n,m) < aaa -- > io 
m=1 m=! (u? + v2)? m,+1 


From this it follows that 


M (n)<2 inedige f. gu my 


- 2°, Taking m; = [Hogs (4 ~ ; logs 5 


ut + 2 





«pnt on 
M (n) Bess t: rg ge B-2 
n 1 


5) _ (n —a log, +) 





Multiplying by two, we obtain the estimate for the mean number of contacts 


n+2 n+2 
Ne le + one 


(+) n—a loge 7 








As n-—> @, we have 
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where € 0, if 


gnt2 
C (n) < ~~ loge (u? -|- r?) oe (1 -!- e). 


Continuing these consideration, we obtain the estimate for the (p, 1)-pole as n + « 


p-2"t 
C (p, n) << — logs (u? + vr?) a+ Inu (1 +-e). 
, 7 
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Current types of automatic relay systems for the statistical quality control of serial production 
are considered. Methods are given for finding the limiting distributions of the regulated para~- 
meter occurring with the use of these systems under conditions of linear time variation of a 
deterministic perturbation. The possibility of optimal choice of parameters in these systems 
is demonstrated. 


The quality of serially produced objects usually depends on a large number of physicochemical and techno~ 
logical factors having differing degrees of variability in the process of production. The variation of these factors 
causes deviation of the product quality from the prescribed requirements and special measures are required to 
eliminate these variations, termed further corrective adjustments of the equipment. 


In many serial production processes gauging of the results of a given technical operation on a given sample 
(part) can be realized only after the operation on the given part is completed. Under these conditions the correc- 
tive adjustments to the equipment to obtain the necessary quality of the given part in a given series can be realized 
on the basis of the results of gauging previously finished parts of the same series.° The possibility of realizing such 


* By series we understand the ensemble of identical parts fabricated one after other on the same installation in the 
absence of explicit changes in the production conditions (for example, without change in the quality of raw ma- 
terials, without change of gauging instruments, etc.). 














a corrective adjustment is based on the fact that in many production processes, together with random quality fluc- 
tuations, characteristic only of a given sample, there exist quality variations constituting a regular (deterministic) 
function of n, where n is the ordinal number of the part in the given series. The presence of such regular components 
is caused, for example, by progressive wear of the equipment, slow temperature variations, etc. 


The magnitude of the regular component cannot usually be given in advance, but its value for a given nth 
detail can be estimated on the basis of a statistical treatment of the measurements on the preceding parts, belong- 
ing to the same series, On the basis of such an estimate it is possible to realize corrective adjustment of the equip- 
ment for the purpose of reducing the quality deviations of the nth part from the prescribed nominal value. Follow- 
ing [1], such a procedure is termed statistical quality control, widely used in mass production (cf., for example, 
[2-4]). However, in these applications, statistical quality control, constituting essentially a nonautomated process, 
has its principle purpose not in the maximum stabilization of production quality but in reduction of the work involved 
in sampling quality control. Rising requirements on manufacturing precision have made it necessary to establish 
automatic statistical quality control systems. The first examples of such systems appearing in the recent years have 
been employed principally for the corrective adjustment of metal-working machines [5-11], as well as for other 


production processes [12, 13]. 


The overwhelming majority of currently known automatic statistical quality control systems are in principle 
of a single type, which we shall term relay-type statistical automata (RTSA). 


1. Relay-Type Statistical Automata 

To concretize the further discussion, let us consider RTSA employed for 
the corrective adjustment of a certain technological aggregate, with the aim 
of stabilizing the geometrical dimensions of the parts produced on it. 





Let the controlled dimension of the given (nth) part in a given series have 
a derivative Y(n) from the prescribed nominal value, which may be represented 
in the form 


¥ (n) = X (n) + re(n) =r (n) + 8(n) + ra (n), on 





Fig. 1. 


where X(n) is the deviation of the controlled parameter in the absence of correc- 
tive adjustments, r(n) is the regular (deterministic) component of X(n), 6 (n) is 
the random component of X(n), r,(n) is the correction introduced by the adjustive mechanism. 


Since the quantity Y(n) can be measured only after processing of the nth part is completed, the process of 
stabilizing the controlled dimension about the value Y = 0 consists in an estimate of the expected value of Y(n) on 
the basis of the measured values Y : y(n—1), y(n—2), . . ., and introduction of the corrective adjustment rj,,(n), in 
the direction of decreasing the expected value of Y(n). In RTSA this process is carried out in the following way. 


In a certain neighborhood of the point Y = 0, corresponding to the nominal value of the regulated parameter 
(dimension), a zone of admissible positions of the center of scattering (i.e., of permissible values of the mathematical 
expectation Y) is defined. The position of the center of scattering with respect to its permissible zone is estimated 
by a special instrument—an analyzer (Fig. 1). The results of measuring already finished parts are applied to the in- 
put of the analyzer A from the gauge G. The analyzer distinguishes three states of the center of scattering: positions 
within the zone, and excursions beyond its positive and negative limits. In this regard, after the part with number n 
is gauged, a signal &(n) appears at the analyzer output which can have one of three values: 0, +1, -1. The signal 
#(n) = 0 corresponds to the absence of the need for corrective adjustment [r,(n + 1) =r, (n)]. The signals +1 and -1 
cause single corrective adjustment steps rj,(n + 1) = r,(n) + a or r,(n + 1) = r,(n)—a, respectively, where the absolute 
value of the corrective step a is a given constant. Thus, the corrective effect of the adjustment mechanism C is 
equal to 


Tp (n) = ay ® (7). 7" 


j=1 
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The signal (n) is a function of a certain set of magnitudes ¢(j) G = n, n—1, n—2, .. .), where ¢(j) is a func- 
tion of the comparison of the. dimension of the jth detail with the control limits, established at preassigned distances 
+b from the prescribed nominal value: 


0 for ly ()| <8, (3) 
p(s) ==: | for y(i)< —b, 
“" for yl/)>b. 


Depending on the method of determining (n), i.e., on the method of estimating the position of the center of 
scattering, existing RTSA can be classified among three types. 


I. In the simplest case, where the range of production tolerances greatly exceeds the range of scattering of 

the random deviations, the quantity b may be taken sufficiently large so that, with coincidence of the center of 
scattering and the prescribed nominal value, the probability of the 
dimiensions of individual parts exceeding the negative (I) or positive 
(II) control limit of the zone be practically equal to zero (Fig. 2). From 
this it follows that the first measurement occurring beyond the control 

i limits of the zone should be accompanied by a corrective adjustment, 
i.e., the equality 


® (n) = 9 (n) a 


should be satisfied. 








-b *b The magnitude of the corrective step a in this case can be taken 
y fairly large b = a = b—d, where d is the half-width of the random 


variation scattering field (cf. Fig. 2). 











The above conditions are characteristic, for example, for systems 
used for the control of centerless grinders [8, 9] and in the experimental 
machine ENIMS [5]. 


For higher requirements on the precision of the adjustment, i.e., with decrease in the width of the tolerance 
field to a value close to the width of the random scattering field 2d, the control limits are established within the 
random scattering field b < d (Fig. 2). Here, obviously, even with coincidence of the center of scattering with the 
nominal value, i.e., with total compensation of the regular component, there will exist a certain probability of the 
deviation of the dimensions of individual parts beyond the limits of the control limits, i.e., a probability of in- 
correct adjustment, leading to increase of the scatter in the stabilized dimension. 


It is possible to decrease the influence of incorrect adjustments either by a suitable, choice of magnitude of 
the adjustment step, respecting relationship (4) (system with “small-pulse” adjustment [5]), or by decreasing the 
probability of incorrect adjustments by using systems of types II and II. 


Il. In type II RTSA the corrective adjustment is carried out only when a certain number m of parts in succes- 
sion exceed the same control limit, i.e., 


( a 


for PU) = ™, 


j=n—m-+1 
n 


P(n) = i for > e)=—m, (5) 


fon—m+1 


0 
























Here, if |#(n)| = 1, after the corrective adjustment is completed, the nth part is not measured, i.e., it is 
assumed that ¢(n)= 0. 


Examples of such systems are given in [7] and [11], where m = 2. 


Ill. In type Ill RTSA a regular corrective adjustment is carried out in-those cases where, among the M last 
parts produced since the previous corrective adjustment, the difference between the numbers of times that the nega- 
tive and positive control limits are exceeded reaches or exceeds m, where M and m are preassigned numbers. In 
other words, if the last corrective adjustment was carried out after measurement of the ith part, with n > i we have: 


forn2=i+M 
P (n) = ( signG/(n), if \G(n); > m, r 
iy if |G n)! <m, 
forn<i+M 
(p(n) = 0, 
where 
Gt) y (/). 


joen—M-+1 


An example of this type of RTSA is given in [12]. 


We remark that type I and II systems may be considered as particular cases of the type III system, obtained 
with M = n and M = m = 1, respectively. 


2. Investigation of RTSA Operation in Stabilization of the Regulated Parameter with 
Linearly Increasing Perturbation 






















Let us consider operation of RTSA within the limits of a series consisting of a certain number of parts n,. We 
shall terms the instant of termination of the nth part the nth instant of time (n = 1, 2,... , ny), and we adopt the 
following assumption. 


A. The regular component X(n) varies linearly, i.e., 
X (n) = en +- b(n), (7) 


where c is constant for the given series. 


Equation (7) may be conside,ed characteristic for the operation of a series of metal-cutting machines, where 
the term cn is due to the progressive wear of the cutting instrument. 


B. The distribution of random deviations is independent of n, and there is no correlation between the random 
deviations of two different parts. 


Obviously, the mathematical expectation of the random deviations is equal to zero, since in the contrary case 
it would have to be taken into account in the regular component. 


C. The integral distribution of random deviations is given by a piecewise -continuous function F(§). Here 
there exists a magnitude 2d (width of scattering range of the random deviations), such that 


F(—d)=0, F(d)=1, O<F()<1 for [d8|<d. (8) 


D. The quantities c and d characterizing the perturbation X(n), and the quantities a and b, characterizing 
the parameters of the automaton, are given whole numbers. For rational a, b, c and d this can always be obtained 
by a suitable choice of scale for the regulated parameter Y. 
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We shall demonstrate that with these assumptions the distribution of the set of quantities Y(n) (n = 1, 2,...,n,) 
tends as n, -> © to the function (y), which may be calculated on the basis of prescribed a, b, c and F(6). We 
shall carry out the proof for the type III RTSA with prescribed values of M and m. At the same type all relations 
obtained will be valid for type I and I RTSA with the substitutions M = m and M = m = 1, respectively. 


Considering (1), (2) and (7), we find the mathematical expectation of the quantity Y(n): 


n—!1 
yy (n) = MY (n)] = cn +a >» (7) =cn+ an’, (9) 


j=1 


where n" is equal to the difference between the numbers of positive and negative corrective adjustments carried out 
up to the time n(|n'] < n). 


The quantity yo(n) will be termed the position of the scattering center at the nth instant of time. 
From (9) it follows that all possible values of yg(n) must be whole numbers, multiples of a, where a is the 
largest common divisor of a and c, such that 


a=aA and c= ac, (10) 


where A and C are mutually prime whole numbers. 


In order that lim y9(n) # © , the magnitude of the corrective step should satisfy the relationship 


n+ 
a 2 Me, (11) 
due to the fact that, according to (6) and (9), during the interval between two corrective adjustments the scattering 


center shifts by a value not greater than Mc. 


Where (11) is satisfied, which shall be assumed everywhere below, the possible values of yg(n) are bounded: 
Ymin = Yon) = Ymax(n=1,2,...,0), where Y¥,,;, and Yay are finite quantities whose values can be found 
from the given process parameters (cf. Appendix I). Thus, the number of all possible positions of the scattering 
center in the given system is equal to 


N= *max —*min , 1. 7 


The coordinates of the N possible positions along the Y axis are defined by the quantity 


Y (s) = Mnin + (s — 1), (13) 
where s is the ordinal number of the position (s = 1, 2, ..., N). 
If it be given that y(n) = ¥(s), then, on the basis of (8) and (9), the probability of all possible values of ¢(n) 
and yo(n + 1) are defined by the system of equations (14) and (15), respectively: 
P(g(n) = + 1] = PI6< —b—F¥ (s)] = Fi — 6 — Yl, 
P(g (n) =0)= P(—b— ¥ (s) C8 Cb—¥ (8) = 
= F(b—Y (s) + 0} — F(—b—Y (s)}, 
P (g(n) = —1) = P{6>b—Y(s)]=1— F(b—Y (s) +0) 


(14) 


and 


P [yy (n + 1) = ¥ (6) + c] = P(@(n) = 0), 
P (y(n + 1) = Y (8) 4+ e+ al = P(M(n) = 1}, 
P [yg (n + 1) = Y (8) 4+, ¢—a]=P[O(n)= —4). 











The state of the system at the instant n will be characterized by a vector with components y9(n) and t(n), where 
t(n) is defined by the difference between given instants of time n and the instant of carrying out the last corrective 
adjustment i(i > n): 


i(n)=n—t : {<n—igM -+-1, 
(16) 
i (n) = M -+ l for n— i> M =“ ‘. 


Thus there exist in all k = N(M + 1) possible states, to each of which we assign its ordinal number h(h = 1 to k). 


Let the state of the regulated process at the instant n be given: t(n)= T, yo(n) = Y(s). Then the right-hand 
sides of Eq. (15) are independent of n and can be found as functions of Y(s), T and the prescribed process character - 
istics. Considering (15) and (16), and the fact that corrective adjustments do not occur in the RTSA for T < M, we 
obtain 

for T< M 

Plyo(n +1) =Y(s) +e, t(n+1)=M +4] =1, - 
for T=M 


P [yo (n + 1)== Y (s)-+-c, t(n-+ 1) = M + 1] =2,(s, T), (18) 
P [yo (n + 1) = Y (s) +e-+a, t(n-+1)=1] =2,(s, 7), 
P [yo(n + 1) = Y (s) + c—a, t(n +1) = —1] =%(s, 7). 


Here 2s(s, T) + z2(s, T) + 2(s, T) = 1, and, consequently, no other states except those given by (18) can exist. 
The method for determining the right-hand sides of the system (18) is given in Appendix II. 


On the basis of (17) and (18) a quantity Pen can be found for arbitrary pairs of possible system states with 
numbers h = € and h= n, independent of n. This quantity is equal to the probability that the system existing at a 
certain instant of time n in the state h = & pass in the following (n + 1)thinstant of time to the state h= 1. Thus, 
the process of variation of the states of the RTSA forms a simple homogeneous Markov chain with a finite number 
of states k and discrete time. The law of formation of this chain [14] is defined by the matrix 


PuPi2. . - Pax 
(Pex) -= PaP22 ar a Pek as) 
PriPra -+ ~. Prk 


It can be shown that the matrix (19) has no isolated subgroups of states, i.e., that all states existing in the 
given chain can be found in a single nondecomposable submatrix of matrix (19). 


Let observations be carried out on the process satisfying (19) in the time interval from 0 to n,. Then, as is 
well known, the ratio of the number of times the process is incident on a certain state h(h = 1 to k) to the total 
number of observation nj, tends as nj, - © to the limit f},, independently of the initial state of the system. The 
quantity f,, can be calculated on the basis of the given matrix (19). 


Let us denote the characteristic determinant of matrix (19) as 
P(i)=[(E— Pj, (20) 


where E is the unit matrix of order k, while p is matrix (19). 


The principal minor of determinant p(1), obtained by suppressing the line and column with index i, we denote 
as Pi; (1), and the above-mentioned quantity fj, is calculated from the relationship [14] 
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Fig. 3a. c/s =0.25; I—a/o=0.5, 2—a/s = 0.75, 3— ajo = 4 
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. 3b 


‘c/s = 0.5; 1—afs = 0.75, 2—a/s =1, 3—a/s = 15, 
4—ajs=2, 5—a's=2.5 
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Fig. Sc. c/s = 1; 1— af = 1.5, 2—a/o=2, 3—aj6=3, 4—a/o=5 











Fig. 3d. c/s = 2; 1—a/s = 2.5, 2—a/o=3, 3—afo=4, 4—alo=6 





LV (21) 


where Pyy,(1) = Pig(1) where i = h. 


It is obvious that for sufficiently large nj, the quantity f,, may be considered as the probability of finding the 
system in the state h. Summing all fj, relating to a state with a given value Y(s) (summing over all T from 1 to 
M + 1), we can find the quantity f,(s = 1 to N)—the probability of finding the scattering center at the point with 











(21) 






























































ela V My —M [y))) | My) Yu es 4. E V My —M [y]))|M [y} VM iY) 

a o oa a fe] a 6 6 oa oa 
0.25 10.25 1.00 3.25 | 3.40 11.00)1.00) 1.00 4.00| 4.413 
0.25 |0.50 1.13 0.85 | 1.42 4.00)1.25 1.20 2.25) 2.56 
0.25 10.75 1.16 0.60 | 1.30 [4.00)1.50 1.25 2.05} 2.40 
0.25 |1.00 1.20 0.53 | 1.34 (4.00/2.00 1.37 2.00} 2.43 
0.25 11.25 1.23 0.50 | 1.33 14.00)3.00 1.58 2.35| 2.84 
0.50 |0.50 1.00 3.50 | 3.65 #1.00/5.00 1.92 3.30} 3.81 
0.50 10.75 1.15 1.35 | 41.77 }2.00)2.00 1.00 5.00} 5.40 
0.50 |4.00 1.20 1.12 | 1.64 {2.00)2.50 1.45 3.70} 3.98 
0.50 |41.50 1.27 1.05 | 1.65 12.00/3.00 1.59 3.65} 3.99 
0.50 |2.00 1.35 1.15 | 1.77 }2.00/4.00 1.76 4.00} 4.37 
0.50 |2.50 1.42 1.35 | 1.92 P: 6.00 2.14 5.00; 5.45 

coordinate Y(é). On the basis of the values of f, obtained we find 
MP) the limiting (for sufficiently large n),) distribution law for the regul- 


lated parameter Y 


N 
¥ (y) = >) /.F ty — Y(s)}. (22) 
s=1 


From relations (8), (13), (22) and the piecewise continuity of 
the function F(6), it follows that ¥(y) is piecewise continuous and 
¥(Ymin-4) = 0, ¥Y%max+d)= 1. From this the mean-square value 
of Y is defined as 











i i = 
. oe ee er 2 ee ¥max +4 
M |Y*| = \ yt" (y) dy, (23) 
Fig. 4.4. 1-— eo = 2, 2—e/o = 1, ade 
3—c/o=0.5, 4— e/6 = 0.25 in 


where ¥'(y) is the generalized derivative of ¥(y) [15]. 


Since with variation of the RTSA parameters the shape of the curve ¥(y) can vary, the quantity M[Y *) is the 
most convenient characteristic for comparing the qualities of regulation obtained for various values of the RTSA 
parameters. 


From obvious physical considerations we shall assume the continuous dependence of the function ¥(y) on the 
system parameters: For given quantities M and m and arbitrary fixed y, the increment to ¥(y) due to variation in 
the parameters a, b and c can be arbitrarily small, if only these variations are sufficiently small. The validity of 
this hypothesis is confirmed by numerical examples presented below. 


Since ¥(y) depends continuously on the parameters, the quantity M[Y 24, as is easily proven on the basis of (23), 
is also a continuous function of the same parameters. Starting from this, it is possible to construct, by graphical- 
analytic methods, the curve M[Y?} as a function of the system parameters, and, on the basis of this curve, find the 
optimal parameters providing minimum mean-square Y. An example of such a construction is given below. 


We note that the continuity hypothesis permits extending the results obtained to a more general case of varia- 
tion of X(n) than that given by (7): 


X (n) = cy + cn + b(n), a4) 
where cy is a given constant. 


In effect, let ¥(y) and ¥;(y) be the limiting distributions of Y with cp = 0 and cy 0, respectively, and the 
system parameters a and c correspond to (10). Then, if cg is a multiple of a, c( coincides with one of the possible 
positions of the scattering center obtained for cy = 0. 








Consequently, the behavior of the system for X(n) given by Eq. (7) and Eq. (24) differs only in the initial states 
which, as we know, have no influence on the limiting distribution and, therefore, ¥(y) =¥,(y). 


Let now Cy not be a multiple of a. We then consider a new system, differing from the previous one only in 
that the quantity a is substituted by a, = a + 4a. It is easily shown that the quantity Aa can be taken such that cy 
coincide with one of the possible positions of the scattering center in this new system, while at the same time Aa 
will be less than an arbitrarily small positive number ¢. From this, on the basis of the continuity hypothesis, we 
come to the conclusion that for arbitrary cp the limiting distribution of Y coincides with ¥(y). 


3. The Operation of Type I RTSA 





Type I RTSA may be studied on the basis of the above general relations, substituting M =m = 1. The general 
relations can be substantially simplified, since the right-hand sides of (18), as is easily seen, do not depend on T in 
the present case. Thus, the state of the system at the time n is fully characterized by a single quantity —the position 
of the scattering center yo(n). We ascribe the index s to this state, corresponding to the value of s in the relationship 
yo(n) = Y(s) (s = 1 to N). 


Let at the instant of time n the system be in the state with index s = €. Then, on the basis of (4), (14), and 
(15) we have 


i. e\ 


yon + 1)= ¥(84+5=*)|=F|—o—¥ @ |, 
Plyo(n+1)=¥ (+ -<)|=Fib—¥ &) +01—F (—0—-V@. se 


Plyo(n +1) =¥ (§ + *S*)|=1—Fb—-F @+ 01. 


P 











On the basis of (25) we calculate the probabilities of passage Pen (E, n = 1. to N), constituting the elements 
of the matrix (19), having in this case order k = N. 


The limiting probabilities f, in this case coincide with f}, and are calculated directly from (21). 


We present below the results of applying this method to a type I RTSA with normal distribution of the random 
deviations, approximated by the function F(6): 


5 5 


| Tae ’ + Xn 
F (6) = oVx Pe ds for_=— [8| << 3s, (26) 
F (— 3s) = 0, F (35) = 1. 


The limiting differential distributions of Y calculated* for various values of the parameters a and c, with the 
condition b = a, are given in Fig. 3a - 3d. 


The unit of measurement for all quantities was taken as the root-mean-square deviation of the random com- 
ponent o. As is evident from Fig. 3a-3d, the curves V¥(y) for relatively small values of a/o are very close to the 
norinal distribution, With increase of a/o the shape of the curves ¥'(y) gradually changes. 


On the basis of the curves of ¥'(y) found, values of the quantities VM[Y2], were calculated and tabulated, 
and used for the plotting of the curve family in Fig. 4. The curves terminate at the left at the locus of minimum 
admissible values, whose equation 





od = Via+ av + 3 (27) 





is obtained from the fairly obvious circumstance that with a= c the curve ¥(y) is obtained from the curve F(8) 


The quantities f were calculated on an electronic computer. 
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[Eq. (26)] by shifting the latter along the axis of abscissae by the magnitude a + 30, The presence of the line (27) 
simplifies interpolation of the curve family of Fig. 4 for intermediate values of c. 


The curve family in Fig. 4 can be used both for analysis of the system operating quality with given a, and for 
choice of optimum values of a under given operating conditions. 


Let, for example, the quantities c and o, constant in each given series, vary independently of each other from 
series to series in the limits: 6, <6 < 26,,~— 1,06, Cc < 6,, where oy is a given constant, It is required to find 
the value of a giving a minimax of the mean-square of Y, i.e., a minimum value of M[ Y*} for the worst series. 
From an examination of the curves of Fig. 4, it is evident that other conditions being equal, M[ Y*} increases with 
increase of |c| and o. Thus, the worst series corresponds to o = 20; and |c| =1.50,;here |e|/ oo = 0.75. 
From a consideration of Fig. 4 we conclude that the minimum M[Y*]with c/o = 0,75 corresponds to a/o = 1.3, 
Thus, agpt = 1.30 = 2.60, and the value of YM [Y*) obtained is equal to 404. 


At the same time, / M [Y?] <4o, for given c and a, lying in the above limits, 


APPENDIX 
Determination of the Maximum and Minimum Possible Positions of the Scattering 








Center Y¥ and Yuin 
By Ymax and Ymin we understand effective positions of the scattering center, i.e., such as can be reached 
from any position Y(s) (Yin < Y(s) < Yynqx) in a finite time segment. 


max 





We denote by Y, and Yz the maximum values Y(s), satisfying the conditions: 


P (p(n) = 0/ yo(n) = ¥(s)] #0, (28) 
P(g(n) = 1/4 (") = ¥ (s)] #9, (29) 
respectively. 
By Ys; we denote the minimum value of Y(s), satisfying the condition 
P (p(n) = —1/ y(n) = ¥ (s)] #9. (30) 
Let d+b=sa+ 4, d— b= sea + 2, 
(31) 
where $1, %, €y, €2 are integers, and 1 =e,e = @ 


Now, on the basis of (8) and (14) we obtain 


Y, = 92 =d+b— 4, VY, = 991 = d— bd — e2, Y; = —s2 = b —d + eg. (32) 


The maximum value of Y(s) reached from the position Y; up to intervention of a negative corrective ad- 
justment will be denoted as Ymax ,. On the basis of (6), (9), (16) and (28), we obtain 


a m (33) 
Yo (") = Ymaxy = Yi + Me for it(n)=M, 


The maximum value of Y(s) reached by a positive corrective adjustment will be denoted as Ymax. Obviously 
this will be reached at a certain time instant n = i + M + 1, such that 


i+m 
Yo (i +m) = Yo, i(it+m)=m, > y ({) = m. 


j=i—-1 





My (a) = a = Y,-1-(M 4-1 --m)ec+ a. oa 


Taking into account that b = a always, on the basis of (32) - (34), we obtain Ymax1> Ymaxz» and therefore 


st ) a sd -+- b -+{- Me- ~&. (35) 


The position Yin reached by a negative corrective adjustment at the time instant n = i + m is such that: 
n—t 


M% (i) == Va, i (i) <= M--m-+- I, Y' @ (J) => me. 
i on—Al 


Yy () == Ymin Y,-{.me--a for t(u) =. 1. (36) 


From this, taking into account (32), we obtain 


Ymin oa - ae d +- b +- mC -;- Cs. (37) 


APPENDIX II 


Method of Calculating the Probabilities of Passage for the General Case (Type III 
RTSA) 








Let us consider the behavior of the system at the time n such that the last corrective adjustment occurred at 
the time if(i < n), while there were no adjustments at the times j(i < j =n). 


Let yg(n) = ¥(s), then, taking into account (9), we obtain 


me 38 
Ye (i) = Ye(n-- t) = ¥ (#) — te, ” 
where ~ =n-j (i<j =n). 
Let us introduce the following notation: 
P;= Pie(/) = 4), 9;= Pig(i)=—1], r5= P(e i) = 9) = 1— p;—9;- @9) 
On the basis of (15) and (38) we can calculate 
P;= Pils, t), q;= 1s, t}, r,=rl[s, T}. (39a) 


We find the probability Q(v, w), such that in the ensemble of M parts with indices j(n—M + 1 = j =n) there 
are just v parts with g = 1 and just w parts with ¢ = -1 


Q (v, w) = >) I] a. Ga,-ee, a,,) >> Ile. Be, *@-, B..) A (a, B). (40) 
a B(a) 

Here II] (ay, a, ..., 4»), is the product v of numbers pj with indices j = a, %,..., %ys representing one of 
the (Mw) possible combinations of all possible values j (n—M + 1<j <n); II (1, Bs, -.-, By) for given values oy 
Gp, .. . , Gy, is the product w of numbers q; with indices j= By, Bz,..., By» tepresenting one of the (,,” .) 
possible combinations of all j, except j = ay, %,..., &y; 


A(@, B) =A (oy, ae,..., @,, Br, Be---+ Buy) = Py" x °° "+(M—r—w)’ 
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37) 
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where ¥3, Ya+--» WM —v— wycumecpons to all possible values of j except j = Gy, Mg, . «5 Hr Byy Barone By 
+ signifies summation over all ( mM) combinations of a; re fe) denotes summation over all (,,%,) combinations of 


B, possible for given a, @, ..., Gy. 


Ce eee) See EE w ae M! terms of the form 
M/\M—v/ > viw! (M —w— v)! 





TT (a1, a2,..., a5) T1 (81, Bes -- ++ By) (2, B)e 
Let us find the probability that the quantity G(n), defined according to (6), is equal to an integer L(0 =|L | =M): 


R(L) = P[G(n) =L). 


It is obvious that the probability R(L) is equal to the sum Q(v, w), taken for all values of w and y satisfying the 
conditions: 


v—w=L, w+v<gM. 





(41) 
From this, with L > 0 
Vmax 
R(L)= >) Q(v, »—L), (42a) 
t=L 
M+L 
where Vinx is the largest integer = M 
With L< 0 
—Wmax 
R(L)= >) Q(w—iL}, w), (42b) 
v=], 


where Wmax = Vmax: 


On the basis of (39a) and (40) the quantity R(L) can be represented in the form of an explicit function of s, 
independent of n 
R(L) = RL, 8). 


On the basis of (6), for n > i+ M, we have 


| G (n) — G (n—1)| =| 9 () — 9 (n— M)| <2. 

From this it follows that, withn > i+M 
|G(n)|<m_ = for = G (n) p(n — M) BO, 
|G(a)|qm+t for &(=)e(n—M) <9, 


since in the opposite case the corrective adjustment would have occurred at the instant j, i+ M sj < n. 
Forn=i+M 
|G(n)| QM. 





in (18) (T = M and T = M + 1). 





Taking into account (6), (42) and (43), we find the values of the probabilities of passage z)(s, T) and z(s, T) 


With T = M(n = i+ M) 


N 
t2(s,M) = P[G(n)>m]= >) R(L, 8), (44a) 
L=m 
—M 
t3(s, M)= P(G(n)q—m]= ) RIL, 5). 
L=—m (44b) 
With T=M+1(n5i+M) 
#2 (8, M +1) = P [G(n) =m] + 9,_y P[G (a) = m+ 1) = | (45a) 
= Rim, s)+q(s, M) R(m +1, 8), 
3(s, M +1) = P(G(n) = —m) + p,_y P[G(a) = —m—1) = 
= R(—m, s)+ p(s, M) Ri—m —1, s}. (45b) 


~ 
om 
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THE QUANTIZATION ERROR IN AUTOMATIC CONTROL 
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A method for determining the required parameter measurement frequency in discrete automatic 
control is considered. A series of assigned curves describing the variation of the controlled para - 
meters and the assigned errors of the time quantization of the function serve as the initial calcu- 
lation data. One of the possible methods for reducing the time quantization error is described. 


In the practice of the automatic control of industrial processes, there are often cases where f(t, + hi) values, 
where i= 0, 1, 2,3, ..., which are arbitrarily taken at certain time intervals h; = const, are measured or recorded 
instead of the total function variation in time. The advisability of such time discretization can be explained by the 
striving for a manifold utilization of individual elements in the apparatus or of the coupling channel in cyclic inter- 
rogation of the parameters as well as by a tendency to use digital recording. 


The replacement of a continuous function by a time-series of instantaneous values leads to irreversible losses 
Al in the information value of the measurements, and, consequently, to an additional error 8, [1]. With a greater 
frequency of readings, this error will correspondingly decrease. However, if the measurements are made too often, 
the possibilities of increasing the time-duty of the equipment are limited, and, for many element types, this leads 
to a reduction of their lifetime. On the other hand, if the readings are not taken often enough, the error 6; attains 
inadmissible values. Hence the necessity of finding a mathematical relationship between the number of readings 
per unit time and the time quantization error. 


In a number of papers concerned with discrete measurements (see, for instance, [3]), linear interpolation 
(Newton's formulas) is used for estimating this error. In this, only the mutual relationship between two adjacent 
values f(tj) and f(t; + hy) of the function is taken into account. It is natural that such statement of the problem leads 
to higher interpolation error values, since the information contained in the totality of discrete values is insufficiently 
utilized. 


For a better approximation of the actual function at the points between readings, it is necessary to take into 
account several neighbouring values of the function, and, if it is possible to follow the autocorrelation function over 
the entire interval in question, to take into account all points of the given interval. 


A notable feature of a number of functions that characterize the variation of many technological parameters 
is their periodicial recurrence. It is understood that we do not consider ideal periodicity (the information obtained 
as a result of measurements would be otherwise equal to zero), but only a periodical repetition of the general charac- 
ter of changes in the function. The cause of this fact, which has been confirmed by the processing of a large number 
of curves, is the periodicity of many technological processes. The advisability of applying harmonic analysis to the 
curve describing changes in the technological process is an important consequence of this situation. 


We shall now formulate the problem thus stated. 


Let a certain arbitrary function be given in the time interval from ty to t,. We shall assume that the chosen 
interval is sufficient for estimating the behavior of the function® and that the readings are taken during time in 


* For a periodic function, the sufficient interval coincides with a single period. The determination of the period for 
actual curve does not involve difficulties if the technological cycle duration is given or if the presence of 24-hr 
(changeable) recurrence is known. 





tervals that are negligibly small in comparison with the hy intervals (Fig. 1). It is necessary to express the error in - 
curred in approximating a graphically assigned curve by a sum of terms of a Fourier series with a finite number of 
terms. It is assumed that the function satisfies the Dirichlet conditions. 
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Fig. 1. Graph of the technological parameter in variation time. 


Let us expand the function into a Fourier series with a number of terms that is generally limited by the highest 
reliably calculable finite difference of the parameter to be measured. It is natural that, in the case of such expansion, 
the function under consideration (strictly speaking, a nonperiodic function) will not be approximated by the obtained 
series beyond the boundary of the section under investigation; however, this is of no importance for subsequent con- 
siderations. 


Let us designate by m the number of series terms obtained in this expansion. If we take the spectrum of the 
curve to be analyzed is limited by the frequency Fm, we have § mp = 0, where &pp is the root-mean-square error 
in approximating a continuous curve by a sum of Fourier terms. 


In this case of such expansion, any actual curve is divided into 12, 24, 36, or 48 intervals in dependence on the 
required accuracy and the behavior of the function. 


In any case, it is not advisable to divide the curve into more than 64 to 72 intervals, since this would lead to 
very cumbersome calculations. 


The coefficients ate calculated by using the well-known Bessel functions:* 


* The complexity of such analysis must not be exaggerated. Even in the case of division into 48 intervals, only the 
following trigonometric function values figure in the calculations: 


sin 0 = cos_* = 0; sin _* = cos°™ = 0, 26; sin -* = cos 7 = . 

5 i 3 cos 3 , 26; sin 5 cos 3 0,50; 

sin * —cos-*. = 0,71; sin *_ = cos 7 = 0,87; sin °*™ = cos * =0,97; 
4 4 3 6 12 12 


~ 9a 
sin —_ = cos0 = 1. 
3 cos 


All calculations can be readily performed by mechanical computers, without mentioning the possibility of 
using harmonic analyzers, electronic evolvers, etc. If special computing devices are not available, one can use 
grouping circuits (circuits designed by Nancy Klock or M. Serebrennikov in the case of 48 ordinates) or patterns 
(for instance, A. M. Lopshits or Zipperer patterns), or one can dispense with Bessel functions and use any of the 
analytical methods, with the exception, naturally, of methods for processing obliquely symmetrical periodic curves 
or other methods the derivation of which is connected with conditions unacceptable in the case under consideration, 
In other words, the practical harmonic analysis should be performed by using any of the possible methods with which 
the designer is familiar. 
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(1) 
= 2kix 
qaj~ = >; ] (tx) cos — 
k=0 
m—1 
ae 7 1 . 2kix 
hz ma ba } (te) sin ‘me 3! 
k=0 
where m, = 12k (k= 1, 2,3, . . . ) is the number of intervals, and i= 1, 2,..., m,/2. 
Fig. 2. 
The approximating equal will be written as 
k=m k=m, 
_ & 2nt . 2mt 
Sm (t)= 3+ > a, cos —— k+ >) by sin | — k. 
k=1 k=1 
If m,> m, then apy 4 y= am +2 =--- = am, = 0 and by + y= bm +2 =... = bm, = 0. 
Therefore, 
- 2nt he 2Qxut 
f(t) Sm(t) = 2 + 5} ax cos——k + S) bx sin ——k. (2) 
k=1 7 k= ‘n 


Nothing that the above method offers possibilities for the trigonometric interpolation of the initial function at 
the points between readings, we shall estimate the root-mean-square error incurred in approximating the function 
f(t) by a sum of n terms of the Fourier series for n < m. 


By definition, we have 


tn tn 


2 i 2 i a ." _, 2nt . — 
bn ==, (d= —(VO- 7-2 a, COS k + >; On sin =— k} dt. 
0 





t 
n n n 
0 k=1 


It can be shown that 





te 


6,2 — i ’ 2 a, i : 2 2 
n= \ /Pdt—2—F> > (ak + b3). 
k=1 


0 


(3) 


Nothing that the addition of a new term will not change the values of the coefficients considered, we shall 
write: 


2 i ° 2 a, bee 2 
on41 = rs \ [7 (¢)]? dt —q-— > >) (ak + Ox). 
0 k=1 


Hence, we conclude that the addition of a single new term to the n terms that have been considered reduces 
1 
the error by 5(an +1+bh+ »- 


Consequently, it is necessary to know the higher harmonics in order to take into account the trigonometric 
interpolation errors. Since &,, = 0 according to the assigned conditions, we have: 


Sina = 5 (am + bn), 


Sra = (ain + Un) + > (Gna + bina) 
(4) 


k—m 
8, =FZ >) (ai, + bi). 
k=n+1 


By comparing expressions (3) and (4), we obtain: 


a. Bos haa 
2 | w@rat= 2 4D} (ak + bi). 


0 k=1 


The last relationship represents the closed-loop equation for a function with a bounded spectrum, which in- 
dicates that the above reasoning is correct. 


Thus, in practice, if we have a curve that characterizes changes in the technological parameter which is being 
telemetered with the time quantization error §;, this curve should be decomposed by dividing it into m, sections. 
After we determine m + 1 coefficients that are different from zero, we successively calculate 62, ,, 83, », . 
up to Sha > ot > e*. by using Eqs. (4). The latter determines in a well-defined manner the number of neces- 
sary readings during the time t, —ty according to the formula 


Nreq = 2n. 


In certain cases, the so-called average error o,, = 0.7979 § of a single measurement is limited instead of 
the root-mean-square error. 


This, naturally, does not change the analytical procedure, with the difference, however, that 6 is calculated 
up to the value 


2 > 2.04 gy > 8 


n+1 = ne 





* The equality must be understood in the sense of the accepted calculation accuracy, while the general accuracy 
estimate can be considered as the root-mean-square error ©, = ing/Vm,, » Where & ing is the error of a single 
f(t) reading (i.e., the diagram recording error), and m, is the number of sections into which the function is divided. 
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It should be noted that, in some cases, a uniform approximation instead of the approximation in the sense of 
least squares is required. 


We assume that the curve characterizing the time variation of the parameter to be measured does not have 
discontinuties (otherwise, the problem would not have a solution due to the Gibbs phenomenon). The continuity of 
the function makes it possible to state that such an approximation exists on the basis of the Weierstrass theorem. We 


h==m m 


shall try to obtain this approximation by means of a similar Fourier series. Since the On and \' },,  Seties 


are absolutely convergent, the series k= @ Ket 


me me 
\' dy COS Jatt es .t 2at 
2D; Ce Oo a k -; »: b,. sin — k 
k-o k--1 a 
are consequently uniformily convergent, and the Fourier series can provide a uniform approximation. 


For calculating the error A(t), we assume that A(t) = 0 for n= m, and 





m m 
2Qnt . om 
A(t)= 5S} axcos x mp > by sin kK for n<m. 
k=n+1 - k=n+1 “| 


Hence, 


pete m 
é i ti ‘ 2nt; 
rep ay ke sin “tk == 2 >) x k cos mi 
k=n+1 fn k =n+1 


where t; is a certain point in the interval under consideration (tg = t; = th), which corresponds to the maximum 
value of A(t). 


It can be readily deduced that, for n = m—1, 
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Further calculations lead to a system of transcendental equations, which can be reduced to high-power al- 
gebraic equations by means of the de Moivre formulas. However, the solution of such equations with respect to tj 
is very difficult, and results which would be suitable for practical calculations have not yet been obtained. 


The case of the coincidence of ordinates with the greatest error with respect to the time axis is of definite 
practical interest. It is natural that the maximum possible time discretization error occurs in this case. The con- 
dition for the coincidence of ordinates can be written as 


mb, +- (m—14)b 


m—tl 





an ma,, + (m —1) Sy—1 


or in a different form: ém-1. 4m hs < 
bn-1 om 








In this, Omax ¢ = Smax 1 + Smax2 +. --» where 


ane ’ . 
Amaxi= Y' a,cosare tg — + YS bysinarctg—™. 
k=m—i m k=m—i om 


In all other cases, it is apparently more convenient to use harmonic analyzers. Thus, for instance, the har- 
monic analyzer designed by Michelson and Straton [2] provides 80 expansion terms and makes it possible to synthesize 
a function from an arbitrary number of these terms by observing the pattern of successive approximations. 


After processing one technological parameter curve, the designer can greatly simplify this calculations in the 
remaining cases by using the following procedure. In the case where the number of readings is found with respect to 
the assigned root-mean-square error, the Fourier coefficients can be calculated only in the expected region of har- 
monics with numbers from i= ntoi=m. In the case of calculations with respect to the maximum absolute error, 
which can be reduced to the solution of a series of transcendental equations, it is most convenient to use the relaxa- 
tion method assuming the first approximation from the calculations of the previous curves. 


As a result of processing a large number of curves, we obtain functions that characterize the probability dis- 
tribution of the errors under consideration for different interrogation frequencies, which makes it possible to obtain 
the average measurement error. 


As far as it is known to the author, devices with At = const. are used almost exclusively in modern systems 
with cyclic interrogation. However, it is advisable to design a rational system in a different way. 


In most cases, there is an excess of information, which is connected with the statistical distribution of the 
derivative of the function to be measured in different time intervals. By taking into account this information, it 
is possible to reduce 6; without a superfluous increase in the parameter interrogation frequency. We shall clarify 
this by considering figure 2. 


This figure shows the curves of relative gas discharge values Q(curves 2 and 4) and pressure values p (curves 
1 and 3) at the Leningrad gas distribution station over a period of two days in winter time. It can be shown that the 
discharge and pressure curves on other winter days will lie in a very narrow region near the above curves while 
preserving approximately the same shape. Only the curves obtained on holidays will be an exception to this. 


Let us divide the entire interval under consideration into sections I to VI. We notice that the function varia- 
tion rate is different for each section. If we calculate the maximum interpolation error in each of these sections 
for an increment h; = const, we can readily see that they are different in each section. Hence, it is convenient to 
chose a separate hy value for each section. For instance, for §¢= 1.5%, we have: hy, = 10 min; hy, = 1 hr; hy, = 
= 10 min; hy, = 2 hr; hy, = 14 min; h,, = 50 min. 


The curves shown in figure 2 pertain to normal operating conditions. Under emergency conditions, it is ap- 
parently advisable to impose on the system the regime of cyclic interrogation of all units or of a certain parameter 
group with the maximum speed. Such a change in hy can be brought about by the emergency signal, which is in 
many cases provided by means of a special data transmitter. 


SUMMARY 
1. Besides errors inherent in instruments and measurement methods, an additional and unavoidable time 
quantization error arises in automatic control processes. 


The root-mean-square value §, of this error can be calculated for a given function with a bounded spectrum 
by using trigonometric interpolation according to the formula 


& = / 4 7 (a; + bi), 


i=n+1 





where a; and bj are the coefficients obtained by expanding the given function intoa Fourier series, and 2n is the 
number of discrete readings. 
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2. The above time quantization error for continuous functions can readily be calculated by means of har- 
monic analyzers and synthesizers. 


3. For periodic functions, the time quantization error can be reduced by a suitable division of the function into 
regions with a specific interrogation interval hy for each region. 
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The paper analyzes the operation of a magnetic amplifier network based on the use of even 
voltage harmonics that appear across the terminals of the output windings of the simplest 
two-core magnetic amplifier when a dc magnetizing field is present. 

The author analyzes the operation of a magnetic amplifier network based on two non- 
identical cores with an output winding connected in series with a symmetrical nonlinear re- 
sistance, 


Magnetic amplifiers based on the use of the even voltage harmonics which appear across the terminals of the 
control or output windings can be subdivided into two groups [1]. The first group includes magnetic voltage am- 
plifiers with a frequency -doubled output. They have a high sensitivity but feature a cumbersome construction and 
involve complex auxiliary equipment, in view of the necessity of connecting high-quality multiple -section band 
pass filters at the output of the amplifier and a blocking filter in series with the excitation windings. One of the 
basic reasons for the appearance of a parasitic ac voltage at the amplifier output is the variation of the character- 
istics and the violation of the identical nature of the cores; this requires the cores to be selected in the most thorough 
way possible. 


The second group includes magnetic amplifiers which are excited by means of a large ac voltage. Pulses 
whose amplitude and duration are functions of the magnetizing signal appear at the terminals of the output winding 
of such an amplifier (across the’ terminals K,—Kg in the circuit shown in Fig. 1). By selective rectification of these 
pulses it is possible to obtain a dc current component at the output of the magnetic amplifier when a signal is present 
at its input; this dc component can be applied to the next stage of a multistage magnetic amplifier. The pulses can 
be rectified by means of a nonlinear resistance R,,, with a symmetrical volt-ampere characteristic. Ideal for this 
purpose is a nonlinear resistance which has a volt-ampere characteristic such as is shown in Fig. 2. 


In this paper we assume that the symmetrical nonlinear resistance has an ideal volt-ampere characteristic 
such as that shown in Fig 2. If we form the symmetrical nonlinear resistor Rpg by connecting a silicon stabilitron 
of the type "D808-813" into the diagonal of a rectifier bridge consisting of semiconductor diodes (cf. Fig. 3), then 





its characteristic will be close to ideal. Magnetic amplifiers of the second group are distinguished by their simple 
circuity when a symmetrical nonlinear resistor is used to rectify the even voltage harmonics, and do not require 
filters for suppressing the odd harmonics in the output circuit. An experimental investigation shows that magnetic 
amplifiers of the type shown in Fig. 1 have a high null stability for changes in the supply voltage and frequency, or 
for temperature changes over wide limits [2]. However, the literature provides no analysis of the causes for the 
null drift and does not explain the physics governing the operation of the magnetic amplifier proper; quantitative 
relationships between the basic parameters of the magnetic amplifier circuit are also absent in the literature. 


The purpose of this paper is to analyze the operation of a magnetic amplifier of the type shown in Fig. 1 and 
to derive quantitative relationships between the basic parameters of the magnetic amplifier. 


Below, we shall demonstrate that the output current reaches its maximum when the control current |, is equal 
to the current value I, for which the core is saturated. Therefore, if we approximate the magnetization curves in 
the manner shown in Fig. 4 we shall lose the portion of the amplifier characteristic corresponding to the values 
I, < I, ; this is undesirable. 


We shall approximate the magnetization curve by means of three rectilinear segements as shown in Figs. 5-9. 
In that case, 


dB 
, a7 = const for |< |/,|, 
dB * 
b= FH 9 for >I Mai. 


The paper cites the results of an analysis for the case where the cores differ somewhat. The case of identical 
cores is a particular case of the more general case studied here. For the circuit of Fig. 1 the presence of a small 
difference between the cores is not only not harmful but is even 

useful from the point of view of sensitivity (1); this follows since 


























Up Sit ———>g in the case of slightly differing cores initial pulses exist across the 
terminals of the output winding when there is no signal, and these 
5 pulses keep the symmetrical nonlinear resistor “open.“ The refore 
w, Wp for the circuit in Fig. 1 we do not require the rigorous identity of 
+ A ae + - * * s . 
A * t cores which is required by magnetic voltage amplifiers that operate 
w\ “4c /4 according to the frequency doubling principle. 
{ 4 at "Hy 
4 al We assume that nonidentical cores are defined as cores in 
We which the inductances of identical windings, and (or) the saturation 
| l e fluxes, and (or) the current for which the cores saturate are not 
I, ~ Ky oF Ky identical. D,,=— Dp, La-~Lp, /ea~/]sR . In order to be 
Py R specific we shall study a sufficiently general case for which 
va Out NL 
i 
o caine 9 Op >On aM Ly < Ly. a) 
Fig. 1. The circuit diagram for the mag- Here Ip > Ia, since 


netic amplifier. 





wo ‘D 
La = Wa, a 10-8 and Le = UB Tr 10°°. 


For simplicity in notation we shall assume that all the windings have an identical number of turns: 


a ee aa, sora (2) 
Wy = WR= UC= Uy = UV. 





®, - P53 a 
Isa —_@ 10°% and IB = WwW co 10°8. (3) 
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The operation of the magnetic amplifier shall be studied for the case of very weak signals when 


I U ay 
c <.T 


where ry is the resistance in the excitation circuit, and U,, is the average value of the supply voltage. 


It should be noted that in ordinary magnetic amplifiers we are interested in the current through the supply 
circuit (which is the load circuit) as a function of the control current. For the magnetic amplifier under study in this 
paper we are interested in the current through the output winding, which is now the load circuit, rather than the cur- 
rent in the excitation circuit. Since a magnetic amplifier of the type shown in Fig. 1 is designed for amplifying 
very weak signals, it follows that the current in the excitation circuit is practically independent of the control cur- 
rent. In order to avoid any effect of hysteresis on the null drift the magnetic amplifier shown in Fig. 1 is in practice 
excited using a large ac voltage; thus, during a portion of each halfcycle both cores are in a saturated state. Under 
these conditions the magnitude of the current in the excitation circuit is limited by the resistance of that circuit. 
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Fig. 2. Idealized characteristics stabilitron 
for the nonlinear resistance. 1-9 
Fig. 3. Fig. 4. 


We shall first study the operation of the amplifier when there is no signal (Fig. 5). For t < ty both cores are in 
a state of negative saturation. At the instant t = t, when the negative exciting current i; (—) decreases to the value 
~Igg in absolute magnitude (the point 1° on the magnetization curve for core B in Fig. 5h), the core B emerges from 
the saturation state. Under these conditions the core A is still saturated (point 1 on the magnetization curve for core 
A in Fig. 5g). The instant t, is determined (when the number of turns is equal) by the condition 
(4) 


U. sin wt, 
eed = — /,pB. 
Ty 





For small resistance ry, | rylgg | «<U,, and t; #0; therefore in our case the current Ip for good permalloy is 
very small when the dimensions of the cores and the number of turns on the excitation windings are taken into account. 


For t; = t = & only the flux through core B varies; an emf is induced in the output winding, and this emf has a 
polarity opposite to the polarity of the exciting voltage (Fig. 5e). 


At the instant t = t at which the current in the excitation circuit reaches the value i;}(—) = —I,q the states of 
the cores A and B are characterized by the points 2 and 2', respectively (Figs. 5g and 5h). The core A emerges from 
the saturation state. The instant t is determined (neglecting the voltage drop rjiy) from the equation 


fs (5) 
\ Um sin wt dt = — w 10° (MD, — Dy) = w 10 Dap — Lal ca. 


ty 


Fort St=t, an emf eo = ea —ep is induced in the output winding; this emf is produced by the joint 
variations of the fluxes a and #3. The emf egy = €4 —ep is not equal to zero, since the inductance of core A is 
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greater than the inductance of core B (La > Lg). The emf eoyt has the same polarity as the voltage across the ter- 
minals of the excitation windings for the core A. 


At the instant t = t, the stages of the cores A and B are characterized by the points 3 and 3° on the corresponding 
magnetization curves. The core A saturates. The instant t, is determined by the formula 


{ 


\ © sin wt dt = w 1078 (20, ,- Og. Dy) = w 1078 (20,4 + Og) +- Lela. 


ty 


(6) 


For ts =t <t, only the flux @, changes. A pulse appears at the terminals of the output windings, since in 
that time interval the value of the exciting voltage is large and this value is transferred to the output winding by 
means of the variation of the flux ®,. The instant t, is determined from the equation 

ty 
\ U » sin wt td = w 107° (2M,, + 20 


i 


7 
wa) ™ 


At the instant t = t, the state of the core B is expressed by the point 4' and the core B saturates. For ts =t = t; 
both cores are in the positive saturation state. 


At the instant t = ts the current in the excitation circuit i(+) decreases to +1,p, and core B emerges from 
positive saturation. The instant ts is determined from the formula 


eI eal Ri (8) 


ri 


From the instant t = ts the fluxes @, and 4, begin to vary in reverse order. 


Assume for simplicity that the magnitude of the emf €,,,; induced in the output winding is greater than the 
threshold voltage Uy of the nonlinear resistance only in the case when €,,,; is in the form of a pulse. Then we ob- 
tain the following formula for the current in the output circuit: 


I erik RU 


out f 2 r 4 (9) 
: +n () 


ou Wy 


where 


“but " 
r= Tow = ry (—S = Tout 7] r;. 


when the equality of the turns w,,,, and wy, is taken into account. 
When there is no signal the average value of the current is 

T 
1 


Tr \ lout dt =. 0. 


0 


The mmf produced by the control current I,, and the excitation current i, add in core A and subtract in core B. 
When the control current I, shown in Figs. 6g and 6h is present, only the flux ©, varies at first from -®,5 (1") to the 
point 2" on the magnetization curve. 


The values of the current i; at these points will be : 
for the point 1’, 
is (—) = — (Usp — 4), 


for the point 2" 


i, (—) = (sa + k ). 





The point 2” for the core B corresponds to a value of the current i, for which the core A emerges from the 
negative saturation state (point 2). After that the joint variation of the fluxes @, and @, begins. The flux 4p varies 
from point 2" and 3", and the flux , varies from the point 2 to the value +, (point 3). After that only %p varies 
from the point 3" to +sp (the point 4"). 


A variation of the flux 4g from the point 3" to the point 4° induces a pulse in the output winding. A reverse 
variation of the flux 4, from the point 2" to the point 1' during the next negative half cycle also induces a pulse. 
Since the flux increment (@, — @,-) > (My, — @,,) and (O,,— M,-) < (M,— Dy) __, it follows that the 
duration of the pulse for t, = t = t, becomes greater when the signal is present while the duration of the pulse for 
ty St < tg becomes less (Fig. 6e). 

We shall determine the value of control current for which the duration of the pulse goes to zero for ty St <t, 
(Figs. 5e and 6e). It is evident that for the equation (Fig. 7g and 7h) 


Tea t+ Ie = 153 — Io, 
which corresponds to the formula 


To Ba “Tr 
2 


the flux 4 and g begin to change simultaneously: @, changes from- sa, and $p changes from— ,,. After the 
equation 4 = +, is reached only the flux pg varies from point 2' to the point 3' (Fig. 7h), and a pulse is induced 
in the output winding. During the next negative half-cycle the fluxes @, and 4, vary in reverse order; first only the 
flux 4p varies from the point 3° to the point 2'. After that the flux @q and 4p vary jointly and simultaneously reach 
~®sq and -®,,. Thus, during the negative halfcycle no pulse is induced in the output windings. 


For yp (153 ~ 1g) = Jos 1), (Isp + 14) two pulses of emf appear with only negative polarity (Fig. 8e). We shall 
demonstrate that for I, = */3(1g3—1,q) the output current reaches a maximum (Fig. 9). 


Whereas at first only the flux @, varies from -@,4 to +4,, and after that the flux 4, varies from -@,, to +®,5, 
it follows that the output current io, reaches a maximum as shown in Fig. 9. Under these conditions the magnitude 
of the current.i, determined from the ampere -turns equation for the point +#sa is equal to the magnitude of the 
current i, determined from the ampere-turns equation for the point -,, (Figs. 9g and 9h): 


For the point +®,, 
i, (—) - Ta 2 fae 
For the point -O.n 


i(—)=—/, +i, 
al 1 B c 


1B + Vey (11) 
— 2 ‘ 
After the control current I, reaches the value 1).(Igp + 1,4), 4 further increase in I, cannot cause an increase 


in the de component of the output current ipy;. In other words, for I, = "/2(Igp + Iq) the de component of the output 
current reaches a maximum. 


We shall derive the relationships between the dc component I,y Of the output current and the signal 1, which 
are of interest to us for the case of a weak signal I, « */2(Ign + Iga); here it is possible to assume that the time in- 
terval during which the output voltage pulse exists is sufficiently small. Then we can assume that the emf eo), is 
constant during the interval over which the pulse exists. The duration At of the pulse depends on how much the flux 
4 © varies: 


UO mAt = w 10°°AQ. (12) 
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For the average values Thos and Iheg of the output current during the positive and negative half cycles, respec - 
tively, we have 


g Ml 5 Oe (13) 
l hos = T r — Atnos 
~ f a Us 
hex: 7 ab ciehi Ateg (14) 
Substituting (12) into (13) and (14), we obtain: 
Uo 
9 U - (15) 
hos = FF WIAD 
iy 
iy Ti (16 
—_ m™ wit 
leg ae) cy Mamet 10 AP eg ) 
The dc component of the output current is equal to 
Uo 
{ ~, (17) 


» Um i 
Tay = 1 p98 — !neg= “7 ——- © 10* (AM 95 + ADyed. 


First we shall study the relationships between the control currents and the dc component of current in the output 
winding for 


O< Ie [Ye (le—n— Ja). 


As demonstrated in Fig. 5, for I, = 0 the flux variation is -A Pneg = O4'— 23", ~ A ®no5 = ;'~@, and 
A®nos + 4 Sneg = 0; therefore, on the basis of Eq. (17) the dc component of the current i,,,, is equal to zero. 


When the signal I, is present |, (0 < I¢< '/2(/s5 — /.a)) the flux ©, varies by the amount ,:—,* instead 
of by the amount ©,'— 4" during the positive half-cycle (Fig. 6h); during the negative half-cycle it varies by the 
amount ®4'— 4," instead of by the amount @4—3:. The flux increments will be 


— ADpos = (Py — By) — (Dy — Vy) = — Oy + Oy, ee 
iad ADeg = (Dy — Dy) — (Oy — Dy) = — Dy + Dy. (19) 


The flux increment over the entire cycle will be equal to 


— A®M pos - A®,., = — (0. — Dy) — (— Dy + Dy) = — 2(@, — Dy). (20) 


neg 
due to the presence of the signals. 
From Figs. 6g and 5h it is evident that 


-O,.B 
2 (Dy — Dy) w 10 = 2 [(Is4 + Le + Ie) Lg — Sealy) = 4s Tew 10°°. (21) 
Substituting (21) into (17), we obtain 
‘tion 
met ®.p ' Ia—l 
lay=—p— 7 wi0* (472) le for OSI By. (22) 
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From formula (22) it is evident that the dc component of the output current ipyr for 0 < [¢ < */g (/4B — I 44) 
is proportional to the signal I. 

We shall study the relationships between the control current and the dc component of current in the output 
winding for 

/e (1 —Ta)< cS "le (158 + Isa). 

During the positive halfcycle the flux @, varies by the amount ®,'—@,', and during the negative halfcycle 
the flux @, varies by the amount 4,— 4. 

From the magnetization curves in Figs. 8g and 8h we determine flux increments: 


— AGhjeg w 10 = (Dy — Dy) w 10 = [(Dy — Dy) — (Dy — Oy) w 10 = 





(23) 
_ (UZ, + lg)Lp — (Usa — 1.) Lg) _ Ls (27.+ (/ 5p = T,a)] 
and 
—AD jog w 10" = — (O, — D,) w 10 = — [((O, — ®,) — (D, — O,)] w10* = (24) 
- (sa + I.) La ad (1p ve [,)] La _ La (21, — (158 — T,a)). 
From (23) and (24) we find that the flux increment over the entire cycle is equal to 
(ADpos -- Ae g) w 10 = —L, (2), — (Teg — I54)] — Lp (20e + (sp —/sa)]- (25) 
It should be noted that by substituting |, = "/9(Isp —1,q) into (21) or (25) we obtain a flux increment equal to 
2 oa (1.8 —Jsa)- in both cases. 
3B 


Substituting (25) into (17), we obtain the expression for the dc component of the output current: 


Us 
4—->— 
Um (26) 
Iyy = — Ge ———™_ [(Lg — La) (lag — Sea) + 2 (Ly + La) [el 
The de component I,, is proportional to the control current I,. 


Equation (26) will not be valid if the inequality 1, << '/,(/,3 + /,,4) is not valid, since under these con- 
ditions the magnitude of the pulse can no longer be considered constant. But such an operating mode of the mag- 
netic amplifier is of no interest to us. 


For the case of identical cores we obtain 





(a 


U ®, 
lay = — Fe w10*-4 Fe. (27) 


from (22) or (26). 


From (22), (26) and (27) it is evident that a slight difference between the cores (Lg * La, Igy » I,,) has only 
a slight effect on the current gain K, = I,,/1,. 


SUMMARY 

The analysis which we have performed for a two-core magnetic amplifier shows that for the case of non- 
identical cores initial pulses appear across the terminals of the output windings when there is no signal. When the 
signal is applied the pulses decrease in one direction and increase in the other direction. 








2A) 


| to 


6) 


g- 


7) 


nly 





We have derived the basic relationships between the dc component I,,, of the output current, the control 
current I, and the magnetic amplifier parameters. 


We have demonstrated that a slight difference between the cores has very little effect on the current gain, 
and therefore no special selection of cores is required. 


For the magnetic amplifier circuit shown in Fig. 1,a small difference in the cores (and therefore a variation 
in the core characteristics) has no effect on the amplifier sensitivity. The lower sensitivity threshold of the mag- 
netic amplifier is limited to a considerable extent by the instability of the symmetrical nonlinear resistance (this 
is studied in [3]). 


The facts that pulse voltages arise across the terminals of the output winding of the amplifier and that the 
input-output characteristic is linear within certain definite limits are verified experimentally. 


The author expresses his appreciation to M. A. Rozenblat for his valuable remarks and advice. 
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A study is made of the effect of pulse noise with a logarithmically normal amplitude dis- 
tribution on a remote control unit. The probability density distribution is determined for 
the duration of the peaks at the output of the threshold unit for a low-pass filter of the RC 
type and for a filter with a frequency response corresponding to a Gaussian curve. 


INTRODUCTION 
The investigation of noise in a low-voltage power line rated at 0.4/6 kv [1], as well as other papers on the 
investigation of noise which is of industrial origin (cf. for example, [2]), have already made it possible to conclude 
that a rather broad class of industrial noise in a passband of the order of several kilocycles is of a pulse nature with 
a logarithmically normal amplitude distribution law. Therefore it is of interest to study the effect of this type of 
noise on remote control units which operate in a wide frequency band over a power line and over a radio channel. 
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Fig. 1. 
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In order to estimate the noise immunity for a number of remote control systems along with the amplitude 
distribution it is also necessary to make use of the probability distribution density for the duration of the peaks at 
the output of the receiver. A somewhat simpler problem arises in finding the magnitude of the error caused by the 
noise in pulsed telemetering systems when it is possible to limit ourselves to determining the probability that the 
peak duration exceeds a specified value; in a number of cases this problem reduces to finding the average duration 
of the pulses at a specified level. 


It is of interest to determine the probability distribution density for the duration of the pulse noise peaks at 
the output of a remote control unit (Fig. 1) consisting of a wide-band input filter 1, a linear inertialess amplitude 


detected a, a low-pass filter 3, and a threshold device 4. We shall perform our analysis for the case where the input 
filter has a frequency response such that when its input is subjected to a brief noise pulse, one pulse whose shape is 
determined entirely by the frequency response of the filter appears at its output. 


Since the distribution of the amplitude envelope for the noise is usually measured at the output of the selec- 
tive metering device, it is natural to assume that the logarithmically normal amplitude distribution of the pulses is 
determined at the output of the filter 1 and can be written as 


af (in x — Ina)* 


fief (1) 
p(inz) = Vins exp [ >? 


| (0< 2 < 0), 


where In a and ¢ respectively represent the average and mean-square values of the normal distribution for In x. 


We shall also assume that the speed of response for the system is such that the noise signals at the output of 

the receiver remain of a pulse nature (i.e., superposition of their transient responses does not occur). Such an assump- 

tion applies to multichannel and fast-response remote control systems 

when the mean repetition frequency of the pulses is sufficiently small. 
pie) For simplicity we shall assume that the input filter has a band so much 
15 wider than the output filter that the noise pulse spectrum remains con- 
stant within the limits of the passband for the filter 3. This does not 
reduce the general nature of the results, since they may be extended 
a:/ to include the case of commensurate pass bands for both filters when 
the mathematical computations are made somewhat more complex. 
For such sufficiently general assumptions the problem reduces to find- 
10 ing the distribution of the peak durations for a number of threshold 
voltage values when the filter 3 is subjected to pulses of identical 
shape governed by a logarithmically normal amplitude distribution 
law. 








a:0 1. A Low~-Pass Filter Consisting of a RC-Circuit 
As we know, the frequency response for a RC low-pass filter 
as-2 can be written as 
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C (ja) = i+ jetR’ | (2) 








where a = Aw is the passband width for the filter at the 0.7 level. 


If the input of such a filter is subjected to a brief pulse of am- 
plitude x whose spectral density is constant within the limits of the 
pass band Aw and proportional to its amplitude x, then the shape of 
the pulse at the filter output can be determined from the formula 
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Fig. 2. 


Uout (t) = CypAwaxe—4“", ” 
The proportionality coefficient Cy depends on the shape of the frequency response of the input filter and re- 
presents the numerical value of the spectral density for a pulse of unit amplitude at the output of the filter 1 at its 


turned frequency wo. 
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The time r during which the noise peak will exceed the voltage Uy, corresponding to the operating threshold 
can be determined from the relationship 

Un = CoAwzre—4*", (4) 
whence 


Aor = Inzg—InW. (5) 


Here U, = Uy,/Co A w is the value of the threshold voltage normalized to the filter input. 


It follows from formula (5) that the distribution for Awrt coincides with the distribution for In x with a shifted 
distribution center. Therefore for the probability density it is possible to write 


p (Aer) = pe 





Ke | (Aor —iIna+inU,¥ 
ee 


2c* 


| (0 < Awt < 0) 
(6) 


VU (A@t < 0). 















































Fig. 4. 


This formula represents the truncated norma! distribution in which the coefficient K_, is determined from the 
condition governining the normalization of the probability density curve. 


In order to achieve a graphical plot and to determine the numerical distribution parameters it is convenient 
to study the normalized distribution for z = 4 w r/o and to express the magnitude of the threshold voltage in terms 
of the coefficient 


In U, —Ina 


a=, 7) 

































by representing it in the form 


Un => ae*’. (8) 
As a result we can write 

_ (2+a)* 

2 9) 
p(z) = ——— (z>0), ( 
2n V( 

iin V ma V(a) 

. wey (10) 
V (a) = ve | e *dt. 


Figure 2 shows the family of distribution curves for the duration of the peaks at the output of an RC filter for 
a number of values of the parameter a. It is evident from the curves that as the cutoff voltage decreases, the 
probability distribution density for the duration of the peaks approaches a normal distribution with an average value 
close to -a. 


The average duration of the peak at the level corresponding to U,, is determined from the expression 


a? 





fi ape Ba (11) 
(= Tv@ 
whence we obtain 
Ee aes sia) 12 
T (a) 8. z es (12) 


by taking into account the relationship between z and r. 
The probability that the duration of a peak at the level Uy, will exceed the duration rT» is equal to 
V (zo + 13 
P(t >t) = Pe >n)=—-Fo. mm 


2. A Low-Pass Filter of the Gaussian Type 
“The frequency response for such a filter is written as 


Cja)=e ™ O<ecer). (14) 
The passband of the filter at the 0.7 level is equal to Aw = 0.84 y. 





For conditions identical to calls which were applied to the RC filter the pulse at the output of the filter is 
defined by the relationships 


—_—— (15) 
C. 
Uout (t) = Vix e* (wo <t<~). 





The time r during which the noise peak with an amplitude x exceeds the threshold voltage Uy, is found from 
the equation 


8 (stedes y (16) 
Uy=0.48CAwze © * * , 
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By solving this equation we obtain 
(0.42A@rt)? = Inz—InU,, (17) 


where Un = me ou is the value of the threshold voltage normalized to the filter input. 
. oh@ 


Thus the quantity (0.42 Awr)} has a distribution that coincides with the distribution for In x but with a shifted 


distribution center. Just as in the preceding case (8), it is convenient to express Uj, as 


U, = ae" (18) 
2 
and to study the normalized truncated normal distribution for v = feet : 
- (otse as 
p(v) = (v >0). 


" ¥ 2xV (8) 


In that case the distribution for the quantity ¥V=u can be written as 


_ weer (20) 


Figure 3 shows the probability distribution density curves for the duration of the peaks in the case of a low- 
pass Gaussian filter for a number of values of the parameter 8. The maximum value for p(u) is obtained when 


nn V $4 Rod 


here the minimal value for p (Upjgx, 8) corresponds to 6 = = ; 





The probability that the duration of a peak at the level Uy, will exceed the specified value T9 is equal to 
(just as in the case of the RC filter) 


V 
Pit >) = Pv > m) = Ga. a 


In order to determine the average duration of the peak at the level corresponding 6 we study the expression 


"B= sp | u“e u. 


Pa After performing certain transformations involving the use of the substitution uw’ = x, formula (23) reduces to 
form 


(24) 




























































The integral (24) can be expressed in terms of confluent hypergeometric functions (cf. for example, [3]) 


ax* _ htt ket 


(ae? deme ao7 r(*t*) (= i: oO 
0 


2 i ae oe 


oY Erb +ipn(dtt. 2, 2) 


where T is a gamma function, and ,F; is a confluent hypergeometric function specified in the form of curves [4]. 
Based on (25), we can write the following expression for u(B): 








(25) 


oO = ace? ‘(a(S ee F)—V2e(G)A(G. FF). ew 


The values of u(B) determined from formula (26) are given in Fig. 3. 
The average duration of a peak at the level corresponding to 8 is equal to 


+... Swe (27) 
_ u(B) Vs 
+ (B) = “0.42 Aw * 


For comparison of results Fig. 4 shows two pairs of distribution curves for the duration of the peaks at the 
outputs of RC and Gaussian filters with identical passbands at the level 0.7 for identical values of the threshold volt- 
age Up. It can be noted from formulas (8) and (18) that the cutoff coefficients a and 6 are related by the following 
simple expressions for identical Aw and U;: 





0.74 (28) 


In order to simplify the plotting of the curves in Fig. 4 the value of o is assumed equal to 0.74. For this case 


B=a+i. (29) 


The analysis which we have performed of the distribution of the duration of the peaks for pulse noise with a 
logarithmically normal amplitude distribution for two types of frequency responses can easily be extended to filters 
with different frequency responses. Moreover, it is possible in analogous fashion to find the probability distribution 
density for the duration of the peaks which directly follow a detector which isolates the pulse noise envelope. Under 
these conditions the distribution of the durations will be determined entirely by the shape of the frequency response 
of the wide-band input filter. 


The results obtained above make it possible to compute the actual noise immunity of remote control systems 
that operate over channels with pulse noise of this type when the limitations imposed on the pulse repetition fre - 
quency and the passband width of the filters are satisfied. 


The problem of finding the optimal frequency response for a low-pass filter so as to insure the optimum noise 
immunity of a remote control system is of especial interest. This problem shall be treated in a special paper. 
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SOME CASES OF SIMILARITY OF TRANSIENT PROCESSES 
IN SINGLE-LOOP NONLINEAR CONTROL SYSTEMS 
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Many nonlinear control systems can be represented in the form of single-loop systems that consist of two 
portions—a linear portion with the transfer function W(p) and a nonlinear element NE (Fig. 1). 


An analysis of such control systems shows that, in the case of nonlinearities of a certain given type, similar 
transient processes can occur in such systems and the results of a particular experiment can be generalized. 


We shall assume that the characteristics of the nonlinear element are described by the equation 
z= f(zx,z, 62), (1) 


where z is the nonlinear element output quantity, x is the controlled variable, and §x is the quantity determining 
coordinate of the characteristic point of the NE nonlinear characteristic. 
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If f(x, x, 6x) is a uniform first-power function of x, x, and §x, i.e., the following equations hold: 


f(a, 2, d2) = dxf (=, =.) @) 


oR ee 





. (3) 
f(z, 2,02) _4/ 2 = 
oo) =7(£. 4.4), | 
then, the general eugqation of the system | 
W (p) = = f (x, 2,2) a) | 


' 
can be written in the following form: | 


Whit (de> wet): 

































Such a transition from the variables x and x to the variables x / 6x and x/ §x makes it possible, in these cases, 
to consider a control system with a single value of the characteristic nonlinearity parameter (Fig. 2) instead of a 
single -loop control system with a nonlinearity having different values of the characteristic parameter §x. 


2-f(z,2 62) 




















Fig. 3, 


Consequently, in single-loop control systems consisting of a linear portion and a nonlinear element that is given 
by f(x, x, 6x), the separation of a group of similar transient processes is possible if the nonlinearity is described by a 
uniform first-power function of x, x, x. Here, the simi- 
larity criteria are the ratios of the controlled variable x and 
of its derivative x to the coordinates §x of the nonlinearity 
characteristic point. 


In particular, uniform functions are used to describe 
such nonlinearities as Coulomb friction, backlash, non- 
linearities with a relay characteristic, and also nonlinearities 
with an insensitivity zone. The characteristics of the above 
nonlinearities have two symmetrical characteristic points 





















































(Fig. 3). From among the indicated nonlinearities, only 
\¢,62’ Coulomb friction and backlash have characteristics that 
a rm Kote depend on x and are described by the equation 
: rd . ‘ 
oa ol XN z =f (z, sign z, dz). (6) 
Gr” 
+e, Oz" — or -6r’ 
-Gfr\— ~-—Ot +62" Nonlinearities with an insensitivity zone and with a 
relay characteristic are described by the equation 
SSE nrs & F§ t 
Fig. 4. z= f(z, dz). 


Under specific conditions, uniform equations can also 
be used for describing nonlinearities that have a relay charac 
teristic with an insensitivity zone as well as a characteristic with limitation or with limitation and an insensitivity 
zone. The sufficient uniformity condition for these dependences is the simultaneous proportional variation of all the 
characteristic points of a nonlinearity, which secures the invariability of the nonlinearity in switching over to the 
x / 6x variable. 


In the cases where x/6x and x/6x are used as similarity criteria, the transition to new variables, which are 
formed by dividing the controlled variable x and its derivative x by the value 6x of this quantity at the nonlinearity 
characteristic point, makes it possible to arrive at the following general conclusions that are valid for the conditions 
considered. 
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1. Within a group of similar transient processes and during equal time intervals, the values of the controlled 
variables x' and x” are expressed by the §x'/ 6x” ratio, while the time interval remains unchanged within the group 
of similar transient processes (Fig. 4). 


A consequence of this conclusion is the fact that the self-oscillation amplitude is proportional to the §x value 
and the fact that the oscillation period is independent of the §x value. 


2. The character of the sequence of limiting cycles and the character of the trajectory of the point in the 
space between two limiting cycles of the phase diagram are independent of the 6x value. Hence, it follows that 
the stability boundaries of the control system in the domain of the linear portion parameters are independent on the 
&x value and that they are determined only by the nonlinearity type. 


3. In analyzing transient processes in similar single-loop nonlinear control systems of the type considered 
here by means of phase diagrams, it is advisable to use, as the coordinates, the ratios of the controlled variable and 
of its derivative to the coordinate 6x of the nonlinearity characteristic point. 


The author extends his thanks to Prof. V. S. Vedrov for his advice, which made it possible to present this com- 
munication in the above form. 
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The present article emphasizes the fact that a large number of arithmetical operations must be 
performed in using digital computers (DC) for solving the problem of optimum smoothing - out 
and for effecting linear transformations of sequences of signals containing random disturbances. 

It is demonstrated that, in a number of cases, the necessary number of arithmetical oper- 
ations can be considerably reduced by using the recursion formula. 


In using digital computers in automatic control systems, it often becomes necessary to perform smoothing -out, 
differentiation, or any other linear operation on the sequence of values of a signal containing a disturbance. The 
desired transformation can be brought about in the following manner by using sliding summation. 


We shall assume that the signal x(t) values are read at equal time intervals Ty : x, = x(kT,), where k= ..., -1, 
0,1,2,3,... . By using a linear combination of the last N + 1 readings 
. qa) 
y= b> Wn7k—n 
n=0 


and for a suitable choice of weighting coefficients w,, we can approximate any linear transformation by the useful 
component of the signal. 


By using the criterion for the minimum of the root-mean-square error, we arrive at a system of linear equa- 
tions for determining the optimum values of the weighting coefficients w,: 
































N 


> R(jTo—nTo)w,=F (i) (i = 0,1, 2,..., N), 
n=0 . 


and also at some additional conditions [1-3]. 


In Eqs (2), F(j) is the linear combination of certain functions, and R(t) is the autocorrelation function of the 
steady-state random signal. We shall assume that the signal random components are such that the spectral function 


corresponding to the autocorrelation function R(r) is a rational-fraction function: 


1M= >) Rayo *= 40 





al B (fC) ° 
m I 

A= > a0, B= >} oe. 
ji=—m j=—l 


In this case, as was shown in [5], the sequence of weighting coefficients w, can be expressed by 


2m I—m-1 
“_= ? (n) + > Bo, (n) + > IC, Mi DS, N—s l, 
=] e=0 


where @,(n) is the complete system of solutions of the homogeneous difference equation 


m 


A(2)@(n)= >) a@(n+))=0, 
=m 


and 


1 for »=5, 
b= 4 0 for n=s. 


In Eq. (5) and below, z designates the unit shift operator: 


zz, = Frit: 


(2) 


(3) 


(4) 


(5) 





The y(n) function in Eq. (4) is a linear combination of certain functions. If the useful component of the signal 


is approximated by a polynomial with the power r, y(n) is a polynomial of the same power: 


9 (*) = Yo+ 18+ ...+ 74,2", 


(6) 


If the optimum values of the weighting coefficients wn are known, the sequence y, (k=...,1,2,3,...) 
can be calculated by using Eq. (1). However, for large N values, a large number of arithmetical operations (N + 1 
multiplications and N additions) must be performed at each calculation step. This would impose severe requirements 


on the DC speed. 


As will be shown below, the number of arithmetical operations necessary for calculating y;, can be reduced at 


the expense of a certain increase in the utilized volume of the operative memory (OM) in DC. 
Assume that ¢(n) satisfies the difference equation 
Pp 


E(2)9(n) = >) pint h=0 (#0, ¢, + 0). 


k=0 


(7) 



































Then, if 
p+m 
G(z)=A (z) & (z) = > g,2°. (8) 
(2) wane 
where 
Pp 
‘| 
i> >) a, 4%, (9) 
the —_ 
tion 
(a, = 0 for k > m and for k < -m), 
G(z)@(n)=0, G(z)®, (a) =0. (10) 
We shall apply the z P ~™ G(z) operator to both sides of Eq. (1); as a result of this, we shall have: 
(3) 
n 
z P-™G(z)y,= > w, G(s)?" 2, 
or viene 
p+m N p+m 
By substituting the summation index >) 5 eVk-+-e—p—m = >) > 85" n* k—n—p—m+s° Q1) 
=——m n=0 s=—m 
(4) 
we find n=s-+a, —m—pqaqNn+m, 
p+m N+m p+m 
} 8 Yk—p—m+s — >) 7 k—p—m—a p> be" s40° (12) 
(5) s=——™m a=>—m—p s&>—_ 
Let us denote 
p+m 
d.= >) be%ete (13) 
s=——™ 
and divide the region where a changes into five intervals and then calculate d,, for each of these intervals: | 
p+m : 
signal 1) —m—pqaqm—ti; d.= > 8ssi.0) | 
s=—a | 
I—m—1 
2) mgagl—t; dg= > Cy 
6 a 8° s—a’ 
(6) pad (14) 
3) l<agqNn—p—l; d,=0, 
J i—m—1 
+1 4) N—l—p+igaqn—p—m; d,= >) Dgy_e—o 
ments 3x 
N—a 
5) N—p—m+1qQagQN+ m; qd, = >) Ss%epa 
ed at s=—m 
The calculation of d, in the first and the fifth interval does norrequire an explanation. In the third interval, 
dq = 0 because of condition (10). 
(7) In the second and the fourth interval, dg is expressed in terms of the coefficients in front of the 6 -functions 


in Eq, (4). In calculating the sums, it must be borne in mind that g,= 0 fors < -m ands > p+m. Finally, we 
obtain: 





Pm ' i—t N-+m 


Bees FS Went lft | (15) 
o=-—m ax=—m—p a=-N—!—p+l 


If we divide all terins by g,, +p&m +p = 4m&p * 0) and transfer all terms with the exception of y,, to the right- 
hand side, we find: 








p4+m—1 g 1} d N+m d 
- > ’ 
fiat > fe ey Ye—p—m4s + 2 2 7 k—p—m—e + > 2. e—p—m—a 
eam OM e=—m—p &m+p anN—i—pyi mtr 


(16) 


Equation (16) provides the possibility of finding y,, by using a recurrent process suitable for programming. If 
Eq. (16) is used, the number of calculations does not depend on N, and, for large N values, it can be considerably 
smaller than the number of calculations where Eq. (1) is used. 


In the particular case where m = 1 = 0 and E(z) = (1-z)P, Eq. (16) is transformed into the equation which was 
obtained in a different manner in [4]. 


Example. We shall assume that the autocorrelation function can be represented in the form 


R (vt) = e~*!"! feos Br + ¢ sin B | tI], (17) 


while the “regular” signal can be expressed by a second-power polynomial. In this case, the spectral function is 
given by A (2) 
f (2) os Oo “Biz)’ 





where (18) 
06 == sh 2aT', (ay — a3), A (z) =2ch p—\z + +), 
f - 4 
B(t)= [2 chy T,—\z + +) [2 chy, — (: + —)|. 
_, , Gad + crt) 
eT aa (19) 
7 a = COs BT, — _ sinBT, 
yae—B. JoeB = sy. a 
We further find: 


m=1, 1=2. D(a) = e™, O(n) = e™, 
a.j=—i. a4,=2chp, a,=—1, 
(A) = Po + Git + Gon’, 
E (z) = (1 — 2), 


p= 3, €, = i, e; = —3, ¢, = 3, és = —1, 
g-1 = —1, 66 = 3+ 2chp, g, = —4—6chp, 
62=4+6chp, g.=—3—6chp, 4 =1, 
d_4 = ga, 
d_s = g3Wy + gar, 
d_s = {ay + £9, + geld, 
d_y = g1Wy + Bai + Bsa + gards, 
dy = Soy + Ei + Baa + Bows + gar, 
dq,= g-1Cp, 
dy_= «Dp, 
dy_.= BoM y_5 t+ BM y_. + S2y_, + Esp + 8-1%y_4> 
dus = Bo y_. t+ Biy_. t+ G20y + 8-1 y_5> 
dys = Gy, + 81% y + -1y_,, 
dy = 6-1My_, + Boy: 
dn41 = 8-1¥y- 
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Equation (16) assumes the following form: 


Vx = — [6-19 pg + oY p—g + Way + G2Yu—9 + SsYy—y] + dat, + Coty, + dist, + 
Hg My gt doty ig + dit y_y + dy_ytp_y + Iy_g%e—w—r + Uy—2%—n—a + Sy + (20) 
+ fp_nigt Eytp_p__ + ly 43% e—n—s 


The latter equation is suitable for programming. In calculations, N + 6 OM cells are required for storing the 
Xy> X; -N -5 Quantities and six cells are required for the y+; .5 quantities; in all, N + 12 OM cells are required. 


Moreover, 17 coefficients must be stored in the memory: g-1 + gs; d+di; dy_,+dy4,. For the perfor- 
mance of each calculation cycle, 17 multiplications and 16 additions are necessary, while in calculations where 
Eq. (1) is used, for instance, for N + 1 = 100, 100 multiplications and 99 additions must he performed. 
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The present article is concerned with a method for investigating transient response conditions 
in a magnetic amplifier with internal and external feedbacks, which is encompassed by flexi- 
ble feedback and which has an inductive load that is connected through a rectifier. 


The problem of transient response conditions in magnetic amplifiers without feedback for an inductive load 
that is connected through a rectifier bridge has been treated in sufficient detail in [1, 2]. These investigations were 
concerned with amplifiers using cores and rectifiers with idealized [1] as well as nonidealized [2] characteristics. 


Amplifiers with a more complicated circuit are very often used in automatic control systems. Therefore, it 
would be of definite interest to determine to what extent can the method proposed in [2] for investigating transient 
response conditions be extended to the more general case where a magnetic amplifier has several control windings, 
internal or external positive feedback, flexible feedback, and an inductive load that is connected through a rec- 
tifier bridge. We shall limit our consideration to the weak signal regime in the magnetic amplifier, i.e., to the 





regime of such signals for which the currents and voltages acting in the amplifier's working windings do not exceed 
the limits of the rectilinear portions of the magnetization characteristics under transient as well as steady -state con- 
ditions, in connection with which the amplifier parameters can be considered as constant. 


1. Transient Response Conditions in a Magnetic Amplifier with Positive External 
and Flexible Feedbacks and with an Inductive Load that is Connected through a 
Rectifier 

Consider a magnetic amplifier which has two control windings W,; and Wc2, to which control signals with 
different polarities are supplied, a flexible feedback winding Wc3, and a positive feedback winding Wg. The am- 
plifier circuit is shown in figure 1. Here, ry and Cy are the resistance and the capacitance in the flexible feed- 
back circuit, respectively, and r, and L, are the load resistance and wevereven respectively. The excitation win- 
ding of a dc motor represents the load. 








By using the assumptions made in [2, 3], we shall compose the equivalent amplifier circuit, which is shown 
in figure 2. In this circuit, in correspondence with the designations and terminology used in [2, 3], tcy, Tce2, Tes, and 
Ifh are the resistance of the Woy, Wc2, Wcs, and Wg, windings; r, and rm, are the forward and the reverse resistances 
a ; Ry is the amplifier's internal resistance; S’ = ma S, where 
S is the amplifier response <p $ ana p are cooper which take into account the effect of commutation in 
rectification; Lj =1,/p; B= Wg /W.; BR, =We/W.i Woy =Wop=Wei Aigy, Aig, Aipy, Ai, Ai, are 
the current increments in the corresponding amplifier circuits “with betpeie to the initial current values. 





of the rectifier arm; R,= 


Let us consider the amplifier transient response equations for the case where the control voltages change in 
steps. Considering that, in core magnetization, the increments of the voltage at the rectifier input terminals and 
of the emf of the ac windings are equal and have opposite signs, the expression for the increment in the magnetic 
induction constant component for the equivalent circuit under consideration can be written as 


W 
AB,= toc (Aig — Aig, + BA — Be Aig,) + 
oe 





y Tl Ai) + Woe 4+, —+* 
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where Q is the cross-sectional area of the two cores, / is the average length of the magnetic lines of force, 
H, = OBJOH, for B_=const, y=—4B/dB_ for H,=const,and w is the alternating current angular frequency. 
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On the basis of the equivalent circuit for the amplifier output circuit, we can write the following equation: 


‘p’ " : . ‘d : . , 2 ‘ dAi d B 
SRE (Aig — Ain, + BAi, — BL Al,,) = Rj (Ai) + Aig) + ry Aiy + Ly a + Wp Q y 1 (2) 


The equations for the control winding circuits will have the following form: 








dAB, 

Ava =o Ai Wl @) 
adAB,, 

Auc, = Pop Micg -f W.Q 7 (4) 


and for the flexible feedback circuit, we have: 


S’R, (Aigg — Ai,, + BAI] — BLA) — Ri, (Ai, + Ais) = 
(5) 
{ dAB, 
= ("y+ Pos) Aicg + z\ At.,dt + W,,Q oy, ae 


Equations (1) - (5) completely describe the transient response conditions in the magnetic amplifier under con- 
sideration. After some transformations of Eqs. (1) - (5), we obtain a single equation for the transient response condi- 
tions in the magnetic amplifier under consideration: 
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2. Transient Response Conditions in a Magnetic Amplifier with Internal Feedback and 
an Inductive Load Connected to a Rectifier 

At the present time, magnetic amplifiers with internal feedback are extensively used in automatic control sys- 
tems, Therefore, it is of definite interest to determine to what an extent can the method for investigating transient 
responses that has been proposed in [2] be extended to the more general case where a magnetic amplifier with non- 
idealized characteristics has intemal feedback, external voltage feedback, flexible feedback, and an inductive load 
that is connected through a rectifier bridge. 
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The circuit of such an amplifier is shown in figure3. 





The magnetic amplifier has voltage feedback, which SR’ E Dic ahh hind A 
is realized by means of the Wg, winding. The amplifier is =< "A, 44° 


‘fyb “AySi3) 





encompassed by flexible feedback, which comprises resistor 

rg, capacitor Cs, and the winding Wr of the magnetic am- 

plifier. Moreover, the magnetic amplifier also has a mag- tm nd 
netization winding, Wr), which is connected to an addi- -@ 








tional supply source. The excitation winding of a dc motor, ia 
which has resistance r, and inductance L,» serves as the {Sw 





amplifier load. 


By using the assumptions made in [2, 3], we arrive at | We 9 


a 











the amplifier equivalent circuit that is shown in figure 4. 








Here, the effect of the magnetization winding W,, on the 
transient response conditions in the amplifier is not taken Br - 








into account. In the circuit of figure 4, 


Fig. 4. 


1= H+ Far Mg = "E+ Tw Cg=F, Ly =1)'P, Bi= WW, Be= We/W., B= W/W. 










The expression for the increment in the magnetic induction constant component for the equivalent circuit 


under consideration can be written as 





HeoW, ° 4 Vv , Pe , di 
AB, =O" (Ai, + BiAi + Bei — Pedi) + a ( i Ai + 4 1 4w_e 


We shall now write the equations for the amplifier circuits: 


output circuit equation: 
S'R,& (Aig + BiAij + Bad: — BsAis) = 


’ , aA 
= Ri (Ai; + Bin + Ai) + Fd +41 +W_ aE 


control circuit equation: 


aAB, 
Aue = reAic + WoQ—s-?: 


equation for the flexible feedback circuit: 
S’R,E (Ai, + BA + “ — BsAis) — Rj (Ai + Ais + Ai) = 
=f git +7- Cc = | Aint + We zr ore 


equation for the voltage feedback circuit: 


SRE (Ai, + BiAi; + Bais — Boi’) — Rj (Ai, + Ain + Ais) = 


dB, 


adAB 





dt 


a). 


(7) 


(8) 


(9) 


(10) 


(11) 


The system of equations (7) - (11) completely describes the transient response conditions prevailing in the am- 
plifier under consideration. After performing 2 number of transformations with these equations, we obtain the tran- 


sient response equation: 
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The expressions for the coefficients in Eq. (12) are given in the table. 


We used Eqs. (6) and (12) for calculating the transient responses in actual amplifier specimens, after which 
we obtained oscillograms of the transient responses in the amplifiers under investigation for a check of the calcu- 
lation data. 


A comparison between the theoretical and experimental data confirmed the correctness of our calculations 
and conclusions. 


SUMMARY 
1. The performed calculations and experiments confirmed the possibility of using the method for investigating 
transient response conditions that has been proposed in [2] for the majority of magnetic amplifier circuits with in- 
ternal feedback. Regardless of the great amount of calculation involved in determining the coefficients of the equa - 
tions, this method is sufficiently simple and clear. 


2. The obtained equations for the calculation of transient response conditions in fairly complicated am- 
plifiers with feedback make it possible to estimate the effect of the parameters of amplifiers and rectifier bridges 
on the duration of transient response conditions in an inductive load for weak signals. 
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An All-Union conference on contactless magnetic elements which are used in automation and computer 
techniques was held on February 20-23, 1961 in Minsk. This conference was organized by the Department of Solid 
State Physics and Semiconductors, Academy of Sciences, Belorus, SSR, the Scientific and Technological Division 
of the Belorussian Republic, The State Scientificotechnical Committee of the Council of Ministers, Belorus, SSR, 
and the Commission on Automation Contactless Magnetic Elements, Academy of Sciences, USSR. The aim of this 
conference was a discussion of problems in the theory, design, application, and development of contactless mag- 
netic elements which are used in automation and computer techniques. 


Over 400 representatives from scientific research organizations, higher institutes of learning, and industrial 
enterprises participated in this conference. Nearly 40 reports and communications, which aroused great interest 
among the participants at the conference, were read and discussed. 


The president of the organizational committee of the conference, N. N. Sirota (Minsk), opened the conference 


The reports by B. S. Sot-skov and M. A. Rozenblat were devoted to a survey of the contemporary state of the 
technology of contactless magnetic elements. In his report, entitled "The scope of electromechanical magnetic 
and semiconductor elements,” B. S. Sot-skov (Moscow) discussed the main characteristics with respect to which 
elements of different types can be compared. These characteristics include the following: the input and output 
power, the actuation time, the input and output resistance, the multiplicity factor, the type of current, the number 
of circuits, the operating temperatures, the maximum allowable mechanical acceleration, reliability (the factor 
of breakdown danger), and the limits with respect to the number of operating cycles. 


This comparison indicated the natural division of the fields where these elements are used. Besides the 
enumerated characteristics, among the other important characteristics are also power per unit volume, the life- 
time before breakdown, the complexity of servicing, and the frequency with which a device breaks down. The 
technical and economic advantage in application depends on reliability, and finally electromechanical devices 
should be used in this order. 


In his report “Modern contactless magnetic elements," M, A. Rozenblat (Moscow) proposed a possible classifics- 
tion of magnetic elements with respect to their operating principle: 1) Parametric elements, where the nonlinear 
character or the magnetization curve is utilized (ferromagnetic stabilizers, frequency multipliers and dividers, 
ferroresonance devices, etc.); 2) magnetic-modulation elements, i.e., elements with the action of two fields that 
differ with respect to their frequency (magnetic modulators, probes, magnetic amplifiers, magnetic relays, etc.); 
3) hysteresis elements (memories, magnetic recording devices, counters, core integrators, etc.). 











This report also indicated the main lines of the development of magnetic elements: increase in operating 
speed, improvement of efficiency and reduction of over-all dimensions and costs, improvement of technological 
efficiency, and dispensation with rectifiers and transistors in discrete -action magnetic elements. 


The reports by N. N, Sirota and K. D. Mart'yanov were devoted to the present state of the development of 
magnetically soft materials and the physical processes occurring in them. In his report “Ferrites with a rectan- 
gular hysteresis loop and magnetic reversal processes in ferrites," N. N. Sirota presented the results obtained in 
investigating the dependence of ferrite properties on their composition. He also discussed the general laws govern- 
ing the variation of characteristics in binary and ternary hard alloys. Besides, the composition, temperature, and 
the magnetizing field strength, the shape of the hysteresis loop also depends on structural uniformity, density, and 
the presence of cavities. Additional annealing broadens the loop, however, thermomagnetic treatment in the 





















































nce 


si fict - 


it 











region of the Curie point and the addition of small components make it possible to vary the characteristics of 
ferrites in the desired direction. The author briefly dealt with the mechanism of magnetization reversal processes 
in connection with the switching time and the composition of ferrites. 


K. D. Mart'yanova (Moscow) reported on new developments in the field of soft ferromagnetic iron-nickel - 
cobalt alloys with a rectangular hysteresis loop and higher temperature stability (in the temperature range from 60 
to 300° C). 





The report by V. G. Baranovskii and I, A. Petrusenko (Leningrad) was devoted to the development of a push- 
pull magnetic bridge amplifier with positive feedback for even harmonics and with a common ac output. The opera- 
tion of the proposed circuit, its basic properties, and the conditions for securing even harmonics in windings en- 
compassing two cores were discussed. The results obtained in testing this circuit under steady-state and transient 
conditions when the temperature changed from -40 to +60°C were described. 





This circuit can be controlled by direct as well as alternating current. 


As an example of application, a servo system with a two-phase asynchronous motor, where a feedback cir- 
cuit with even harmonics is used as the output stage, was discussed. The application of this circuit made it possible 
to reduce the size and weight of the output stage. 


A number of reports were concerned with the results obtained in developing and investigating highly sensitive 
magnetic amplifiers and magnetic modulators. A. 1. Moskalev (Moscow) reported on the development of a push- 
pull de magnetic amplifier circuit with a single core, which made it possible to organize the industrial production 
of amplifiers with a zero stability of 1 x 107"w at 60°C. 





In their report "A two-phase magnetic modulator with a low sensitivity threshold,” F. 1, Kerbnikov and M. A. 





Rozenblat (Moscow) described the application of the higher frequencies of magnetic -transistor oscillators in developing 


magnetic modulaters with two-phase supply. The modulators were based on ferrite cores with mutually perpen- 
dicular excitation and signal fields. Two-phase feed modulators make it possible to dispense with the use of filters 
in the supply, control, and output circuits. For a supply frequency of 500 cps, an experimental model had a lower 
sensitivity threshold of 10-" w, which was reduced to the 1-cps passband. M. A. Rakov and L. A. Sinitskii (L'vov) 
described similar modulators where two-phase supply is provided by means of two magnetic frequency dividers. 


In their report "Magnetic pulse amplifiers," V. G. Anikin and 1, A. Ikorokov (Moscow) dealt with the problem 
of improving the sensitivity of magnetic amplifiers with self-saturation by varying the supply voltage shape, Due to 
pulsed supply, the sensitivity can be increased to 10-“-10-" w. Methods for designing pulse sources with magnetic- 
transistor elements were also considered, 








E. I. Gurvich and L. B. Shchukin (Moscow) described methods for measuring the characteristics of miniature 
memorizing transformers (with a diameter of 2-7 mm), which are to be used in contactless elements in automation 
and computer techniques. 





The report by G. D. Kozlov (Moscow) was devoted to a method for calculating the basic core characteristics 
(the basic magnetization curve and the limiting and the particular magnetization reversal cycles) by approximating 
hysteresis loops with parallelograms. 





The effect of the core shape on magnetization reversal processes was discussed in a communication presented 
by_K. E. Volkovitskii (Moscow). 





V. V. Sarv (Tallin) presented a report entitled “Acceleration of the growth of a magnetic set by means of 
accumulation and repeated utilization of the energy of a decreasing magnetic field." 


E. L, L'vov (Moscow) discussed a method for determining the transfer function of choke-coupled magnetic 
amplifiers with resistive -inductive de loads. 


The author of this report obtained an expression for the transfer function of an idealized magnetic amplifier 
in terms of its circuit parameters and the input-putput characteristic. The expressions for the transfer and the tran- 
sition functions were compared with the expressions obtained in other papers, and the parameter regions where the 
equations yield results different from experimental results were determined. 


V, V. Kokhanoy (Frunze) presented a communication on the development of a contactless automatic pulse 
generator which is designed for operation in combination with a data transmitter capable of converting nonelec- 
trical quantities into ac voltage. 














Many reports presented to the Conference were devoted to the results obtained in theoretical and experimental 
investigations and in the application of magnetic logical elements for industrial control systems. 


N. P. Vasil’'eva (Moscow) and I. Gaskovec (Prague) compared the operation of magnetic logical repeater cir- 
cuits under steady-state conditions and under maximum -output conditions. 





It was shown that a stable characteristic of circuits with decoupling by means of controllable diodes can be 
obtained by using a special supply voltage shape, while in circuits with decoupling by means of resistors, this can 
be achieved by a suitable choice of the circuit parameters. 


It was also shown that, for a frequency of 50 cps, the decisive factors in determining the maximum number of 
elements are the conditions of the limiting magnetization reversal and nonreversal of magnetization in cores, while 
in the case of higher supply frequencies, the main factor is the heating of the control element core winding. 


With respect to the above characteristics, the best circuit is that with decoupling by means of resistors. 


V. S. Podlipenskii (Kiev) presented a communication entitled "Some theoretical relationships for logical- 
action magnetic elements.” Elements based on the Reimi circuit were considered. 





L. V. Shopen (Moscow) presented a calculation of magnetic logical elements (on the basis of push-pull resistors 
and a transformer circuit) for the case where the power is supplied by a sinusoidal voltage source. 





I, A. Baranov and N. L. Prokhorov (Moscow) described the application of the magnetic logical elements that 
have been developed at the Institute of Automation and Remote Control, Academy of Sciences, USSR, in control 
circuits for controlling the electric and hydraulic drives of machine tool units and automatic assembly lines. 





The report was devoted to methods for designing units and control circuits by using magnetic logical elements, 


K. R, Lebedev (Kalinin) reported on logical element types whose mass-production is planned at the Kalinin 
plant for electrical equipment. The author described the supply sources and the control signal transmitters, and 
he provided information on the parameters of the elements as well their design. 





B. V. Serov (Petrozavodsk) presented a report on the application of magnetic elements in devices for the 
automated sorting of heavy timber. 
M. G. Kogan (Sverdlovsk) presented a report on control devices which are based on magnetic logical elements. 


The report by V. S. Matorina and O. A. Sedykh (Moscow), entitled “Lag element for systems with magnetic 
logical elements,” was devoted to various types of time lag elements, the operating principle of which is based on 
the dependence of the core magnetization reversal time on the voltage applied to the core winding. 





The report also dealt with the possibility of using integrating cores for designing time lag elements. 


M, M. Lotosh (Moscow) reported on the development of large smoothly controllable time lags which are based 
on magnetic cores. 





Circuits which make it possible to obtain lag times of up to several hours were considered. 


The author also provided a comparative analysis of circuits with low-frequency oscillators and of circuits 
with frequency dividers. Much attention was paid in the report to the development of a divider for dividing a num- 
ber of pulses by a large factor. 


A lag line consisting of a ferrite -triode cell was considered in the report by N. A. Yarmosh (Minsk). 





The report by A. L. Pisareva and G. N. Derim-Oglu was devoted to the development of output magnetic am- 
plifiers for contactless automation circuits. Contactors and ac electromagnets are most widely used as loads, which 
involves great difficulties in developing magnetic amplifiers due to sharp changes in the load in dependence on the 


magnetic system's position. 


A series of magnetic amplifiers for the control of ac contactors with a power consumption of up to 350 va as 
well as for the control of ac apparatuses (TU series for 50 cps and TUAS series for 400 cps) have been developed. 





The report “Transformation of digital codes into analog codes by means of magnetic amplifiers," which was 
presented by M. A. Rozenblat and G. V. Subbotina (Moscow), was devoted to decision devices consisting of mag- 
netic elements. This report was concerned with the theory and the results obtained in the development of a mul- 
tiplace converter. The report by V. I. Semenova and £. F. Stepura (Leningrad), entitled “Application of magnetic 
amplifiers for the multiplication of two or several electrical signals" was devoted to the same subject. 
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In their report "Magnetic amplifier with combined feedback,” V. G. Anikin and I, A. Okorokov (Moscow) de 
scribed the development of an amplifier with large power gain and with high stability of the input-output response. 





A. M. Bamdas, S. V. Shapiro, Yu. A. Savinovskii, and I. V. Blinov (Gor'kii reported on their work on logical 
frequency converters and on voltage stabilizers with magnetizable shunts, on the development of single -phase and 
three-phase voltage stabilizers with magnetizable autotransformers, and on new static frequency multipliers. 





In their report "On the theory of parametrons,” O, N. Sitnikov and G. Ya. Karasik (Sverdlovsk) analyzed the 
parametron operation by taking into account changes in the dynamic permeability of cores. 





The authors demonstrated that the detuning of the parametrons is determined by the nonlinearity oi the mag- 
netic permeability and that a definite operating frequency range exists for a given material. 


V. S. Toropov (Moscow) presented a report on a new method for information recording in magnetic memories. 








A. B. Gorodetskii (Moscow) presented a report entitled "A promising direction in the development of mag- 
netic amplifier technology,” in which he offered a comparative survey of the existing designs and the parameters 
of domestic and foreign makes of magnetic amplifiers; he also set forth a number of proposals concerning the ter- 
minology and the methods of graphical representation of magnetic amplifier circuits. 


V. M. Kurotchenko and N. I, Babanov (Frunze) reported on the development of the BKU-60 contactless code- 
pulse device for the telemetering of rotation angles. 





Stefanovich, Shvarts, Gorodetskii, Isaev (Moscow), Atlasov, Petrusenko (Leningrad), Il'in, Chaikovskii (Minsk), 











Chizhukhin (Penza), Pashkovskii (L'vov), and others participated in the discussion which followed the reading of the 


reports. It should be mentioned that great interest was shown in the discussion concerning the problems in the tech- 
nology and mass-production of standardized magnetic elements and units as well as the problems connected with 
perfecting the terminology. 


At the conclusion of the conference, the participants adopted a resolution emphasizing the prospects offered 
by the use of contactless magnetic elements and the necessity of further scientific research and experimental work 
in this field. The resolution also proposed a number of measures to be taken with the aim of accelerating the in- 
troduction of the most advanced contactless elements into the national economy. 
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